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Abstract 


Electronic  analogue  computers  and  additional  circuitry 
are  employed  to  simulate  two  interconnected  hydro-power 
plants  acting  under  the  influence  of  a  load- frequency  time 
controller 0  The  controller  which  was  developed  by  Calgary 
Power  Coo  Ltd,  is  unique,  to  the  author's  knowledge. 

The  units  in  the  individual  plants  are  assumed  to  act  as 
one  due  to  synchronizing  power  between  machines.  It  was 
further  assumed  that  the  generators  of  the  two  plants 
have  identical  characteristics, 

A  study  was  made  to  obtain  controller  settings  for 
optimum  plant  operation  for  different  hydro  and  governor 
characteristics , 


ACKNOWLEDGEMENTS 


The  author  wishes  to  express  his  sincere  thanks  to 
Professor  Y0  J,  Kingma  for  the  assistance  given  throughout 
the  project  and  to  Calgary  Power  Limited,  who  permitted 
the  analysis  of  their  controller.  Special  mention  is  also 
due  to  Professor  J,  B,  Nuttal  of  the  Civil  Engineering 
Department  and  G.  Gaskell  of  Calgary  Power  Limited  for 
their  helpful  suggestions  in  determining  the  hydro  and 
controller  transfer  functions  respectively,  A  sincere 
thanks  to  Mr,  J,  A,  Ash  for  his  assistance  in  the 
construction  of  the  analogue  simulation. 

Acknowledgement  of  financial  support  given  by  the 
University  of  Alberta  is  made. 


1.  «.  ‘ <'''  V  '  '  ' 


:  ‘  )qc  "  f>»  :.1  .r.  r  i  >  i  r.\  ;  >3.  >;0;  'vioA 


. 


GENERAL 


OF  E  Q  U  I  P  M 


VIEW 


ENT 


Hii 


iii 

TABLE  OF  CONTENTS  Page 

ABSTRACT  i 

ACKNOWLEDGEMENTS  i;j_ 

TABLE  OF  CONTENTS  iijL 

CHAPTER  I 

SECTION  I  — 1  General  Introduction  1 

Synopsis  5 

SECTION  1-2  The  Method  of  System  Simulation  6 

CHAPTER  II 

SECTION  II-l  The  Generator  and  Voltage  Regulator 

Section  II— 1-0  Introduction  12 

Section  II-l-l  The  Stator  Equivalent  Circuit  13 

as  seen  by  the  Stator  Phase 

Section  II-1-2  The  Stator  Equivalent  Circuit  18 

as  seen  by  the  Rotor 

Section  II-1-3  The  Comparison  of  the  Generator  21 

Equations  as  seen  from  the  Stator 
and  Rotor 

Section  II-1-4  The  Rotor  Per  Unit  Quantities  22 

Section  II-1-5  The  Representation  of  Generator  25 

Saturation 

Section  II-1-6  The  Voltage  Regulator  26 

Section  II-1-7  The  Generator  and  Voltage  27 

Regulator  Equations 

Section  II-1-8  Conclusion  29 

SECTION  II-2  The  Speed  Governoring  Model 

Section  II-2-0  Introduction  30 

Section  II-2-1  The  Power  Speed  Relationship  31 

Section  II-2-2  The  Governor  Transfer  Function  34 

Section  II-2-3  Conclusion  40 


' 


IV 


Page 


SECTION  II- 3 

The  Hydro  Model 

Section 

II-3-0 

Introduction 

40 

Section 

I I- 3-1 

The  Idealization  of  the  Circuit 

42 

Section 

II-3-2 

The  Idealization  of  the  Piezometric 
Head 

43 

Section 

II- 3-3 

The  Relation  Between  the  Pressure 
Variations ,  Change  in  Conduit  Area 
and  Fluid  Compressibility 

43 

Section 

II- 3-4 

The  Relation  Between  the  Flow  Rate 
Variation  and  the  Change  in  Volume 

45 

Section 

II-3-5 

The  Relation  Between  the  Pressure 
Variations  and  Water  Inertia 

46 

Section 

II-3-6 

The  Approximation  of  the  Equations 

47 

Section 

II-3-7 

The  Solution  and  Interpretation  of 
the  Equations 

49 

Section 

II-3-8 

The  Transient  Hydraulic  Performance 
Transmission  Line  Analogy 

57 

Section 

II-3-9 

The  Hydro  Equations  for  Systems 
with  a  Surge  Tank 

62 

Section 

II- 3-10 

The  Idealization  of  the  Surge  Tank 

67 

Section 

II-3-11 

The  Hydro  Equation  for  Systems 
with  a  Surge  Tank 

68 

Section 

II- 3-12 

The  Effect  of  Higher  Harmonics  in 
the  Pressure  Variations 

75 

Section 

II-3-13 

The  Effect  of  the  Turbine  Wicket 

Gate  Non-linearity 

76 

Section 

II-3-14 

Conclusion 

77 

SECTION  II-4 

The  Load  Model 

Section 

II- 4-0 

Introduction 

78 

Section 

II-4-1 

The  Generator  Representation 

78 

Section 

II-4-2 

The  Current  Equations 

81 

Section 

II-4- 3 

Conclusion 

86 

' 


v  >-  1*0 


. 


■ 


SECTION  II-5  The  Controller  Model 


Section  II-5-0  Introduction  86 

Section  II-5-1  The  Description  of  the  Controller  86 

Section  II-5-2  The  Approximations  Made  of  the  9  4 

Controllers  Operation 

Section  II-5-3  The  Simulation  of  the  Controller  96 

Section  II-5-4  Conclusion  99 

CHAPTER  III 

SECTION  III-l  The  Analysis  of  Results 

Section  III-1-0  Introduction  100 

Section  III-l-l  The  Sustained  Synchronizing  101 

Torque  and  System  Frequency 
Osci llations 

Section  III-1-2  The  Time  Error  Resulting  front  a  112 

Load  Change 

Section  III-1-3  The  Optimum  Controller  Parameters  123 

for  a  10  Second  Sampling  Duration 

Section  I I I- 1-4  The  Optimum  Sampling  Duration  135 

Section  111-1=5  The  Effect  of  Dead  Zone  and  142 

Backlash 

Section  III-1-6  Conclusion  152 

SECTION  III-2  The  Computer  Results  155 

SECTION  III-3  General  Conclusion  169 

REFERENCES  174 

APPENDIX  A 

SECTION  1-2  The  Symbols  Utilized  in  the  Method  176 

of  System  Simulation 

SECTION  II-l  The  Symbols  Utilized  in  the  176 

Generator  and  Voltage  Regulator 
Models 

SECTION  II-2  The  Symbols  Utilizing  the  179 

Governoring  Model 


' 


S-  ■  T  .{  'ApilT'J.  ' 


■ 


1  o  ■  '  ■' 


VI 


SECTION 

II-3 

The  Symbols  Utilized  in 
Hydro  Model 

the 

180 

SECTION 

II-  4 

The  Symbols  Utilized  in 
Model 

the  Load 

184 

SECTION 

II-5 

The  Symbols  Utilized  in 
Controller  Model 

the 

185 

SECTION  III-l 

APPENDIX  B 

The  Symbols  Utilized  in 
Analysis  of  Results 

the 

186 

SECTION 

II-l 

The  Generator  and  Voltage 

Regulator  Parameters 

189 

SECTION 

1 1-2 

The  Governor  Parameters 

189 

SECTION 

II-3 

The  Hydro  Parameters 

190 

SECTION 

II-  4 

The  Load  Parameters 

19  0 

SECTION 

II-5 

The  Controller  Parameters 

190 

APPENDIX  C 

THE  COMPUTER  SCHEMATIC 


191 


CHAPTER  I 


Section  I  - 1  Introduction 

The  two  major  aspects  of  power  system  regulation  are 
frequency  and  load  control.  Due  to  the  advent  and  popularity 
of  synchronous  clocks,  it  has  become  of  increasing  importance 
to  supply  correct  time  to  the  consumer.  The  system  frequency 
is  adjusted  by  raising  or  lowering  the  permanent  speed  droop 
curve  as  required.  The  load  apportionment  between  units  is 
0 r^" q^inal ly  set  by  the  slope  of  the  above  curve,  however,  it 
may  be  altered  by  applying  a  bias  signal  to  frequency  control 
function.  Correct  time  can  also  be  incorporated  by  temporary 
over- regulation  of  the  frequency. 

The  frequency  on  modern  power  systems  is  held  within 
0.1  cycles  per  second.  This  exceeds  consumer  requirements 
in  regard  to  frequency.  As  already  indicated,  this  in 
itself  does  not  provide  for  accuracy  of  synchronous  clocks 
as  their  time  is  the  integral  of  frequency. 

The  control  unit  devised  by  Calgary  Power  Limited 
is  capable  of  providing  frequency  tolerances  considerably 
better  than  0.1  cycles  per  second,  load  sharing  between 
units,  and  correct  time  to  within  5/8  of  a  second,  based 
on  the  assumption  that  rapid  system  load  changes  do  not 
occur.  Synchronous  time  rather  than  frequency  is  the 
basis  of  operation.  If  a  time  error  exists,  the  governor's 
permanent  speed  droop  is  adjusted  to  correct  for  the 
deviation  in  frequency  and  a  second  signal  is  applied  to 
cause  the  unit  to  overspeed  to  correct  the  time  error. 
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Load  control  is  provided  by  applying  a  bias  signal  to  the 
controller  such  that  every  generator  may  operate  on  a 
different  time  datum. 

It  is  the  purpose  of  this  analysis  to  determine  optimum 
controller  settings  for  system  operation,  with  respect  to 
frequency ,  time,  load  distribution  and  power  disturbances. 

To  accomplish  this  it  was  necessary  to  simulate  a  system. 
Therefore  a  method  of  simulation  was  first  decided  upon, 
followed  by  the  development  or  acceptance  of  the  required 
building  blocks:  generator,  voltage  regulator,  hydro- 
system,  load  and  controller.  The  system  was  then  tested  by 
applying  consecutive  load  changes  of  25%  and  15%, 

To  facilitate  the  interpretation  of  the  results,  the 
controller0 s  action  on  the  governor  is  converted  to  equivalent 
frequency  error  signals.  To  generalize  the  results,  all 
simulated  quantities  except  time  and  synchronizing  torque 
angle  are  expressed  in  per  unit  quantities.  However,  in 
the  development  of  the  power  system  building  blocks,  it  is 
sometimes  necessary  to  represent  the  quantities  in  their 
respective  units.  Thus  to  avoid  confusion  in  the  forth¬ 
coming  analysis,  all  quantities  will  be  assumed  to  be 
expressed  in  their  units  unless  it  is  specifically  indicated 
that  they  are  expressed  in  per  unit. 

The  following  tests  were  made  on  the  system. 
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The  quantities  represented  by  the  symbols  are  defined  as  follows: 
8p  the  permanent  speed  droop  setting  in  per  cent 

k  the  effective  speed  droop  signal  expressed  as  a  per  cent 

of  the  permanent  speed  droop  settinq 

*• 

V  the  dashpot  frequency  signal  expressed  as  a  per  cent 

of  the  effective  speed  droop  signal 
^  the  percentage  load  setting  of  the  controller 

rw  the  water  starting  time  constant 

m  the  sampling  duration  in  seconds 


The  bracketed  tests  are  shown  in  Section  III-2  of  Chapter  III, 
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Synopsis 

The  required  system  building  blocks;  generator,  governor, 
hydro-system,  load  and  controller  models  were  simulated  on  the 
computer^,  The  voltage  regulator  model,  developed  in  Reference  4, 
was  accepted  and  incorporated  with  the  generator  in  Section  II- 
1-6  of  Chapter  II e 

The  operation  of  the  controller  was  not  found  to  adversely 
affect  the  voltage  variations  in  the  generator  and  voltage 
regulator  models  for  the  sampling  periods  tested  (2  to  10 
seconds)  „  Further  it  was  found  that  the  controllers  did  not 
affect  the  power  oscillations  between  the  generating  units. 
However,  it  was  found  that  the  sampling  periods,  of  the 
individual  generators ,  should  be  equal  to  or  greater  than 
the  natural  period  of  oscillation  of  the  governoring  loop 
(10  seconds)  (20).  The  governoring  loop  consists  of  the 
generator,  turbine  and  the  governor^ 

It  was  found  that  the  controllers0  signals,  between 
generating  units,  should  be  related  through  the  permanent 
speed  droop  curves 0  Further,  it  was  found  that  the  magnitude 
of  the  individual  controller0 s  signals,  were  dictated  by 
the  requirements  for;  a  minimum  system  disturbance,  maximum 
controller  response  and  accurate  load  division  among  the 
generating  stations. 
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Section  I  -  2  The  Method  of  System  Simulation 

Two  methods  of  system  simulation  are  considered. 

First*  the  representation  of  the  generator ,  voltage  regulator* 
governor  and  hydro-system  on  the  analogue  computer  with  the 
load  being  represented  on  the  network  analyzer.  Second*  the 
complete  system  is  represented  on  the  analogue  computer. 

Method  I 

The  network  analyzer  operates  on  a  frequency  of  500  cycles 
per  second.  It  is  mainly  used  for  transient  stability  studies 
where  constant  voltage  behind  the  direct  axis  synchronous 
reactance  and  constant  input  power  over  the  duration  of  the 
disturbance  can  be  assumed.  As  the  desired  studies  are  of 
the  steady  state  nature*  it  is  necessary  to  automate  the 
analyzer  as  the  former  assumptions  are  no  longer  valid. 
Therefore  two  possibilities  exists  operate  the  system  at 
500  cycles  or  represent  the  steady  state  conditions  by  D.C0 
and  convert  between  A 0C3  and  D.Ch  when  going  from  the  computer 
to  the  analyzer. 

A. Co  -  Ao Co  representation  (refer  to  Figure  I-2-I) 

The  unit's  direct  and  quadrature  axes  currents  and 
voltages  are  represented  by  A.C.  quantities.  The  power 
output  is  a  DoC.  signal  directly  proportional  to  VqIg+  vdId. 
This  is  developed  as.  follows: 


- 
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MULTIPLIER 


A. Co  -  A. Co  REPRESENTATION 
Figure  1-2-1 
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Unit  #1  Unit  #2 

Figure  1-2-2 


The  power  output  of  the  machines  may  be  obtained  as  a 
DoC.  signal,  from  the  product  of  the  A0C0  currents  and 
voltages,  by  utilizing  a  low  pass  filter.  This  may  be 

realized  from  the  following  mathematical  statements  (refer 
to  Figure  1-2-2) „ 

Unite  #1 

pi  +  Harmonics  =  2 (V  cos  wt  ^  cos  (wt  -  ©1) 

+  Vdl  sin  wt  I1  cos  (wt  “  ep)  ) 

^  cos  2wt)  cos  0^  +  sin  2wt  sin  0^) ]] 

+  vdl  cos  2wt)  sin  +  sin  2wt  cos  9^)] 


(1-2-1) 
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Unite  #2 

P2  +  Harmonics  =  2  [Vg2  i2  cos (wt+S12) cos (wt-©2+&12) 

+  Vd2  I2  sin (wt+$!2) cos (wt”@2+^12^1 
=  Vg2  I 2  ! l+cos2 (wt+5^2) cos92+sin2 (wt+S^2) sin02 ] 

+  Vd2  I2  t1+cos2 (wt+ £  ) sine2+sin2 (wt+Si2) cos92] 

(1-2-2) 

Utilizing  a  low  pass  filter  one  can  eliminate  the  1000 
cycle  components ,  producing  equations  1-2-3  and  1-2-4,  which 
are  equal  to  the  developed  electrical  power  at  the  air  gap 
when  the  equations  are  expressed  in  per  unit  quantities. 


P1  =  Vql 

h 

cos©! 

+  vdi 

I 

1 

sin  ©! 

(1-2-3) 

P2  "  Vq2 

*2 

cos02 

CN 

> 

+ 

X2 

sin  ©o 

(1-2-4) 

Thus  frequency g  torque  angle g  input  and  output  power 
are  represented  by  D0C0  quantities.  It  is  necessary  to  phase 
modulate  the  voltage  of  either  unit  with  the  D.C.  voltage 
depicting  the  synchronizing  torque.  The  two  major  problems 
encountered  in  the  above  simulation  are  the  phase  modulation 
and  the  representation  of  generator  saturation. 

A.C.  -  D.Co  Representation  (refer  to  Figure  II-1-3) . 

The  portrayal  of  saturation  may  be  included  in  the 
construction  of  a  D.C0  generator  model.  This  approach  also 
requires  phase  modulation  and  provides  further  problems  as 
it  is  necessary  to  make  A.C.  to  D.C.  and  D.C.  to  A.C. 
conversions . 

Although  no  evidence  of  A.C.  -  A.C.  simulations  have 


S6V  ♦ 


. 
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A®  C 


Do  C®  REPRESENTATION 

Figure  1-2-3 
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been  found  in  literature „  references  1,  2  and  3  have 
successively  attempted  A.C.  to  D.C.  conversions „  Corres¬ 
pondence  with  the  authors  of  references  2  and  3  was  made» 

Dr»  James  E0  Van  Ness,  reference  2,  was  sceptical  with  regard 
to  the  accuracy  of  his  method  and  recommended  Dr.  Kimbark8s 
approach.  Dr.  Kimbark  stated  that  high  accuracy  had  been 
attained,  however  the  circuitry  had  to  be  very  carefully 
balanced,  thus  tests  of  long  duration  could  produce  errors. 

Neither  of  the  above  references  were  able  to  supply 
diagrams  as  to  how  the  equipment  was  built.  Reference  3 
shows  the  required  circuitry  for  both  approaches  when  con¬ 
verting  from  D.C.  to  sixteen  kilocycles.  It  was  found  to 
be  complex.  In  view  of  the  former,  it  was  decided  to 
simulate  the  complete  system  on  the  analogue  computer. 

Method  II 

The  complete  system  is  simulated  on  the  analogue 
computer.  The  generator  may  be  placed  on  an  A.C.  or  a 
D.C.  basic.  The  former  encountered  the  problems  of  generator 
saturation  and  phase  modulation.  As  the  A.C.  representation 
offers  no  advantage  over  the  D.C.  representation  where  the 
steady  state  values  are  depicted  as  D.C.  voltages,  the 
latter  method  of  simulation  is  utilized  in  this  study. 
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CHAPTER  II 

Section  II-l  The  Generator  and  Voltage  Regulator  Model 
Section  I I - 1 - 0  Introduction 

The  stator  of  the  generator ^  under  balanced  three  phase 
conditions  can  be  portrayed  on  a  per-phase  basis  by  the  use 
of  Park 8 s  Equations „  The  rotor  operates  on  a  D.C.  basis  and 
the  equations  relating  field  voltage  and  current  are  known. 

The  relation  between  the  A.C„  per-phase  voltages  and  the 
D.C.  voltage ,  applied  to  the  rotor,  are  determined  by  the 
unit 8  s  voltage  regulator 0  The  voltage  regulator  equation, 
as  derived  in  reference  4,  was  accepted  as  correct  and 
utilized  in  this  study.  However,  to  obtain  a  simulation,  it 
was  still  necessary  to  determine  the  rotor  equations  in 
terms  of  the  stator  equations. 

The  derivation  of  the  equivalent  per-phase  stator 
circuit  as  seen  by  the  stator  and  rotor  is  outlined  in 
Sections  I I- 1-1  and  II-1-2.  The  relationship  of  the 
equivalent  circuit  as  seen  by  both  the  stator  and  rotor  are 
then  developed  in  Section  II-1-3. 

Motor  action  was  assumed  to  be  positive  in  the  develop¬ 
ment  of  the  equivalent  circuits  as  the  field  and  direct  axis 
currents  are  in  the  same  direction.  This  makes  it  unnecessary 
to  attach  a  vector  notation  to  them  and  simplifies  the  develop¬ 
ment.  This  assumption  is  accounted  for  in  writing  the  equations 
for  simulation  in  Section  II-1-7. 

The  rotor  parameters  are  transferred  to  the  stator  in 
Section  II-1-4.  The  effect  of  generator  saturation  is 


•  • 
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considered  in  Section  II~l-5„  The  voltage  regulator  equation 
is  stated  in  Section  I 1-1-6 0  The  equations  representing  the 
generator  and  voltage  regulator  are  given  in  Section  II-1-7* 

Section  II-l-l  The  Stator  Equivalent  Circuit  as  Seen  by 

the  Stator  Phase 

The  equations  determining  the  stator  and  rotor  self  and 
mutual  inductances  are  developed  for  phase  A0  (5) 


Axis  of 
Phase  b 


ic  =  I  cos (wt+120-©) 

Figure  II-l-l 

Stator  Self  Inductances  (refer  to  Figure  II-l-l) 

Stator  magnetomotive  force  F_  =  N~i  (II-l-l) 

The  force  along  the  direct  axis  Fda  =  Facos0  (II-1-2) 

The  force  along  the  quadrature  axis  Fga  =  Facos (0  -  90)  (II-1-3) 

The  flux  created  by  the  stator  m0mc,f0 


■ 


' 
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0gda  =  FdaPgd  =  FaPgdcos^  along  the  direct  axis  (II-1-4) 

0rT__  =  F  P  =  -F  P  sin0  along  the  quadrature  axis  (II-1-5) 

gqa  qa  gq  a  gq  ^ 

Flux  linkages  of  phase  A  created  by  the  air  gap  flux 

0gaa  =  ^gdacos^  "  0gqasin0  (11-1-6) 


D 


=  Fa  I  P^cos20  +  P^sin  0 


gd' 


gqk 


(1 1-1-7) 


Naia 


P_.„  +  P  -]  —  P 

gq  gd  gd  gq 

+  cos  20 

- T -  ~ — T - 


(XI-1-8) 


The  self  inductance  of  phase  A  due  to  the  air  gap  flux 


L  =  N.0 
gaa  a 


gaa 


r 

L 


N,^  jP^+P  +  P^-P 
a  *  gd  gq  gd 


gg  cos  20 


=  Lgo  +  Lg2  cos  20  (II-1-9) 

The  self  inductance  of  phase  A 

L  =  Lai  +  L  =  L  +  L  +  L  0  cos  20 
aa  ai  gaa  go  al  g2 


=  Laao  +  Lg2  cos  20  (II-1-10) 

Stator  mutual  inductances  between  phases 
The  flux  created  by  phase  A 


da’gd  ^gda 

along  the  direct  axis 

( *1-1-4) 

'  P  =0 

qa  gq  ^gqa 

along  the  quadrature  axis 

(I I- 1-5) 

-  bpW  -  sbe** 


51  «i8pp<l  ♦  .* 


y*j  :  qsg  ixb  off*  o:t  oi/b  A  sped  lo  vDcwloubnl  Jlss,  9t1T 


y  £  JO  *■  r.  ,  J 
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The  flux  created  by  phase  A  linking  phase  B 


0 


gba  =  0gda  cos  ~  120>  -  0  gqa  sin<0  “  120> 


(II-l-ll) 


0gba  =  Fa  (pgd  cos0  cos(0  -  120)  +  Pgcj  sin0  sin(0  -  120)) 


=  F 


-P  ,  + 
gd  gq 


gd  gq 


cos (20  -  120) 


(II-l- 


The  mutual  inductance  between  phase  A  and  phase  B 


Lgba  =  Na  0gba  =  -°-5Lgo  +  \2  =°s(20  -  120)  (II- 1-13) 


1 

a 


The  mutual  inductance  of  the  rotor  with  phase  A 

Let  the  mutual  inductance  of  the  direct  axis  of  the 
stator  with  the  direct  axis  of  the  rotor  to  be  Laff  then 
the  mutual  inductance  with  phase  A  is  La^cos0o  The  rotor 
self  inductance  is  then  equal  to  where  Lff  =  +  La£. 

The  remaining  phase  inductances  may  be  determined  in 
the  same  manner.  The  complete  set  of  inductances  are  listed 
below: 

Laa  =  Laao  +  Lg2cos20 

^b  =  Laao  +  Lg2cos<20  +  120) 


cc 


aff 


Jbff 


cff 


aao 


f  f  a 


'ffb 


Jffc 


+  L  ~ cos (20  -  120) 

g2 

=  Lafcos0 
=  Lafcos(0  -  120) 
=  La^cos(0  +  120) 


12) 
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Lab  “  Lba  “  -°’5Lgo  +  Lg2cos(20  -  120) 

Lbc  =  Lcb  =  •°-5Lgc  +  Lg2oos(20, 

Lac  =  Lca  =  -°*5Lgo  +  Lg2COS(20  +  120) 

The  per  phase  equations  may  be  greatly  simplified  by  trans¬ 
formation  to  the  direct  and  quadrature  axis.  Thus  equations 
■^“^-—14,  II-1-15,  II-1-16,  II  — 1-17 ,  and  II— 1-18  are  defined* 


Fds  =  Na  l 

ia  cos0  +  ib  cos (0  -  120)  + 

ic  cos (0  +  120) ] 

(II-1-14) 

Fqs  =  Na  [ 

-ia  sin0  -  ib  sin (0  -  120)  - 

-  i  sin  (0  +  120)] 

(II-1-15) 

Jd  -  Ka  [ 

ia  cos0  +  ib  cos (0  -  120)  + 

ic  cos (0  +  120) ] 

(II-1-16) 

I  =  K  f- 
q  q  l 

i  sin0  -  ifa  sin (0  -  120)  - 

ic  sin(0  +  120)] 

( II-1-17) 

Fds  =  »«Id 

F  =  N  I 

qs  a  q 

(II-1-18) 

Kd 

K 

q 

Letting  —  2/3,  it  is  found  that  the  following  inverse 

transformation  produces  correct  per  phase  currents* 


V 

COS0 

-sin0 

ib 

— 

cos (0  -  120) 

-sin (0  -  120) 

1 

• 

l 

c 

cos (0  +  120) 

-sin(0  +  120) 

1 

=  i  =  A  I 

P  q 


Quadrature  axis  current  transformation 

cos0  cos (0  -  120)  cos (0  +  120) 

-sin0  -sin (0  -  120)  -sin(0  +  120) 


I 

.  r 

* 

d 

I 

q 

=  2/3 

i 

L  °J 

■A 


’A 


V- 


• 

i 

a 

• 

■'•c 

*  ■ 

=  I  =  B  i 

q  p 


■ 


17 


Per  phase  flux  linkages,  where  Ifg  and  Xfg  are  the  rotor 
currents  and  flux  linkages  as  seen  by  the  stator. 


Quadrature  axis  flux  linkage 


^d 

~^a 

B 

V 

J 

j 

i 

L  ^  J 

_V 

1*  ^"1 
\ 

a 

1 

f 

Laa 

Lab 

Lac 

Laf  f 

f  "s 

^a 

Lab 

^b 

Lbff 

% 

Lac 

Lbc 

Lcc 

Lcff 

Lc 

>'fs 

> 

■ 

Laff 

Lbf  f 

Lcff 

Lff 

ifs 

<  J 

=  \  - c  S> 


\  =  B^P 


Substituting  the  matrices  Tv 


B  >  and  i  =  A  I  into 

r  r  4 


r 

the  matrices  Tv  =  C  i  and  expanding,  the  following  equations 

lr 


are  obtained s 
”?v ,  =  t  L  _  +  I 


d  LLal 


fs  af 

fsLaf  ’  ‘d  I  "al 


+  3/,  (L„„  +  L„?)] 


=  IfAf  +  la  |L  =  ,  +  3/,  N  =  2P  "! 


2  '^go  T  “g2‘ 
2  Na  Pgdj 


“  XfsLaf  +  Jd 


Lal  +  Lad 


J 


(II-1- 


A  = 


+  3/2  -  Lg2>] 

=  Jq  [Lal  +  3/2  (Na2Bgq)] 


=  I. 


'al 


+  L  1 

aqj 


19) 


(II- 1-20) 
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[' 
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(II  —  1—  21) 


(II-1-22) 


Exnanding  determinant  I 


q 


B  i  ,  it  is  found 
P 


1 3  =  Iasin© 


(II-1-23) 


I  =  I  cos© 
q  a 


(II-1-24) 


zo  =  0 

Thus  it  is  seen  that  the  stator-rotor  mutual  impedances  as 


viewed  from  the  stator  phase  are  equal  to  V2  NaPgd  and 
^ / 2  NaPgq  alon9  t^le  direct  and  quadrature  axis  respectively 


This,  however,  does  not  establish  a  relationship  between 
La£  and  La£.  To  obtain  this  relationship,  the  direct  and 
quadrature  axis  circuit  as  seen  by  the  rotor  is  developed 
in  Section  II-1-2. 

Section  II- 1-2  Stator  equivalent  circuit  as  seen  by  the  rotor 

The  phasor  currents  and  windings  are  120  electrical 
degrees  apart.  This  results  in  an  m.m.f.  travelling  around 
the  stator  at  synchronous  speed.  Under  steady  state  con¬ 
ditions,  the  stator's  m.m.f,  appears  stationary  and  of 
constant  magnitude  to  the  rotor  as  the  latter  also  rotates 
at  sychronous  speed.  This  may  be  realized  by  referring  to 
Figures  II-l-l  and  II-1-2. 


V 


■ 
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Rotor  Quadrature  Axis 


The  stator  current  components,  referred  to  the  stator, 

along  the  direct  and  quadrature  rotor  axes  are  i  cos0,  -i  sin0, 

a  a 

i^cos (0  -  120°),  -ibsin(0  “  120°),  iccos(0  +  120°)  and 
-x^sin h*  x20^y  respectively* 


Expressing  the  relation  in  matrix  form: 


r  i 

cos0  cos (0  -  120)  cos (0  +  120) 

*a 

= 

-sin0  -sin(0  -  120)  -sin(0  +  120) 

*'0 

L*  J 

V2  v2  v2 

•  ^ 

Expanding  the  determinant  I  =  D  i 

~1  ir 


I'd  =  3/2  Ia  sip©  =  3/2  1^ 
I’q  =  3/ 2  la  =  3/2  Iq 


D  i 


P 


( II-1-25) 
(II-1-26) 


Direct  axis  m.m.f*  Fds  =  3/2  NaId 

Quadrature  axis  m.m.f.  F  =  N_I_ 

qs  ^  a  q 


(II-1-17) 

(II-1-18) 


■  »  m  .  ,  mb  ^  #  >  ' 
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Flux  created  by  the  stator  m0m.f» 


0  =  F,  P  , 

gd  ds  gd 

(II-1-27) 

0  =  F  P 

gq  qs  gq 

(II-1-28) 

The  mutual  inductances  as  seen  by  the  rotor 

are  ; 

Lgd  =  Nf^gd  -  NfNa  3/2  Xd  pgd 

If  If 

(II- 1-29) 

Lgq  =  Vgq  =  NfNa  3/2  Iq  Pgq 

1  f  If 

(II-1-30) 

As  the  stator  and  rotor  are  magnetically  coupled,  then  a 
relationship  1^  =  ^ / 2  N&  must  exist.  Thus  equations 

II-1-29  and  II-1-30  becomes 

Lgd  =  Nf2  pgd  =  <Nf/Na> 2  2/3  Lad 

(II-1-31) 

Lgq  ~  Nf  pgq  =  (Nf/Na)  ^  / 3  Laq 

(II-1-32) 

The  equations  as  viewed  from  the  rotor 

^dr  =  Tf  Lgd  +  Na  (3/2  Id)Lgd  (II-1-33) 


^qr  =  Na  <3/2  Iq)Lgq  (II-1-34) 


Xf  "  Jf(Lfl  +  Lgd)  +  Na  (3/2  Id)Lgd  (II-1-35) 


Substituting  equations  II-1-31  and  II-1-32  into  II-1-33 


- 
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to  11-1-35 ,  equations  II- 1-36  to  II-1-jI8  are  obtained 


9\ 


dr  = 


N. 


N. 


/s  IJad 


N. 


N. 


Xd  Lad 


(II-1-36) 


Tv 

o 


0 


(II-1-27) 


Lfl  + 

2 

2/ 3  Lad 

Nf 

+  _  ^d  Lad 

(II-1-38) 

.V 

J 

Na 

Section  II-1-3  The  comparison  of  the  generator's  equations 

as  seen  from  the  stator  phase  and  rotor 
The  rotor  flux  linkage  equations  will  first  be  transferred 
to  the  stator,  via  the  turns  ratio  (Na/Nf) ,  then  the  stator 
equations  will  be  re-written  to  facilitate  the  comparison. 

Rotor  Flux  Linkage  Equations  Transferred  to  the  Stator 


X 


d 


If  Lad  +  Xd  Lad 


(II-1-39) 


Tv  =  X  =  I  LaCT  (II-1-40) 

q  a  qr  q  aq 

%o  =  0 


Xfs  =  Na  ’'f  -  Nf 

2/ 3  If 

_3/ 2 

*Na~ 

2  Lfl+Lad 

iq  fq 

— 

(II-1-41) 


■ 
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Stator  Flux  Linkage  Equations 
\m  xfs  Laf  +  Id-‘Lal  +  Lad> 


=  Xq(Lal  +  Laq> 


■K0  =  o 


=  1  f s  Lff  + 


3  /  J  T 

4  2  Id  Laf 


Comparing  the  equations ,  it  may  be  seen  that  3/2  L  ^ 
Lff  “  Lad  +  3/2  (Na^f*2  Lfl 


(II-1-19) 

(IX-1-20) 


(II-1-22) 


Xfs  =  <VNa>  2/3  If 

Further,  it  may  be  realized  that  Park's  equations  divide 
the  stator's  direct  and  quadrature  currents  and  multiply  the 
impedance  by  a  factor  of  ^ / 2*  As  it  is  desired  to  represent 
the  generator  on  a  per  phase  basis,  this  factor  will  be 
applied  to  the  field  quantities* 

Since  the  stator  quantities  will  be  represented  on  a 
per-phase  basis,  it  is  necessary  to  establish  the  stator- 
rotor  turns  ratio  to  determine  the  rotor  base  quantities. 


Section  I I- 1-4  The  Rotor  Per  Unit  Quantities 

The  rotor  and  stator  are  coupled  together  through  a  mutual 
inductance  La£,  which  may  be  represented  pictorially  together 
with  the  respective  leakage  reactances,  as  shown  in  Fig.  II-1-3. 


Lad 


Figure  II-1-3 
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The  current  may  be  considered  as  the  magnetizing 
current  required  to  produce  Rated  Voltage  on  the  Air-Gap 
line  as  shown  in  Figure  X 1-1-4 „ 


Figure  II-1-4 

Thus  in  accordance  with  the  turns  ratio  established  in 
Section  II-1-3 

( 2/ 3  If0  NfAa)  xads  =  v  peak/phase  (II-1-42) 

.  o  Na/Nf  =  a  =  V  peak/phase 

2/ 3  Ifo  xads  (II-1-43) 

Therefore  the  base  rotor  quantities  become 

VBr  =  v  peak/phase  (II-1-44) 

a 

IBr  =  a  I  peak/phase  (II-1-45) 

This  transformation  results  in  a  fictitious  base  rotor 
field  voltage,  as  it  assumes  that  the  rotor  operates  on  a 
A„C.  basis o  This  factor  is  not  thought  to  be  significant 
however  since  the  transformation,  as  will  now  be  shown, 
affects  only  the  form  of  the  field  equations. 


. 
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Vf  IfRf  +  Lff  d  If 


aVf  = 


/3  If  iia 

/ajj 


dt 
■L2  3 


/  2  R- 


(II-1-46) 


.2  3 


/  2  L 


ff 


d_ 

dt 


*/3  I- 


L 


/; 


aVf  = 


y3  Jf 


r.2  3 


/; 


/ 2  Rf  +  I  a2  3/  x 


2  Afl  +  x 


ads 


W. 


d«  /3  If 


'gF 


Expressing  the  equation  in  per  unit 


aV 


,[2/>  ^ 


V  peak/phase  I  peak/phase 


a2  3/2  xfl  +  x 


ads 


Wo  XBs 


3/ 2  Rf 


X. 


Bs 


d  2  I 


/3  If/a 


dt  I  peak/phase 


WoVf  =  (IfWoRf>  +  <X1  +  Xads'  d_ 

dt 


(II-1-47) 


Where  1^  p,u.  = 


‘/3  If/a 
I  peak/phase 


R  r;  PoU, 


=  *2  3 


/2  Rf 


X 


Bs 


=  a2  3/  x 


2  Afi 


X 


Bs 


Xfl  p*u. 


. 
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Section  II-1-5  Generator  Saturation 


The  voltage  behind  the  armature  leakage  reactance  varies 
along  the  open  circuit  curve  in  the  region  near  rated  voltage. 
Thus  in  this  range  of  voltages,  the  deviation  between  the 
o.c.c.  and  the  air  gap  line  is  not  large  (Figure  II-1-6) , 
therefore  the  former  voltage  is  assumed  to  vary  along  the 
air  gap  line. 


Figure  II-1-7 
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It  is  necessary  to  obtain  the  voltage  behind  the  quad¬ 
rature  reactance  to  represent  the  salient  pole  machine  (6) 
(Figure  II-1-7) 


E 

q 


xd(xd 


I4rX  , 

f  ad 


-  Vxd 


-  V 


(II-1-48) 


As  the  direct  axis  reactance  varies  as  a  function  of  the  unit's 
load,  the  saturated  value  of  is  used  in  the  above  equation. 


Section  II-1-6  The  Voltage  Regulator 

Equation  II-1-49,  describing  the  voltage  regulating 
system  in  per  unit  values,  is  obtained  from  the  thesis 
presented  by  0,  W,  Hanson,  reference  4, 

v’£  =  K^e  "t/T  +  K2Ek  +  K,  ERdt  (II-1-49) 

where  ER  =  Vact  -  Vref,  Kl  =  1,  K2  =  5,  K  =  0.2,  T  =  10 

Due  to  the  difference  in  the  method  of  simulation ,  V'^  is 
equal  to  WQVf  in  the  equations  as  described  by  the  author. 

The  term  Epe"t/T  represents  the  temporary  over¬ 
regulation  of  the  generator's  field  voltage.  This  is  a 
non-linear  term  and  to  represent  it,  equation  II-1-49  is 
approximated  by  equation  II-1-50  (in  the  s-domain) . 


wovf  = 


(K-l  +  K2)  -  K1/T 
s+l/T 


K. 


R 


(II-1-50) 


Applying  a  step  change  (E^)  in  equation  II-1-50,  equation 
II-1-49  is  obtained.  As  the  system  is  tested  by  applying 
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step  load  changes,  which  approximately  result  in  step  voltage 
changes,  the  latter  equation  is  a  valid  approximation. 

It  is  noted,  in  the  reference's  computer  schematic,  that 
this  term  (Ere“ t'/’r)  is  represented  as  a  delay.  It  is  further 
noted  that  Mr.  Hanson  attained  5%  agreement  between  his 
computer  simulation  and  field  results. 

On  the  above  basis,  equation  11-1=50  is  accepted  and 
utilized  in  this  study. 

Section  II-1-7  The  Generator  and  Voltage  Regulator  Equations 
The  damper  windings  are  located  on  the  rotor;  therefore 
they  have  the  same  appearance  to  the  stator  as  dees  the 
rotor  circuit.  The  equations  are  now  expressed  in  per  unit 
and  written  for  generator  action. 


Ad  = 

Zf  Lad  ”  Id(Lal  +  Lad)  +  IkdLad  (7) 

(II-1-51) 

%  - 

^kqLaq  “  "^q^al  +  Laq^ 

( II-1-52) 

A  kd 

=  IfLad  “  xdLad  +  Ikd(Lad  +  Lkd) 

(II-1-53) 

^  kq 

Ikq^Laq  +  Lkq^  ~  IqLaq 

(II-1-54) 

II 

C1 

> 

W  +  s\  ~  I  R 

d  q  q  a 

( II-1-55) 

ll 

> 

— W  ?V  +  -  I,R 

q  d  d  a 

( II-1-56 ) 

Vkd  : 

IkdRkd  +  S^kd  =  ° 

(II-1-57) 

Vkq  : 

=  I  R  +  sX.  =0 

kq  kq  kq 

(II-1-58) 

Vf  = 

If Rf  +  s 7^  ^ 

(II-1-59) 
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Vf  =  ER  !(K1  +  K2>  - 

k1A 

(II-1-50) 

L 

s  + 

1/T  J 

E  =  IX  .  -  I,  fXa  - 
q  fa-4'  d  L  d 

Xq]  + 

Xkd  xad 

(II  -1-60) 

T  =  7\  I  -  T\ -I  _  (9) 

u  q  q  d  d 

(II-1-61) 

vt  =^Vd2  +  vq2 

(II-1-62) 

Approximations  made  to  the  above  equations 

1.  The  terms  s  Tv  and  s  ,  are  neglected  as  the  resulting 

q  Q 

error  is  less  than  one  per  cent.  (9) 

2.  The  armature  resistance  is  considered  to  be  negligible o 


It  is  noted  that  the  term  I.  J  ,  is  added  to  equation 

kd  ad 

II-1-48  to  produce  equation  II-1-60  as  it  acts  to  hold  the 
voltage  behind  the  quadrature  reactance  constant. 

With  these  approximations,  the  following  equations 
are  derived  from  equations  II-1-50  to  II-1-62  and  used  for 
the  analogue  simulation. 
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o  r 


-  Wo(Vf 


IfRf)/s 


=  ^(Xad  +  Xfl> 


IdXad  +  IkdXad 


(II- 1-63) 


E  «  W  A,  -  IfX., 
gp  o  f  f  fl 


=  (I.  -  1^  +  IV^)X 


kd;  ad 


Tv  ,  =  E__  -  I^X  , 

o  d  gp  d  al 


V  =  W  A 


Ikd  Wo^  kd  "  Egp 


X. 


kd 


(II-1-64) 

(II-1-65) 

(II-1-66) 


( II- 1-67) 
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W  /V ,  ^  =  W  (V.  j  —  I  R  ) 
O  kd  o  *  kd  kd  kd; 


=  -W„  I 


kd  \d 


Wo  “Iq(Xaq  +  XalJ  +  XkqXaq 


Vj  =  -W'A 


Xkq  Wo^kq  *  IqXaq 


X  4  X, 
aq  kq 


W°V  -  Wc(Vkq  -  IkqRkq) 


=  Wo\< 


WoVf  =  ER  !<ki  +  K2) 


-  Kx/T  +  K3 

c4-  1  /<V 


E  =  “  I  fx,  -  X  1+  I,  ,x  , 

q  f  ad  d  L  d  kd  ad 

PL  =  WTU  =  Vq  Iq  +  VdId  (8) 

Vt  ='Vvd2  +  Vq2  (8) 


(II-1-68) 

(II-1-69) 

(II-1-70) 

(II-1-71) 


(II-1-72) 

(II-1-50) 

(II-1-60) 

(II-1-73) 

(II-1-74) 


Section  II-1-8  Conclusion 

It  was  found  that  Park’s  equations,  which  represent 
the  stator  on  a  per-phase  basis,  reduce  the  current  in  the 
rotor  and  increase  the  rotor  impedances  by  a  factor  of  3 / 
This  multiplication  factor  does  not  affect  the  equations 
(stator  or  rotor)  as  it  leaves  the  flux  linkages  unaltered. 
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As  it  is  desired  to  represent  the  generator  on  a  per- 
phase  basis,  the  above  multiplier  was  applied  to  the 
respective  rotor  quantities.  They  were  expressed  in  terms 
of  the  stator  equations  through  the  turns  ratio  of  the 
generator  rotor  and  stator. 

Section  II-2  The  Speed  Governor  Model 
Section  II-2-0  Introduction 

The  governor  and  controller  determine  the  turbine 

•  ' 

wicket  gate  position  as  a  function  of  speed.  The  speed  of 
the  unit  varies  as  the  integral  of  the  difference  between  the 
power  input  and  power  output  of  the  system. ( refer  to  Figure  II-2-1) 


Figure  II-2-1 


The  relation  between  the  power  and  speed  deviation  is 
developed  in  Section  II-2-1.  The  governor  transfer  function, 
relating  the  speed  deviation  to  the  wicket  gate  position  is 
developed  in  Section  II-2-2.  The  controller  action  is  super¬ 
imposed  in  Chapter'll,  Section  II-5-1. 
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Section  II— 2—1  The  Power— Speed  Relationship 

The  unit's  speed  varies  as  a  function  of  the 
between  the  turbine  and  load  torque. 

tt  “  tl  =  WR2  d(wm) 
g  dt 

Tt  -  Tl  (2/p)WR2  d ( 2 trf e ) 

Trr  " ^  dt 


Assuming  Rated  KVA  =  Rated  KW  = 
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Therefore  equation  II-2-3  becomes 


f  p.u. 


f  p.u.  2H  df  p.u. 

dt 


difference 


(II-2-1) 


(II-2-2) 


/  T  T_  T|_  "3  \ 


(II-2-4) 


(II-2-5) 


( II-2-6 ) 
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Expressing  equation  II-2-2  in  per  unit 

PT  "  PL  =  f  2H  ££  (II-2-7) 

3t 

The  turbine  and  load  power  are  functions  of  system  frequency, 
refer  to  Figure  II-2-2  which  shows  their  variations  when 
the  input  power  to  the  turbine  is  held  constant.  (10) 


Power 


Figure  II-2-2 

Expanding  PT  and  PL  in  a  truncated  Taylor  series  about 
synchronous  speed  equation  II-2-8  is  obtained  which  is 
expressed  in  per  unit  quantities. 

Pt  -  P-L  +  r®t  dPi~1  df  =  f  2H  M  (II-2-8) 

-  j  dt 

The  magnitude  of  the  term  dPfc  is  contingent  on  the 

“df 

synchronous  speed  of  the  turbine.  For  low  and  high 
synchronous  speeds,  its  values  are  approximately  equal 
to  “Id  and  -0.6  respectively. 

The  magnitude  of  the  term  dP^  depends  on  the  nature 

“cFF 


of  the  load  and  the  voltage  frequency  relationship  of  the 
voltage  regulator.  The  load  contribution  to  the  above  term 
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varies  between  0  and  2  depending  on  wheather  the  load  is 
principally  ohmic  resistance  or  turbo-machines,  respectively. 

In  the  simulation,  the  damping  action  of  the  load  is 
absorbed  in  the  term  dP^  and  then  the  load  is  considered 

df~ 

to  be  independent  of  frequency.  The  voltage  regulator,  however, 

is  simulated  directly,  thus  its  contribution  to  the  term  dP^ 

is  not  considered. 

Therefore  dP  dP'^ 

-0  o  6  >  -  >-3.1 

dF~  “dF" 

The  term  f  2H  df  in  equation  II-2-8,  which  is  equal  to 

dt 

2H  df  +  Af  2H  df ,  is  approximated  by  2H  df.  This  is  justified 
dt  dt  dt 

as  -A  f  2H  df  is  small  when  considering  speed  deviations  of 

dt 

less  than  .04  per  unit.  Small  deviations  in  speed  are 

considered  as  it  is  desired  to  analyze  the  system  under 

normal  operating  conditions.  Thus  equation  II-2-8  becomes 

(P  -  P'  )  -  P,  Af  =  2H  df  ( II-2-9 ) 

1  "  d  3t 

Where 

pd  =  dPt  -  ae-j 

af  ~ar 

is  the  developed  turbine  power  at  rated  speed 
P  •  2_  the  load  output  power  at  rated  speed  plus  the  change 
in  power  resulting  from  the  voltage  regulator  frequency 
relationship. 
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is  the  self  regulation  of  the  turbine 

is  the  self  regulation  of  the  load 

Section  II-2-2  The  Governor  Transfer  Function 

The  governor's  function  is  to  hold  the  unit  at  synchronous 
speed.  It  accomplishes  this  by  controlling  the  wicket  gate 
position,  and  therefore  the  developed  power,  as  a  function 
of  speed.  The  operation  of  the  governor  as  depicted  in 
reference  II,  is  shown  in  Figure  II-2-4.  Simplification  of 
the  governor  is  then  made  to  the  extent  shown  in  Figure 
II-2-6.  In  the  development  of  the  equations,  the  wicket 
gate  movement  is  assumed  to  be  slow  enough  that  the  governor 
is  capable  of  sustaining  its  oil  pressure.  The  block  diagram 
of  the  governor  action  is  shown  in  Figure  11-2-3. 

The  sensing  and  detection  of  speed  error  is  achieved 
in  the  ball  head  assembly  by  comparing  the  centrifugal  (11) 
force  of  the  flyweights  against  the  force  exerted  by  the 
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d F“ 
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df 


Figure  II-2-3 
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speeder  spring  (Block  #1) »  The  movement  of  the  flyweights 
is  transmitted  to  the  turbine  wicket  gates  through  mechanical 
linkages ,  pistons,  valves  and  an  oil  pressure  system  (Block  #2). 
The  stabilizing  features  are  represented  in  Blocks  3  and  4. 

Block  #3  represents  the  permanent  speed  droop.  This  is  the 
proportional  reduction  of  the  governor  speed  setting  in  a 
given  ratio  as  the  turbine  gate  opening  increases.  It  is 
effected  through  the  restoring  mechanism  and  is  required  to 
produce  stable  load  sharing  between  units.  As  the  permanent 
speed  droop  is  decreased  the  load  change  increases  for  a 
given  speed  change.  Block  #4  represents  the  temporary  speed 
droop.  Due  to  the  water  inertia,  a  delay  is  incurred  between 
the  wicket  gate  opening  and  the  developed  turbine  power. 

As  the  former  is  proportional  to  the  integral  of  the  pilot 
valve  position,  it  is  necessary  to  provide  a  large  permanent 
droop  or  another  means  of  compensation  to  counteract  over 
speed  regulation.  Therefore  compensation  is  provided  in 
the  form  of  temporary  speed  droop.  This  results  in  an 
initial  reduction  of  the  governor9 s  speed  setting  which 
gradually  decreases  to  zero.  Block  #5  is  discussed  on 
page  39 0 

Transfer  Function  of  the  Ball  Head  Assembly  (13) 

This  transfer  function,  B(s)f  is  represented  as  a 
gain  constant  B.  This  is  justified  as  represents  a 

second  order  system  whose  natural  frequency  is  much  higher 
than  the  natural  frequency  of  the  governoring  loop. 
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Figure  II-2-6 
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Governor  Gain 

The  gate  opening  time  determines  the  gain  between 
the  servo  motor  and  the  pilot  valve .  It  is  assumed  that  a 
frequency  error  equal  to  the  permanent  speed  droop  will 
completely  open  the  pilot  valve  when  compensation  is  not 
applied.,  Thus  the  integral  of  this  gain  times  the  speed 
droop  over  the  interval  of  time ,  equal  to  the  gate  opening 
time,  will  completely  open  the  wicket  gates 0 

In  the  forthcoming  analysis,  all  values  are  expressed 
in  per  unitQ  As  is  shown  by  Figure  II=-2— 3,  the  position  of 
the  pilot  valve  is  determined  by  the  frequency  deviation 
and  droop  signals „ 

Frequency  deviation  signal  BAf 


Permanent  speed  droop  signal  G 


Temporary  speed  droop  signal;  refer  to  Figure  II-2-7  where 


Pilot 

Valve 


f 


y 


Figure  II-2-7 
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Cs(G  -  y)  =  Dy 


where  d  =  1/Tr 

The  pilot  valve  position 


xp  =  B 


Af  -  sS9  G 

BTs+I/Tr) 


The  servo-motor  valve  position 


X  =  K 9  X 
v  P 


(II-2-10) 


(II-2-11) 


(II-2-12) 


Let  S  !t  S 
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6'  S  ,  K 9  B 

P  =  P 

~B  ■ 
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A 


(II-2-13) 


The  turbine  wicket  gate  position  is  a  non-linear  function  of 
the  servo-motor  valve  position  (refer  to  Figure  II-2-3) , 
(Block  #5) » 


G  =  F(XV) 


(II-2-14 ) 


The  non-linearity  of  the  wicket  gate  position  as  a 
function  of  the  servo-valve  position  was  found  to  have  little 
effect  on  the  stability 0  An  analysis  is  made  on  this  effect 
after  the  development  of  the  hydro  transfer  function  (see 
Chapter  II,  Section  II-3-13) 0 
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Thus  the  simulated  equations  representing  the  speed 
governoring  models  expressed  in  per  unit  are 

pt  -  Pd  Af  -  P"  =  2H  df  ( II-2-9 ) 

dt 


( II-2-15 ) 


Section  II-2- 3  Conclusion 

The  per  unit  torque  equations  may  be  replaced  with  per 
unit  power  equations  and  the  rotating  inertia  power  may  be 
assumed  constant  considering  speed  deviations  of  less  than  4%« 


Section  II-3  The  Hydro  Model 
Section  II-3-0  Introduction 

The  transfer  function  relating  the  turbine  output  power 
to  the  wicket  gate  position  is  determined  in  this  section 
for  systems  with  and  without  a  surge  tankQ 

Turbine  output  power  PT  =  Q  H  '*7  (I I- 3-1) 

Qt  =  Q(H,G,Wc,Ec)  =  GCd~\f2gKt  550 

Flow  rate  ( II— 3—2 ) 


Ht  =  H<Qt'V 


Piezometric  head 


G 

Wc  -  W(H) 
Ec  =  E (H) 


Efficiency 

Effective  length  of  the  water 
column  as  seen  by  the  turbine 
Wicket  gate  position 
Water  compressibility 
Elasticity  of  the  conduit 
Coefficient  of  Discharge 
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As  may  be  seen  by  the  above  functional  relationships, 
the  determination  of  the  turbine  output  power  as  a  function 
of  the  gate  position  is  a  complex  problem  as  it  involves 
a  large  number  of  interrelated  variables 0  Therefore,  to 
obtain  a  solution,  the  problem  is  sectioned  in  the  following 
manner c 

Ae  Systems  without  a  surge  tank  are  considered  first® 

10  Idealizations  are  made  to  the  problem  in  Sections 

) 

II-3-1  and  II-3-2. 

2,  The  relationships  between  the  variables  are 
established  in  Sections  II-3-3  to  II-3-5® 

30  Approximations  are  made  to  the  equations  in 
Section  II-3-60 

4a  The  solution,  to  the  approximated  equations,  is 
outlined  and  interpreted  in  Section  II-3-7® 

5,  A  transient  hydraulic  performance  -  transmission 
line  analogy  is  made  in  Section  II-3-8® 

6 0  The  transfer  function  is  obtained  for  the  fundamental 
component  in  the  pressure  and  flow  rate  variations  in 
Section  II-3-90  The  higher  harmonic  components  are 
considered  in  Section  II-3-12* 

Bo  Systems  with  a  surge  tank  are  considered  next, 
lo  Idealizations  are  made  in  Section  II-3-lOo 
2.  The  transfer  function  is  obtained  in  Section  II- 3-11. 

The  effect  of  the  turbine  wicket  gate  non-linearity  and 
the  higher  harmonics  in  the  pressure  and  flow  rate  variations 
are  considered  in  Section  II-3-13  and  Section  II-3-12  res¬ 
pective  ly  o 


' 
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Section  II-3-1  The  Idealization  of  the  Conduit 

The  conduit  from  the  reservoir  to  the  turbine  is  assumed 
to  lie  in  a  straight  linec  If  the  conduit  length  remains  the 
same,  this  assumption  has  no  effect  on  the  problem  as  seen  by 
the  turbine  (refer  to  Figures  II-3-1  and  II-3-2) „ 


Figure  II-3-1 
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Section  II  3-2  The  Idealization  of  the  Piezometric  Head 


Assuming  no  friction*,  and  steady  state  conditions^  the 
pressure  and  flow  rate  in  the  conduit  may  be  determined 
from  Bernoulli 's  equationQ 


=  P/-TT+Z  +  Up2/2g 


H 


(II-3-3) 


o 


The  conditions  at  the  turbine  are  given  by  equation  II-3-4 
Ho  =  p/^  +  Up2/2g  =  Hfc  +  Up2/2g  (II-3- 

The  water  velocity  in  the  turbine  has  been  found  to  follow 
the  relation  as  expressed  by  equation  II-3-20 

Qt  =  ^d^^^t  (ii _ 3__ 

Further*,  the  velocity  head  between  the  reservoir  and 
turbine  is  assumed  to  be  negligible  in  accordance  with 
reference  #  17 „ 


Change  in  Conduit  Area  and  Fluid  Compressibility 


The  change  in  area  of  the  conduit  as  a  function  of  the 
change  in  pressure  may  be  determined  by  cutting  a  section 
out  of  conduit  (refer  to  Figures  II-3-3  and  II-3-4)0 


T 


dx 


Figure  II-3-3 


3  T 

\rU  "to  nox.lQA]  j3  b  &:  j/,.  .t:oj 
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dx 


AT  =“&dH  R  =  "&dH  D  /2 

P 


AL  =  TT  D  dH 

Jr 


(II-3-5) 


(II-3-6) 


Ad  =  =  Dp2  "S’  dH 

P  _____ 

2  t'E 


AAp  =  ^/4  ^(Dp  +ADp)2  -  Dj 


=  ty4 


2ADpDp  +ADp2] 


(II-3-7) 


t*/2d  Ad 

p  p 


Aa„  =  D  3  ~5dH 
P  p 

4  t'E 


(II-3-8) 

(II-3-9) 


The  change  in  volume  over  the  element  dx  is  given  by  equation 
II-3-10 

dvp =  DP3y>dn 

_  dx  (II-3-10) 

4  t'E 

The  change  of  the  volume  of  the  fluid  is  given  by  equation 
II-3-11, 

K  =  -V  dH  IS  =  -V  dP 
dV  dVc 


dvc  =  -V  ~6  dH 
K 


(II-3-11) 
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The  total  change  in  volume  over  the  element  dx  is  given 
by  equation  II-3-120 

dV^  +  dV„  =  VtfdH  +  TfD^3xdH  dx 
v  c  — p 

~i  ~t '  E 


+  dVc 


(Apdx)lTdH 


+  D 

P 

Fe 


(II- 3-12) 


Section  II-3-4  The  Relation  Between  the  Flow  Rate  Variation 

and  the  Change  in  Volume 

The  equation,  relating  the  flow  rate  change  to  the  volume 
change,  may  be  determined  by  referring  to  Figure  II-3-5. 


dx  +  "b  U 

P 

Tt 


dt 


x2 


dt 

X1 


(II-3-13) 


Change  in  flow  rate  =  A  du  (II-3-14) 

p  P  ' 

The  required  change  in  volume  due  to  the  flow  rate 

change  =  A  du  dt 
P  P 


(II-3-15) 
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Assuming  bu^  dt  is  constant  over  the  element's 
b  x 

the  required  change  in  volume  =  b  n  A  dt  dx 

P  p 


length , 


(II- 3-16) 


b  x 

As  the  change  in  flow  rate  is  the  resultant  effect  of 
the  change  in  volume,  one  may  equate  equations  II-3-12  and 


II-3-16. 


A  dx  dH  15 
P 


1  +  E> 
K 


t'E_ 


-A  "b  U  dx  dt 
P  P 

"b  X 


dH  *6 
dt 


1  +  E) 
K 


b  x 


(II-3-17) 


dH  =  bH  +  bH  dx  (II-3-18) 

dt  bt  bx  dt 


IS 

^H 

bt 

+ 

UD  bH~ 
F  bx 

1 

Q  1 

+ 

-  bUp 

(II-3-19) 

L  tEj 

b  X 

bH 

+ 

up  M 

—  _ 

1 

_  ^>u 

bt 

P  bx 

'I 

+  °P 

% 

tE" 

-J 

b  H 
bt 

+ 

UpbH 

=  -si 

9 

bup 

(II-3-20) 

bx  "6  x 


Section  II-3-5  The  Relation  Between  the  Pressure  Variation 

and  Water  Inertia 

The  equation,  relating  the  change  in  flow  rate  to  the 

change  in  pressure,  may  be  determined  by  referring  to 

Figure  II-3-6.  bH 

H  N-  dx  H  +  dX 


Figure  II-3-6 
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F  =  ’S  dH  A 

P 


=  -  ma  = 


dt  B  t 


=  Ua, 


b  U 


+ 


BH  dx 

+  BH 

dt 

-  BH 

jbx 

B  t 

'x 

X2 

dt 

^pdx) 

dup 

"dt 

’  dx 

BU 

P 

U  7 

+  P 

:  dt 

Tt 

dt 


J  P 

B  x 


(II-3-21) 


(II-3-22) 


(II-3-23) 


Assuming 

equation 


bH  dt 
Bt 

II-3-21 


is  constant  over  the  elements  length, 
becomes  equation  II-3-24. 


rAp||  dx  =  -eApdx  pup  +  up  -8 Up' 

[  Bt  B  X  _ 

g  H  =  “[^up  +  up  BUpl  (II-3-24) 

l  B  t  B  x] 

Section  II-3«-6  The  Approximation  of  the  Equations 

The  two  equations „  relating  the  pressure  and  flow  variations 
in  the  conduit  under  transient  conditions,  are  as  follows: 


b  h 

+  up 

B  H  = 

-'sL  ^up 

dt 

Bx 

g  _ 

(II-3-20) 

B  x 

-  g 

BH 

B  x 

=  aup 

~ 

+  up  ^UP 

(II-3-24) 

B  x 

The  boundary  conditions  of  these  equations 
discontinuities  are: 

At  the  top  of  the  reservoir  H  =  0,  U 

At  the  turbine  Ut  =  G 

At  the  reservoir  conduit  junction  H  =  HQ, 


,  neglecting  pipe 


=  0,  Ar  =  QO 
CdV2g  H  (II-3-2) 

U  =  0 ,  Ar  =  OO 


r  VK  '  'Snr  **Xr •  ^  Trr  ->ff  ►  r  ■ 


' 

i  ..  .  c  j  );■  •  ..  -  i’T 


oo  «-  a  ,o  u  %0 


co  «  ,0  =  J  „  H  nox^onuc  aiubnoo  :iovi»a  ortt  3A 
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The  initial  condition  is  H  =  H  ,  due  to  the  idealization  of 
the  problem  made  in  Section  II-3-2. 

The  problem  may  be  simplified  by  approximating  equations 
if ”3-20  and  II-3-24  and  then  obtaining  two  second  order 
equations »  The  approximations  and  the  justification  of  them 
are  now  shown. 


BU 
"6  x 


an 

bt 


IL 


a»x 


= 

Bt 


BU 

__Z 

B  x 


Bt 


(II-3-25) 


As  b2  >>Up2 


an  = 
at 


-b2  ^U. 

g 


a  x 


( II-3-26 ) 


Substituting  equation  II-3-20  into  equation  II-3-24 

a»H 

9  ~ 


BH 

ax 


r(2V  - 


b‘ 


^  . 

VTsh  +  u 

Sf[ 

at 

b2 

"a  x 

b2) 

g 

- 

at 

(II-3-27) 


BU. 


^  -  g  BH 

at  ax 

Two  second  order  equations  are  now  obtained  by 
differentiating  equations  II-3-26  and  II-3-28  and 
substituting. 


(II-3-28) 


■ 
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Differentiating  equation  II-3-26 


2 

B  ZH 
dt2 

=  ^ 

g 

»up 

'hx^bt 

d2H 
ex  at 

=  =£ 
g 

qJ 

G 

ts 

Bx2 

Differentiating  equation  II-3-28 
-g  B2h 

Bfbx 


B2Up  = 

dx^t” 


-g 


( II-3-29) 


(II-3-30) 


(II-3-31) 


(II— 3-32) 


From  equations  II-3-29  to  II-3-32*  the  following  expressions 
are  derived,, 

b2  =  Vh  (II-3-33) 

”dx2  dt2 


b^  Vup  =  Vup 

x2  Bt2 


(I 1-3-34) 


The  problem  is  now  reduced  to  that  of  solving  equations 
**”3-33  and  II-3-34*  subject  to  the  boundary  and  initial 
conditions  given  at  the  beginning  of  the  section,, 

Section  II” 3-7  The  Solution  and  Interpretation  of  the  Equations 
The  general  solution  of  equations  II-3-33  and  II-3-34 
is  known  to  be  (15) s 

H  -  H0  =  F ( t  +  x/b)  +  f  (t  -  x/b)  (II-3-35) 

Up  ”  Upo  =  g/b  “  x/b)  “  F  (t  +  x/b)] 


■-  -  r 


qis^n  b  it  j  .  -'iri  jS  erfl' 
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QP  ”  Qpo  =  ^  -f(t  “  x/b)  "  F(t  +  x/b>]  (II-3-36) 

“b“ 

In  deriving  the  equations,  x  and  U  were  measured 

XT 

positive  from  the  reservoir.  In  the  application  of  the 
solution,  it  is  more  convenient  to  measure  x  positive  from 
the  turbine  as  this  is  where  the  disturbance  occurs. 

Rewriting  the  equations  on  this  basis 
H  -  H0  =  F(t  ~  x/b)  +  f(t  +  x/b)  =  F  +  f  (II-3-37) 

Qp  -  Qp0  =  Apg  [F(t  -  x/b)  -  f(t  +  x/b)]  =  -Apg  [f  -  f] 

“F*  —5“ 

Although  equations  II-3-37  and  II-3-38,  do  not  satisfy 
the  boundary  conditions,  they  are  still  useful  as  they  enable 
a  physical  picture  of  the  hydraulic  transient  performance  to 
be  obtained.  When  the  hydraulic  performance  is  understood, 
it  is  possible  to  use  a  graphical  method  to  solve  the 
equations  subject  to  the  boundary  conditions. 

To  gain  an  insight  into  the  hydraulic  transient  per¬ 
formance,  the  interpretation  of  the  equations  and  their 
graphical  solution  is  outlined  for  rapid  gate  closures. 
Following  this,  the  graphical  solution  is  outlined  for  slow 
and  oscillatory  gate  movements. 

Physical  interpretation  of  the  equations  for  rapid  gate 
closures  is  shown  in  Figures  II-3-7  to  II-3-12. 
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Initial  Conditions 

f||,j  Instantaneous  Gate 
ill'  Closure  Assuming 
'h!  No  Friction 


II-3-7 


II-3-9 


II-3-10 


II-3-11 


II-3-12 


r 


52 


At  time  t  =  4L+  ,  the  process  starts  over  again. 

a 

At  time  t  =  0,  the  system  is  undisturbed  (Figure  II-3-7) 0 
The  gate  is  assumed  to  close  instantaneously  at  time  t  =  0+ 
(Figure  II-3-8) „  As  a  consequence  of  the  inertial  force, 
a  pressure  increase  ensues,  causing  the  conduit  to  expand 
and  the  water  to  compress.  This  increase  in  pressure, 
called  an  F  wave,  arrives  at  the  reservoir  at  time  t  =  L/b+, 
and  stops  the  flow  into  the  pipe.  Due  to  the  increased 
pressure  in  the  conduit,  water  flows  out  of  the  conduit 
until  time  t  =  2L/b+  (Figure  11-3— 10)  .  The  force  of  the 
inertia  allows  the  water  to  continue  to  flow  out  of  the 
conduit  until  time  t  =  3L/b+ ,  thus  generating  an  F  wave 
(Figure  II-3-11) .  At  this  time,  suction  pressure  is  built 
up,  and  an  f  wave  is  generated  causing  water  to  flow  back 
into  the  conduit  (Figure  II-3-12) .  This  cycle  repeats 
itself  and  results  in  an  oscillation  called  water-hammer. 

During  the  time  intervals,  (0+  +  2nL/b)  <  t  <  (L/b+  +  2nL/b) , 
an  F  wave  travels  from  the  turbine  to  the  reservoir.  During 
the  time  intervals,  (L/b+  +  2nL/b)  ^  t  <  (2L/b+  +  2nL/b) , 
an  f  wave  travels  from  the  reservoir  to  the  turbine.  This 
is  expressed  in  equations  II-3-39  to  II-3-44. 


Ht 

(L/b+  +  2nL/b)  -  Hr  (0+  +  2nL/b) 

=  F 

( II-3-39) 

Hr 

(2L+/b  +  2nL/b)  -  Ht (L/b+  +  2nL/b) 

=  f 

(II-3-40) 

Qt 

(L/b+  +  2nL/b)  -  Qr  (0+  +  2nL/b)  = 

-V  F 

(II-3-41) 

°r 

(2L/b+  +  2nL/b)  -  Qt  (L/b+  +  2nL/b) 

b 

=  -Apg  f 

(II-3-4?) 

F 


•-  1  ) 


. (SI-C-II 

til 


(d\J tl  +  +d\J)  H  -  (d\Jn£  +  d\+J£) 


(lfr-£-II)  1  -  p^A-  *  (d\JnS  +  +0)  XD  -  (d\J n£  +  +d\J) \,Q 
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Ht  (L/b+  +  2nL/b)  -  H  (0  +  +  2nL/b) 


=  ~  b  [Qt (L/b+  +  2nL/b)  -  Q  ( 0  +  +  2nL/b)l 

Hr  (2L/b+  +  2nL/b)  -  Ht(L/b+  +  2nL/b) 

=  ~b  [Qr(2L/b+  +  2nL/b)  -  Qt  (L/b+  +  2nL/b)] 

Apg 


(II— 3-43) 


(II — 3—44) 


As  the  conditions  are  the  same:  between  L/b+  and  2L/b , 
2L/b  to  3L/b,  equations  II— 3— 43  and  II— 3— 44  are  rewritten 
as  II- 3-4 5  and  II-3-46. 


Ht2(n+1)  L/b  -  II  (2n+l)  L/b  =  -b_ 

r  Apg 


Hr(2n+1)  L/b  -  Ht(2n)  L/b  =  -b_ 


A  q 


[Qt  2(n+l)  L/b  -  Qt(2n+1)  L/b)] 

(II-3-45) 

[Qr(2n+1)  L/b  -  Qt(2n)  L/b)] 

(II-3-46) 


The  graphical  solution  of  the  equations  II-3-45  and 
II-3-46,  subject  to  the  boundary  condition  Qt  =  CdG V29H / 
is  shown  in  Figure  II-3-130  As  the  gate  is  assumed  to  close 
instantaneously,  the  flow  Qt  =  0. 


ter  e 
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The  resulting  pressure  head  oscillation  at  the  turbine  is 
shown  in  Figure  II-3-14, 


time 


Figure  II-3-14 


The  graphical  solution  of  the  equations  for  slow  gate 
closures.  The  gate  closure  time  is  divided  into  n(2L/b) 
intervals  and  a  step  approximation  to  the  closure  is  made 


(refer  to  Figure  II-3-15) 


Gate 

Position 


The  graphical  solution  is  shown  in  Figure  II-3-16, 


The  head  and  flow  rate  variations,  as  a  function  of  time, 
are  shown  in  Figures  II-3-17  and  II-3-18. 
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Figure  II-3-17  Figure  II-3-18 

The  graphical  solution  of  the  equation  for  oscillatory 
gate  movements 0  The  assumed  gate  oscillation  is  shown  in 
Figure  II-3-19o  The  graphical  solution  is  shown  in  Figure 
II-3-20  and  the  pressure  and  flow  variations,  as  a  function 
of  time,  are  shown  in  Figures  II-3-21  and  II-3-22. 
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Figure  II-3-20 


Figure  II-3-21 


Figure  II-3-22 

The  form  of  the  pressure  and  flow  rate  variations  for 
different  types  of  gate  movements  has  now  been  established. 
Utilizing  the  results,  for  the  case  of  oscillatory  gate 
movements,  a  transient  hydraulic  performance  -  power  trans- 
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mission  line  analogy  is  made  in  Section  II-3-8.  From  this 
analogy ,  the  relationship  between  Q  and  H  is  established. 

Section  II-3-8  The  Transient  Hydraulic  Performance  — 

Transmission  Line  Analogy 

By  sectioning  (17)  the  conduit  (Figure  II-3-23)  ,  equations 
■^“3“26  and  II-3-28,  may  be  written  as  finite  different  equations. 


Figure  II-3-23 


dH 

at 


r L  (n+1)  +  L  (n)  ,  t 

=  -b2 

Tq  L (n+1) 

,  -  Q  L (n)  ,t 

L  2  J 

A  g 

L  L  (n+I) 

~  L  (n)  J 

(II-3-47) 

dH 

at 


[■ 


L ( n+ 1 )  +  L(n)  ,tl  =  -b2 


J 


V 


aqp 

A,  L 


(II-3-48) 


L (n+1)  +  L (n) 

- - 2 - 


A^g  TH  L (n+1) ft  -  H  L(n) 
1  L  L ( n+ 1 )  -  L  (n) 


(II-3-49) 


dQ^  L (n+1)  +  L (n) ,t 

_ _  2 

dt 


-Apg  All 
AL 


(II-3-50) 


The  pressure  and  flow  variations,  over  the  Section  ,  are 
lumped  and  represented  as  capacitance  and  inductance  respec¬ 
tively  (refer  to  II-3-24) . 
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AQ 


T 


AHLn  -*l 

^TOTryrjp - 
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Figure  II-3-24 


The  turbine  sees  the  circuit  shown  by  Figure  II-3-25, 
as  the  reservoir  offers  zero  impedance  to  the  variations  and 
appears  as  a  D0C„  source 0 


-Turbine 


Figure  II-3-25 

To  make  the  transient  hydraulic  performance  analogous 
to  the  power  transmission  line,  it  is  necessary  to  have  an 
undamped  source  equivalent  to  voltage »  This  is  attained 
by  considering  sustained  pressure  oscillations  resulting 
from  undamped  gate  oscillations. 


Equivalence 


Hydraulic  System 

H  (head) 

-AQ(flow  rate) 
1/Apg 


Apg/b2 


w 


Power  Transmission  Line 

E (voltage) 

I ( current) 

L' (inductance) 

C (capacitance) 

$ ( frequency) 


■ 


. 
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Transmission  Line  Equations  (16) 

Es  =  ErcoshmL  +  IrZQsinhmL  AHfc  =  HrcoshmL  +  AQrZ0sinhmL  (II-3-51) 

Is  =  lrcoshmL  +  ErY0sinhmL  -  AQt  =  QrcoshmL  +  HrYosinhmL  (II-3-52) 
where 


Es 

sending  voltage 

AHt 

turbine  pressure  deviation 

Is 

sending  current 

AQt 

turbine  flow  rate  deviation 

Er 

receiving  voltage 

A  Hr 

reservoir 

pressure  deviation 

Ir 

receiving  current 

AQr 

reservoir 

flow  rate  deviation 

zo  =' 

~\lj  f  L’  =VL'/C  = 

’  j  f  C 

b/Ap 

9 

Yo  = 

b/®o 

m  = 

VTj  )2L'C  =  j 

Vb 

The  equivalent  'W  circuit  (refer  to  Figure  II-3-26) 


_  ^ 

— nnnno - 

8 

Tf 

tJ 

Figure  II-3-26 

8^  =  80sinhmL 


^  L j  7 

~Tr~  (II-3-53) 


8^  =  8QsinhmL 

coshmL  -  1 


b  j 
Ap9 


p  4 

1  —  (  L/b  )  t  (  L/b  )  t  oe«ooao*ua«see 

Jl  §1 

-TL/bT - + — l  T75) — =H~E7bT5"+  — - 

Tr~  31  61 

( II-3-54)  J 
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The  natural  frequency  of  the  governoring  loop  (  ^  ) , 
is  normally  less  than  0o7  radians  per  second.  Assuming  L/b 
to  be  less  than  0 „ 6 ,  the  impedances  offered  to  the  funda¬ 
mental  components  of  the  pressure  variation  are  given  in 
equations  II-3-55  and  II-3-56. 

%  =  b  j  (0.42  -  0.0122  +  0.0001  - + - ) 


3J.il 

Ar,S 


(II-3-55) 


1-0.0294  +  0.00026 

-0.21  +  0.0039 


( II-3-56 ) 


The  impedances  offered  to  the  third  harmonic  component 
of  the  pressure  variation,  assuming  =  0.7  and  L/b  =  0.6, 
are  given  by  equations  IT-3-57  and  II-3-58c 

3  =  b  3  (1.26  -  0.333  +  .0263  + - ) 

3  =  j  0.927  _b_  =  1.22  j  (3  "?  )  _L_  (II-3-57) 

Ap9  Ap9 


b 

"l  -  0.264 

+  .0209  - 

V 

f-0.63  + 

TOTH  - 

1. 37 

b 

2.18 

3 

A  g 

3  LApg 

(II-3-58) 

As  the  reservoir  offers  zero  impedance  to  the  A.c. 
variations,  the  capacitor  at  the  reservoir  is  shorted  out. 
Thus,  the  relationship  between  the  fundamental  component 
of  the  pressure  variation  and  the  flow  rate  variation  is 
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given  by  equation  II-3-59 
r 

t 


=  (-AQfc) 


(“AQt)  3^ 


*« 

"7 


+  i 


AHt  =  (-AQt)  j  fL/A  g 


1  +  (p)2  |L\2 


AHt  «  (-AQt)  j  >  L  ,  as  ( j  l  ) 2  IL)2  <.088  (II-3-59) 

Ap9  2  lb/ 


Expressing  equation  II-3-59  in  operator  notation,  one 
obtains  equation  II-3-60e 


£>H. 


/  L 

W 


(II- 3-60) 


The  relationship  between  the  third  harmonic  component 
of  the  pressure  variation  and  the  flow  rate  variation  is 
given  by  equation  II-3-61. 

AHt3  =  bK  3  3  ^Qt3* 


^Ht3  “  3  0.927  _b 

“l  -  6.^77 


Ap9 


(-  AQt3) 


*Ht3  -  (.aot3) 


[3.4j  (31)  _L_1 

l  v  j 


(II-3-61) 


The  impedances  offered  to  the  fundamental  and  third 
harmonic  components  of  the  pressure  variations  are  con¬ 
siderably  different  as  may  be  seen  by  equations  II-3-59  and 
II-3-610  Equation  II-3-59  shows  that  storage  elements  do 
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not  effect  the  relationship  between  &Q  and  for  the  fun¬ 
damental  frequency  component  Equation  II-3-61  shows  that 
the  inertia  and  storage  elements  act  as  a  tuned  circuit 
to  the  third  harmonic  frequency 0 

The  significance  of  the  water  storage  elements  and 
the  non-linearity  of  the  tt  circuit  impedance  depends  on 
the  magnitude  of  the  harmonic  components  in  the  pressure 
variations.  Thus,  it  is  pertinent  to  determine  the  source 
and  magnitude  of  the  harmonics,  as  there  is  no  point  in 
attempting  to  represent  the  above  factors  if  the  harmonic 
components  are  not  present  in  the  simulation. 

As  the  relationship  between  AH  and  6Q  in  the  conduit 
is  linear,  the  harmonics  are  created  by  the  non-linear 
boundary  condition  (equation  II-3-2)  .  If  equation  II-3-2 
can  be  linearized  for  small  gate  movements,  the  harmonic 
generator  is  removed,  and  equation  II-3-60  yields  a  correct 
model  of  the  hydraulic  performance. 

The  boundary  condition  is  linearized  and  the  hydro 
transfer  function,  based  on  a  rigid  conduit  and  an  incom¬ 
pressible  fluid,  is  developed  in  Section  II-3-9,  Further 
considerations  of  the  third  harmonic  component  are  made  in 
Section  II-3-12. 

Section  I I- 3-9  The  Hydro  Equation  for  Systems  Without  a 

Surge  Tank 

The  relation  between  pressure  and  flow  rate  in  the 
turbine  is  given  by  equation  II-3-2. 
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Qt  =  CdG  V 2 gH t 


dQt  =  CddGV2gHt  +  CdG  V^gH  dHt 


(II-3-2) 


2  Wt 


AQt  =  cd  AG  V2gHt  +  cd  G\2g^H 


(II-3-62) 


2  VHt 


Placing  equation  II-3-62  on  a  per  unit  basis 


&Qt 


(II- 3-6 3) 


Equation  II-3-63  has  time  varying  coefficients.  If  the 


m  \  can  be  assumed  constant 


over  a  given  range,  the  equation  is  linearized  over  this 
range .  This  may  be  accomplished  by  assuming  small  gate 
movements*  (17) 

The  change  in  pressure  for  large  gate  movements  was 

investigated  by  the  author.  It  was  found  that  the  head 

deviation  did  not  exceed  30%  for  a  system  with  2,000  ft. 

of  conduit.  Thus,  for  small  gate  movements,  it  is  justified 

to  assume  head  deviations  of  less  than  10%  U  «,  This  means 

that n\/h^/h"0  does  not  vary  more  than  plus  or  minus  5%. 

Therefore  equation  II— 3— 63  may  be  assumed  linear  for 

small  gate  deviations  as  the  coefficients  VHx./H0  and  (G/G  )x 

^  o  o 

(ViVHt)  are  almost  constant. 
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As  equation  II-3-2  is  linearized,  equation  II-3-60  is 
assumed  to  describe  the  hydraulic  performance  for  the 
reasons  stated  in  Section  II-3-8,  Placing  equation  II-3-60 
on  a  per  unit  basis 

Q0L  ^Qt 

w  s  “ 


AHt 

T 

o 


-T  s 
w 


(II-3-64) 


Substituting  equations  II-3-63  and  n-3-64 


=  2 


G  H 
o 


t 


G  H 

o 


&  G 


(II— 3-65) 


(II- 3-66) 


The  equation  relating  the  turbine  output  power  to  the 
head  and  flow  rate  variations. 
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pt  =  ^tHt  ^ 

~“55  JT 

dPt  =  [dHtQt  +  HtdQtJ 


550 
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Expressing  equation  II-3-67  in  per  unit 


A  p 


^HtQ. 


HtAQfc 


HoQo 


Ho  Q0 


(II-3-1) 

(II-3-67) 


(II-3-68) 


Substituting  equations  II-3-65  and  II-3-66  into  equation 
II-3-68 


Equation  II-3-69  has  now  been  obtained  in  the  Laplace  domain. 
It  will  now  be  expressed  in  the  time  domain,  by  taking  the 
Inverse  Laplace  transform,  to  consider  the  effects  of  the 
coefficients  assuming  a  step  change  in  gate  position. 

-1 
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(II-3-70) 


Tne  following  deductions  are  made  from  equation  II-3-70. 
me  power  output  first  decreases  and  then  exponentially 
rises  to  its  steady  state  value.  The  decrease  in  the  pressure 
head  tencs  to  de-stabilize  the  system  as  it  acts  tc  decrease 
the  power  output  and  increases  the  exponential  decay  time. 

The  term  GQ/G  acts  to  stabilize  the  system  as  it  decreases 
the  exponential  decay  time. 

According  to  the  rigid  conduit  and  incompressible  fluid 
theory,  as  may  be  realized  from  equation  II-3-66 ,  the  smaller 
gate  position  the  larger  is  the  head  deviation 
for  a  given  change  in  gate  position.  This  effect  becomes  of 
greater  siqnifi  car'-e  ns  ~  ~  ~  -  j i _ _ _ 

(Tw  increases)  .  However  as  the  length  of  the  water  column 

increases,  the  storage  capacity  of  the  conduit  also  increases 

and  as  a  result  the  pressure  variation  at  the  turbine  is  less 

se\ ere ,  On  this  basis,  coupled  with  the  fact  that  the  xe cove ry 

time  and  the  magnitude  of  the  head  deviations  are  inverselv 

related  for  a  given  gate  position,  it  is  assumed  that  the 

resultant  contribution  of  the  terms  G/G  and  E^/E  may  be 

o  t  o  J 

neglected.  Thus  equation  II-3-69  is  simplified  to  equation 


II-3-71. 


2  (1 


AG 


(II-3-71) 
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+  Initial  Conditions  (11-3-72) 
Go 


Section  II-3-10  The  Idealization  of  the  Surge  Tank 

idealization  of  the  physical  system  is  portrayed 
in  Figure  II-3-27. 


Reservoir 


Surge 

Tank 


Turbine 


Figure  II-3-27 

The  equivalent  circuit  as  seen  from  the  turbine  is  shown  in 

Figure  II-3-28. 
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Figure  II-3-28 
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The  surge  tank  supplies  or  absorbs  water  as  required  by  the 
flow  rate  variations  and  is  represented  as  a  capacitor.  The 
pipe  friction  and  throttling  losses  dissipate  energy  and  are 
represented  by  resistances.  These  losses  are  significant 
as  they  are  the  only  factor  acting  to  damp  out  oscillations 
existing  in  the  section  of  conduit  between  the  reservoir  and 
surge  tank. 

Cs  capacitor  representing  the  surge  tank 
Rt  resistance  representing  the  throttling  losses 

Rf  resistance  representing  the  pipe  friction  losses 

L ' ^  inductance  representing  the  water  inertia  over  distance  L1 

L'2  indictance  representing  the  water  inertia  over  distance  L2 

The  (15)  surge  tank  is  approximately  85%  effective  in  isolating 
the  section  of  conduit  between  the  reservoir  and  surge  tank 
from  the  turbine.  Provided  that  oscillations  or  other  undesirable 
effects  are  not  introduced,  the  surge  tank  may  be  assumed  to  act 
as  the  reservoir  and  the  relationship  between  the  gate  position 
and  output  power  is  given  by  equation  II-3-72. 

Due  to  the  large  capacitance  placed  in  the  circuit  by 
the  surge  tank,  a  very  low  frequency  oscillation  mode  is 
induced  between  the  surge  tank  and  the  reservoir.  The  effect 
of  this  oscillation  is  investigated  in  Section  II-3-11  to 
determine  if  one  is  justified  in  assuming  that  the  surge 
tank  acts  as  a  reservoir. 

Section  II-3-11  The  Hydro  Equations  for  Systems  with  a 

Surge  Tank 

The  purpose  of  the  surge  tank  is  to  reduce  the  effect  of 
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the  water  inertia  as  seen  by  the  turbine „  However,  in 
accomplishing  this,  it  introduces  a  low  frequency  oscillation 
ini_o  the  system.  Therefore,  before  the  surge  tank  can  be 
assumed  to  act  as  the  reservoir  it  is  necessary  to  determine 
the  effect  of  this  oscillation  on  the  turbine  output  power. 
The  effect  of  this  oscillation  is  investigated  by  making 
the  following  assumptions:  the  surge  tank  is  located  beside 
the  turbine,  the  throttling  and  pipe  friction  losses  are 
directly  proportional  to  the  flow  rate  in  the  surge  tank 
and  conduit  respectively Q 


Figure  II-3-29 

The  equations  relating  the  pressure  and  flow  rate  variations 
may  be  determined  by  referring  to  Figure  II-3-29. 

As  ^  aHt  =TQS  +  K2Qs 
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(II-3-73) 
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(O  L  S  A.Q  +  A  Q 


ApTJAHt 


(II-3-74) 


AQ, 


(II-3-63) 


^Qp  =  Qs  +  AQt  (II-3-75) 

Substituting  equations  II-3-73  and  II-3-74  into  II-3-75 

=  7S  Ap  A H t  (II-3-76) 
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(II-3-77) 


Substituting  equation  II—3-—77  into  IX-3-63  and  solving  for 


£Qt  AG 

as  a  function  of 
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Substituting  equation  II-3-77  into  II-3-78 


AH 


H, 


GoHt 

‘  lk2  2 

s  +  s 

"  L  K-,K2  ' 

+ 

K,  I  Q 
+  1  ° 

G  Ho 

LsApTS 

gAp kplS2. 

An  H 

P  J  o 

AG 

G 


o 


(II-3-79) 


Substituting  equations  II-3-78  and  II-3-79  into  equation 
II- 3-68 . 

L 


Ht  As 


G  Q  K,, 
O  2 


S2  + 


G  H 
O  o 


Hj 


H 


GKO 

1  o 


Go^apho 


A  P 

~ 


y 

\ 

Ft 

/K  2 

K  A  \ 

,  1  s\ 

G 

/L 

K 1K2  \ 

0 

a 

\ 

l_ 

i  _ 

+ 

_  + 

'Ho 

VS 

IS  A  / 
P 

Go 

\  gA 

p 

Ap-S2/ 

H 

O  /> 

D  /H( 

2G  i  H-. 
o  \  t, 


V- 


A  G 
G, 


(II-3-80) 

Equation  II-3-80  may  be  simplified  by  approximating  the 
s  and  s°  coefficients  as  follows. 
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The 


numerator 


coe  f f icient 


G  QoL  K2  =  \  /Ht 
^09^  TJ  ^ 


As  G  <  G0,  and  assuming  that 
Then 


G  QCL  K2  L 

c:oHo  ^AP  "CT  >  g 
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Now 


K 


2 


=  Throttling  losses  in  ft. 


The  maximum  throttling  losses  are  in  the  neighbourhood  of 
40  ft.  (reference  10,  pp.  237). 

K2  Us  =  40  ft.  and  Ug  max.  =  (D  /Ds)2  U  max. 

ir  P 


O 


Assuming  U  max.  =  20' /sec.,  then  K2/5  max.  =  10 

Therefore,  if  the  length  of  the  conduit  is  greater  than  322 
2 

ft.  the  s  coefficient  is  always  positive.  Under  normal 
operating  conditions,  Us  is  small,  therefore 
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Assuming  Ht  =  HQ,  (G/Ap)Uc  =  Up ,  Dp  =  8',  Up  max  =  20!/sec., 


ff  =  0.014,  Ap  = 


=  G. 
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Therefore,  if  the  head  H0  is  greater  than  22  ft.  and  L/Ho 
70142,  the  s^  term  will  be  positive. 
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The  s  numerator  coefficient 

Assuming  Ht  =  Ho,  (G/Ap)Qc  =  Qp,  ff  =  0.014,  Dp 
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If  L/Ho<  100,  then 


G  K  Q 
1  o 

G0^p*S  Hc 


approximating  the  and  s°  terms  of  the  numerator  and 

denominator,  equation  II-3-80  is  reduced  to  equation  II-3-81. 

It  has  been  shown  that  the  terms  (G/Gc)  (Q  /Hc)  (L/gA  )  (K  /-$  ) 

p  A 

and  (G/Gq)  (K^/Ap~^  )  (Q0/H0)  may  be  neglected  with  respect  to 

the  terms  (H^/H^)  (LAg/gAp)  and  It  can  also  be  shown 

that  the  former  terms  may  be  neglected  in  the  denominator* 
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The  coefficients  of  the  s^  terms  are  very  small  and 
therefore,  two  poles  and  two  zeros  are  introduced  at 
/g^p/LAg «  The  poles  and  zeros  cancel  each  other,  thus 
the  effect  of  the  low  frequency  oscillation  between  the 
surge  tank  and  reservoir  may  be  emitted* 

The  surge  tank  may  be  assumed  to  act  as  the  reservoir, 
with  equation  II-3-72,  giving  the  relationship  between  the 
gate  position  and  the  output  power„ 


Section  II  3-12  The  Effect  of  the  Higher  Harmonics  in 

the  Pressure  Variation 

The  linearization  of  equation  II-3-2  implies  that  the 
harmonic  content  in  the  pressure  variation  should  decrease 
as  smaller  gate  movements  are  considered*  This  was  checked 
by  graphical  means  and  was  not  found  to  be  the  case*  In 
this  section,  the  harmonic  components  will  be  assumed  to 

and  their  effect  on  system  stability  will  be  considered. 

The  system  is  considered  linear  with  a  third  harmonic 
signal  being  applied  to  represent  the  harmonic  components 
generated  by  the  non-linearity  (refer  to  Figure  II-3-30) „ 
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Figure  II-3-30 


The  pressure  variations  exhibit  half  wave  symmetry 
(Figure  II-3-22) ,  indicating  that  it’s  constituents  are 
mainly  a  fundamental  and  odd  harmonics.  Only  the  third 
harmonic  is  considered  as  it  is  the  largest  component. 

Assuming  the  power  system  does  not  greatly  attentuate 
the  third  harmonic  component  relative  to  the  fundamental, 
the  third  harmonic  will  appear  in  the  frequency  deviation 
(  Af) .  It  will  not  appear,  however,  in  the  gate  movement 
as  the  governor  is  not  capable  of  responding  to  a  frequency 
oscillation  of  0.3  cycles  per  second  (reference  10,  page  94). 
The  oscillation  will  not  have  a  large  effect  on  the  power 
system  as  it's  natural  .frequency  is  approximately  1.0 
cycles  per  second.  (18) 

In  view  of  the  above  and  the  fact  that  the  linearization 
appears  sound  mathematically,  the  harmonic  components  in 
the  pressure  variation  are  ignored. 

Section  II-3-13  The  Effect  of  the  Turbine  Wicket  Gate 

Non-linearity 

The  turbine  gate  position  is  a  non-linear  function  of 
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the  servo-motor  valve  position  as  shown  by  Figure  II-3-310  (10) 


A  step  load  change  of  50%  rated  output  power  was 
applied  to  the  governoring  loop  with  and  without  the  non¬ 
linearity.  The  resulting  oscillations  are  shown  in  Figures 
II-3-32  and  II-3-33. 


As  both  oscillations  are  damped  at  almost  an  equal 
rate,  the  non-linearity  is  seen  to  have  little  effect  on 
the  stability  of  the  loop  and  is  neglected* 

Section  II-3-14  Conclusion 

It  is  found  that  the  compressibility  of  the  water  and 
the  elasticity  of  the  conduit  do  not  effect  the  turbine 
power  output,  and  that  the  surge  tank  isolates  the  turbine 
from  the  reservoir. 


. 
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Section  II-4  The  Load  Model 

Section  II-4-0  Introduction 

To  represent  the  power  system  load  it  is  necessary  to 

determine  the  equations  relating  voltage ,  impedance  and  current. 

To  obtain  these  relationships,  the  generator  is  considered 

as  a  voltage  source  with  an  internal  impedance  looking  into 

a  network.  The  impedances  of  the  network  are  assumed  constant 

as  the  load  damping  action  is  accounted  for  in  equation  II-2-9. 

The  internal  impedance  of  the  salient  pole  generator 

is  a  non-linear  function  of  the  machine's  internal  and 

synchronizing  torque  angle  as  the  direct  and  quadrature 

reactances  are  differento  In  the  interest  of  simplicity,  it 

is  desirable  to  represent  this  impedance  as  a  constant. 

Three  alternatives  are  available:  assume  a  wound  rotor 

machine  with  an  impedance  Xd  or  X  ,  or  obtain  the  voltage 

u  q 

behind  the  quadrature  reactance  and  assume  an  impedance  Xa. 

Si 

The  latter  possibility,  although  not  perfect,  is  the  most 
accurate  and  is  utilized  in  this  study. 

The  generator  representation  is  given  and  the  errors 
are  discussed  in  Section  II-4-1.  The  equations  are  then 
developed  in  Section  II-4-2.  The  current  equations, 
along  with  equations  II-l-Gfe,  II-1-70  and  II-1-73  allow 
the  generator  output  or  load  power  to  be  determined. 

Section  II-4-1  The  Generator  Representation 

To  represent  generator  saliency,  one  must  consider 
the  voltage  behind  the  quadrature  reactance  (refer  to 
Figure  II-4-1) .  (6) 
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Figure  II-4-1 

The  above  representation  yields  a  correct  equivalent  circuit 
if  the  generator  and  its  load  are  considered  to  be  isolated 
from  the  system0  Therefore  the  non-linearity  of  the, internal 
reactance  as  a  function  of  the  internal  torque  angle  is 
eliminated,.  However f  if  two  machines  are  connected  together 
through  a  mutual  load,  the  non-linearity  of  the  impedance  of 
the  second  machine  as  seen  by  the  first  is  in  error.  There¬ 
fore  the  non-linearity  of  the  internal  impedance  as  a  function 
of  the  synchronizing  torque  angle  still  exists®  This  may 
be  seen  as  follows  (refer  to  Figures  II-4-2,  II-4-3,  II-4-4.) 


=E'-in  (Xi-y) 

f  d  o,  q' 


jXql  +  Slext  jXq2  +  S2ext 


Figure  II-4-2 
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Applying  the  theorem  of  superposition  (refer  to  Figure  II-4-3) 


A  B 


Figure  II-4-3 

The  vector  diagram  for  Figure  II-4-3-A  is  shown  by  Figure 
H-4-4,  assuming  =  jxx,  %2  =  jX2  and  83  =  Fa 


E 


I12q2  is  the  current  generated  by  unit  #1  along  the 
quadrature  axis  of  unit  #2 

1 ] 2d2  is  the  current  generated  by  unit  #1  along  the 
direct  axis  of  unit  #2* 

As  the  reactance  associated  with  I-,.,-  should  be  X-, 

12d2  d 

and  not  X  ,  this  representation  decreases  the  impedance  and 
increases  the  synchronizing  power  between  units „ 
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The  magnitude  of  the  error  depends  on  the  values  of 
3  and  8  and  the  difference  between  X.  and  X  .  The  difference 

u.  CJ 

between  Xd  and  Xq  becomes  less  significant  as  and  8~2  increase. 
Also,  the  difference  between  X^  and  X^  decreases  with  machine 
saturation.  These  factors  reduce  the  model's  error  to  very 
low  proportions. 


Section  II-4-2  The  Current  Equations 

Referring  to  Figure  II-4-3,  equations  II-4-1  to  II-4-4 
are  written,  in  terms  of  vectors,  by  inspection. 
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Through  the  use  of  a  resolver,  it  is  possible  to 
represent  E  and  E  as  vectors.  This  makes  it  possible 

qi  q2 

to  obtain  the  direct  and  quadrature  currents  correctly  and 
also  allows  a  vector  notation  to  be  attached  to  them.  However, 
to  obtain  the  power  output  it  is  necessary  to  multiply  these 
currents  by  the  direct  and  quadrature  voltages,  which  are 
represented  as  scalars  (refer  to  equation  II-1-73) *  To 
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allow  this  multiplication  to  be  carried  out  correctly,  zero 
angle  must  be  attached  to  the  currents.  This  is  accomplished 
by  considering  each  generator  to  be  at  zero  reference  and 
the  other  generator  to  be  leading  or  lagging  it.  This  has  no 
deterimental  effects  on  the  simulation  as  it  yields  correct 
generator  output  power,  which  in  turn  means  correct  frequency 
and  synchronizing  torque  angle.  On  this  basis  equations 
11“ 4- 8  and  II-4-9  are  rewritten  as  equations  II-4-10  and 

ii”4-ll,  assuming  generator  #1  to  be  leading  generator  #2 
by  an  angle  £>  . 
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The  direct  and  quadrature  axis  components 
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( II-4-16 ) 


As  the  currents  do  not  change  instantaneously  when  a 
load  disturbance  is  applied  to  the  system,  equations  II-4-13, 
II-4-14,  II-4-16  and  II-4-17  only  represent  the  currents 
under  steady  state  operating  conditions.  To  represent  this 
transient  condition,  the  direct  axis  time  constant  is  assumed 
to  be  equal  to  the  short  circuit  time  constant  and  the  quad¬ 
rature  axis  time  constant  is  assumed  to  be  zero.  This  is 
only  considered  to  be  an  approximation  as  the  time  constants 
associated  with  the  direct  and  quadrature  axis  vary  with  the 
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It  is  necessary  to  divide  the  equations  for  the  direct 
axis  currents  by  -j  to  remove  the  vector  notation  and  place 
them  on  the  same  basis  as  the  direct  axis  voltages. 

Rewriting  equations  II-4-13,  II-4-14,  II-4-16  and 
in  per  unit  and  including  the  above  mentioned 
modifications,  one  obtains  the  following  equations: 
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The  synchronizing  torque  angle  is  obtained  from  the 
frequency  deviation  which  is  expressed  in  per  unit.  As  there 
are  360fQ  degrees  per  second  squared  in  one  per  unit  cycle 
per  second,  the  individual  machine's  datum  angle  is  given 
by  equation  II-4-20. 
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The  torque  angle  between  machines 
5  =  0X  -  02  =  360  (q  -  f2) 


(II-4-21) 
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The  simulated  equations  are  II-4-14,  II-1-16,  II-1-18 
and  II-1-19  expressed  in  per  unit,  and  II-4-21  expressed  in 
degrees  . 

Section  II- 4-3  Conclusion 

The  generator  may  be  accurately  represented  by  considering 

it  as  a  voltage  source  (E„)  with  an  internal  impedance  X  . 

q 

Section  II-5  The  Controller  Model 
Section  II-5-0  Introduction 

To  simulate  the  controller,  it  is  necessary  to  determine 
the  relationships  between  the  variables  of  its  operation. 

This  is  most  easily  accomplished  by  understanding  how  the 
unit  operates,  allowing  approximations  to  be  made  and  then 
developing  the  relationships.  Therefore  a  description  of 
the  controller  is  given  in  Section  II-5-1,  approximations 
are  made  in  Section  II-5-2  and  a  method  of  simulation  is 
decided  upon  in  Section  II-5-3. 

Section  II-5-1  The  Description  of  the  Controller's  Operation 
The  functional  purpose  of  the  controller  is  to  provide 
correct  time  and  load  sharing  between  units  with  correct 
frequency  being  a  by-product.  The  operation  of  the  controller 
is  shown  in  Figure  II-5-1. 

As  the  load  sharing  feature  is  introduced  by  the  use 
of  a  bias  time  signal  and  does  not  directly  effect  the 
controller's  operation,  it  is  ignored  at  this  point  to 
facilitate  the  operational  description  of  the  unit. 
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Power 


Figure  II-5-1 

The  master  clock,  located  at  the  central  control  station, 
sends  out  a  time  signal  of  one  second  duration  every  ten 
seconds.  The  front  and  trailing  edges  of  this  signal  are 
compared  to  the  electrical  time,  thereby  determining  whether 
the  latter  is  fast  or  slow  respectively.  If,  at  the  nth 
sampling  instant*  the  electrical  time  is  in  error,  two 
signals  are  applied  to  the  governor.  The  first  adjusts 
the  permanent  speed  droop  curve  by  changing  the  speeder 
spring  setting  and  the  second  is  applied  to  the  turbine  wicket 
gates  through  the  compensating  dashpot.  When  the  time  is 
correct  at  the  (n  +  1)  sampling  instant,  the  dashpot 
signal  is  removed,  however,  the  speeder  spring  setting  is 
retSifted.  The  applied  signals  are  of  constant  magnitude 
and  not  proportional  to  the  time  error. 

The  purpose  of  these  signals  is  most  easily  seen  by 
considering  the  effect  of  the  first  without  the  second.  A 
time  error  necessitates  the  past  or  present  existence  of 
a  frequency  error,  as  time  is  the  integral  of  frequency. 
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The  purpose  of  the  speeder  spring  adjustment  is  to  correct 
the  frequency  error  and  is  therefore  a  permanent  rather 
than  a  temporary  signal 0  If  no  other  compensation  is  avail¬ 
able,  it  is  necessary  to  over-regulate  the  frequency  to 
eliminate  the  existing  time  error.  In  doing  this,  one  is 
left  with  a  frequency  error  of  opposite  sign,  which  in  turn 
produces  a  time  error.  It  is  therefore  necessary  to  provide 
a  means  of  compensation  to  counteract  this  hunting  effect. 

This  compensation  is  provided  by  the  second  signal  which  is 
applied  over  a  thirteen  second  interval.  The  first  signal 
adjusts  the  speed  droop  curve  to  correct  the  frequency  error, 
while  the  second  causes  the  unit  to  overspeed  to  correct  the 
time  error. 

The  comparison  between  the  master  clock  and  electrical 
time  is  made  through  the  use  of  a  synchronous  motor,  a 
system  of  relays  and  mechanical  linkages.  The  synchronous 
motor  rotates  an  arm  in  the  control  unit  once  every  ten 
seconds.  Assuming  correct  time,  the  arm  will  pass  over  the 
top  position  when  the  signal  from  the  master  clock  -is  received. 
If  the  time  is  slow  the  arm  will  not  have  reached  the  top 
position  and  vice  versa,  when  the  clock  signal  is  received. 

Any  error  is  sensed  through  the  system  of  relays  and  governor 
action  ensues.  It  is  to  be  noted  that  the  relay  contacts 
are  mounted  on  a  quadrant  which  is  not  rigidly  attached  to 
the  motor  frame  (refer  to  Figure  II-5-2) .  This  is  of  prime 
importance  to  the  load  sharing  feature  and  is  utilized  to 
provide  a  degree  of  stabilization. 
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Figure  II-5-2 

Load  sharing  between  units 

The  generator  output  is  measured  by  a  Kelvin  balance 
whose  force  is  counterbalanced  against  a  spring  which  is 
indirectly  connected  to  the  quadrant.  If  the  forces  are  not 
equal,  a  motor  is  energized  and  the  quadrant  is  rotated, 
thereby  making  the  quadrants"  position  directly  proportional 
to  the  load  being  carried  by  the  unit.  The  position  of  the 
quadrant  affects  the  time  error  as  seen  by  the  master  clock, 
as  with  an  increase  in  load  it  advances  the  electrical  time. 
If  one  machine  is  carrying  more  than  its  fair  share  of  load, 
its  time  appears  fast  with  respect  to  the  other  units.  The 
controller,  in  correcting  this  bias  time  error,  forces  the 
unit  to  drop  load  as  its  speed  droop  curve  is  lowered.  If 
it  is  desired  for  one  unit  to  carry  other  than  its  share 
of  load,  it  is  put  on  percentage  load  control  by  adjusting 
the  force  of  the  Kelvin  balance. 

The  controller  converts  the  load  sharing  between  units. 
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from  a  frequency  to  a  time  droop  basis.  This  may  be  seen 
by  assuming:  steady  state  conditions,  that  two  generating 
units  constitute  the  system. 

Load  sharing,  in  the  absence  of  the  controller,  is  in 
accordance  with  the  frequency  droop  (refer  to  Figure  II-5-3)  . 
Assume  the  system  load  equals  M. 

As  shown  by  Figure  II-5-3,  the  system  frequency  decreases 
as  the  load  increases.  The  division  of  the  load  between 
units  is  determined  by  the  relative  values  of  the  speed 
droops  (Sfl,  Sf2). 
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(II-5-1) 


M  =  +  M2 


Figure  II-5-3 

Load  sharing  with  a  controller  and  100%  load  control  (refer 
to  Figures  II-5-4  and  II-5-5) . 

The  controller  shifts  the  speed  droop  curves  in  the  vertical 
direction  to  yield  correct  frequency  in  its  attempt  to 
provide  correct  time.  As  S.^  =  ^1  =  ^2'  an<^  the  curves 

are  raised  proportionally  to  the  droop  settings. 


^f2  =  sf2  m2 
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5/8  sec0  =  Maximum  Time  Droop 


Figure  II-5-5 

Load  sharing  with  a  controller.  Unit  1  and  Unit  2  on  100% 
and  60%  load  control  respectively  (refer  to  Figures  II-5-6 
and  II-5-7) . 


Figure  II-5-6 
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M1  = 
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Figure  II-5-7 

When  the  machines  are  placed  on  percentage  load  control, 
the  droop  curves  are  raised  proportionally  to  the  speed  droop 
and  percentage  load  settings. 
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Stabilization  is  introduced  through  the  quadrant  movement. 

If  a  time  error  exists  at  the  ntn  sampling  instant, 
the  Kelvin  is  unbalanced  and  the  quadrant  is  rotated  very 
slightly  to  reduce  the  chances  of  a  time  error  existing 
the  following  sampling  instant.  When  the  time  error  is 
corrected,  the  Kelvin  is  rebalanced  and  the  quadrant  resumes 
its  position. 

The  quadrant  movement  is  made  very  slow  due  to  the 
destabilizing  effect  created  by  load  rejection.  If  the  load 
were  to  fall  and  quadrant  movement  was  rapid,  the  time  droop 
load  bias  signal  would  decrease  very  rapidly.  The  time  would 
appear  slow  and  the  controller  would  call  for  a  raise  signal 
thus  assisting  overspeed  and  possible  loss  of  stability.  The 
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Figure  II-5-8 
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optimum  rate  of  quadrant  movement  would  therefore  depend  on 
the  day  to  day  load  characteristics  of  the  individual  units, 

Mr.  Gaskell,  of  Calgary  Power  Limited,  tried  different  rates 
of  quadrant  movement.  He  found  that  the  rate  of  movement  was 
not  critical  and  allowing  the  quadrant  to  move  through  its 
time  droop  (5/8  seconds)  in  2-3/4  minutes  produced  very 
satisfactory  results.  As  the  optimum  value  not  only  varies 
from  day  to  day  but  also  from  unit  to  unit,  the  value  deter¬ 
mined  by  Mr.  Gaskell  is  utilized  in  this  study. 

Summarizing,  the  interrelation  between  the  controller, 
governor,  turbine  and  generator  may  be  seen  by  Figure  II-5-B. 

Section  II-5-2  The  Approximations  Made  to  the  Controller 

Operation 

The  speed  droop  curve  is  assumed  to  change  instantaneously 
when  a  time  error  is  sensed  by  the  relay  system.  In  practice, 
it  takes  approximately  one  second  for  the  droop  curve  to 
reach  its  final  value.  As  the  sampling  period  is  ten  seconds, 
this  assumption  is  justified. 

The  signal  applied  to  the  governor  dashpot  is  assumed 
to  have  a  10m,  rather  than  a  10m  +  3  second  duration,  for  m 
raise  or  lower  signals  made  in  succession  by  the  controller. 

The  error  incorporated  by  the  above  assumptions  becomes  small 
when  m  ^  3. 

The  stabilization  introduced  through  the  quadrant  move¬ 
ment  is  assumed  to  partially  counteract  the  instability  effect 
resulting  from  backlash  in  the  mechanical  linkages.  A  small 
amount  of  positive  feedback  (1/50  second)  is  utilized  to 
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increase  the  time  signal  as  seen  by  the  controller,  to 
simulate  the  backlash. 

The  controller  is  assumed  to  have  a  dead  zone  of  1/20 
of  a  second. 

Section  II-5-3  The  Simulation  of  the  Controller 

The  operation  of  the  controller  is  non-linear  as  was 
shown  previously  in  Section  II-5-1  and  II-5-2.  As  an 
accurate  linear  approximation  was  not  evident,  the  non- 
linearities  were  simulated  directly „  The  method  utilized 
to  simulate  these  non-linearities  is  now  given 0  The  repre¬ 
sentation  of  the  time  error  signal  as  seen  by  the  controller 
is  presented  first,  followed  by  the  representation  of  the 
controller  action  on  the  time  error  signal.  The  controller 
action  on  the  governor  is  then  converted  to  equivalent  fre¬ 
quency  error  signals. 

The  electrical  time  as  seen  by  the  controller 


T* 


r 


f 

o 


dt 
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(Load) 
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(Backlash) 


(II-5-4) 


The  time  error  as  seen  by  the  controller 
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(II-5-5) 


The  Load  Bias  Time  Signal 

Assuming  a  step  load  change  from  no-load  to  full  load, 
the  quadrant  will  move  to  its  extreme  position  in  2  3/4 
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minutes 0  The  quadrant  is  driven  by  a  motor  and  its  position 
is  directly  proportional  to  the  time  droop  (Figure  II-5-9) , 


Time 

Droop 


Figure  II-5-9 

The  quadrant  movement  is  represented  by  a  time  delay, 
with  hard  limiting  applied  after  t  >  T^  (refer  to  Figure 
II-5-10)  where  =  2  3/4  minutes. 


Time 

Droop 


Figure  II-5-10 

The  percentage  load  control  is  represented  by  increasing 
the  magnitude  of  the  generator  power  by  the  factor 


100 


%  Load  Setting 


The  Backlash  Signal 

The  backlash  signal  is  assumed  to  increase  the  time 
error  signal  as  shown  by  Figure  II-5-11. 
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1/50 


Figure  II-5-11 

The  former  accounts  for  the  effective  time  error  signal 
(  T '  )  as  seen  by  the  controller.  The  operation  of  the  con¬ 
troller  on  this  signal  is  shown  in  Figure  II-5-120 


Governor  Governor 

Dashpot  Speeder 

Spring 


Figure  II-5-12 

The  time  error  (AT')  must  be  greater  than  the  dead 
zone  for  the  controller  to  respond.  The  input  to  block  b 
is  zero  if  |£>T'|  <  dead  zone.  The  controller's  output  signals 

are  of  constant  magnitude  and  not  proportional  to  the  time 
error.  The  time  error  signal  is  converted  to  one  of  con¬ 
stant  magnitude  in  block  b.  The  outputs  signals  of  block  b 
are  sampled  every  ten  seconds  and  accumulated  in  block  c, 

Comparators  provide  the  required  circuit  components  in 
blocks  a  and  b.  The  sampling  and  accumulating  functions 
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required  the  use  of  circuitry  designed  and  built  by  Mr0  J.A. 

Ash  of  the  Electrical  Engineering  Department  at  the  University 
of  Alberta. 

The  feedback  signal  fed  from  block  b  to  block  a,  shown 
in  Figure  II-5-11,  represents  backlash. 

Controller  Output  Signals 

It  is  to  be  noted  that  the  magnitude  of  the  controller 
signals  are  very  carefully  set  and  held  constant.  This  is 
necessary  as  any  error  in  the  controller  simulation  is 
accumulated  as  the  signals  are  summed  in  an  arithmetic 
progression.  The  magnitude  of  the  signals  as  seen  by  the 
governor  may  be  varied  by  the  use  of  potentiometers. 

The  signals  fed  from  block  b  and  c  represent  the 
signals  applied  to  the  governor  dashpot  and  speeder  spring 
respectively.  From  the  viewpoint  of  governor  operation, 
both  of  these  signals  are  equivalent  to  speed  errors  as 
they  act  to  vary  the  position  of  the  pilot  valve. 

The  equivalent  speed  droop  frequency  signal  applied  by  the 
controller  is  equal  to  Ek'B  =  Ik  (II-5-6) 

The  equivalent  dashpot  frequency  signal  applied  by  the  con¬ 
troller  (refer  to  Figure  II-5-13) . 


Wicket 

Gates 


Figure  II-5-13 
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Assuming  that  the  wicket  gate  movement  results  from 
the  pilot  valve  movement  rather  than  from  the  applied  dashpot 
force,  the  pilot  valve  movement  resulting  from  this  force  is 
given  by  equation  II-5-7„ 

S',  F9  F 

v  =  g  =  __  ____  =  ___  (n-5-7) 

C  1  +  Ts  1  +  T  s 

r  r 

Section  II-5- 4  Conclusion 

The  operation  of  the  controller  was  found  to  be 
highly  non-linear®  As  an  accurate  method  of  linearization 
was  not  evident,  the  non-linearities  will  be  simulated® 
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CHAPTER  III 


Section  III-O  The  Analysis  of  Results 
Section  II-l-Q  Introduction 

The  general  observation  of  results  accentuates  two  querie 
The  first  being  the  oscillations  superimposed  on  the  frequency 
and  synchronizing  torque  angle.  The  second  is  the  large 
deviation  in  the  time  error  when  a  step  load  change  is  applied 
The  cause  of  the  f orementioned  oscillations  are  investigated 
and  discussed  in  Section  III-l-l.  The  time  error  resulting 
from  a  load  change  depends  on  the  response  of  the  controllers 
which  in  turn  depends  on  the  effective  speed  droop  and  dashpot 
frequency  signals.  The  response  of  the  controller,  for  an 
isolated  generator  and  load,  is  discussed  in  Section  111-1-2, 
from  the  viewpoint  of  minimizing  the  time  error  following  a 
load  change.  At  this  stage,  the  dashpot  frequency  and  time 
droop  bias  signals  are  grouped  together  and  are  defined  as  a 
damping  signal,  since  they  both  act  to  decrease  the  time 
error  seen  by  the  controllers  during  the  load  change. 

The  optimum  operation  of  the  controllers  directly  depends 
on  the  settings  of  their  parameters;  the  effective  speed  droop 
signals,  the  dashpot  frequency  signals,  the  time  bias  signals 
and  the  sampling  duration.  For  the  individual  controller  and 
generator,  these  parameters  are  interrelated  among  themselves 
and  are  also  related  to  the  permanent  speed  droop  setting. 
Further,  the  parameter  settings  of  the  controllers  are  related 
from  unit  to  unit  through  the  permanent  speed  droop  curves 
of  the  respective  units. 

To  simplify  the  analysis,  the  relative  optimum  and  then 

the  optimum  magnitude  of  the  effective  speed  droop  and  damping 
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signals,  assuming  a  sampling  duration  of  10  seconds,  are 
investigated  in  Section  III-1-3.  By  determining  the  optimum 
relative  magnitude  of  the  controllers’  signals,  the  problem 
is  reduced  to  finding  the  optimum  magnitude  of  the  signals 
for  a  single  controller,. 

From  the  former  analysis „  the  optimum  magnitude  of  the 
effective  speed  droop  signal  is  established,  as  related  to 
the  permanent  speed  droop  setting  and  the  percentage  load 
setting  of  the  individual  unite  Further  the  optimum  magnitude 
of  the  dashpot  frequency  signal  is  determined,  as  related  to 
the  effective  speed  droop  signal. 

The  optimum  value  and  effect  of  the  sampling  duration 
is  determined  in  Section  III-1-4,  as  related  to  the  natural 
period  of  oscillation  of  the  governoring  loop.  The  governoring 
loop  comprises  the  turbine,  the  generator  and  the  governor. 
Following  this,  the  effect  of  dead  zone  and  backlash,  which 
are  inherent  characteristics  in  any  mechanical  device,  are 
discussed  in  Section  III-1-5.  A  summation  of  the  total 
analysis  is  then  made  in  the  conclusion.  Section  III-1-6. 
Section  III-l-l  The  Sustained  Synchronizing  Torque  Angle 

and  System  Frequency  Oscillations 

The  sustained  synchronizing  torque  angle  and  system 
frequency  oscillations,  shown  in  the  computer  results,  will 
now  be  investigated.  The  generators  and  transmission  lines 
supplying  the  load  are  assumed  to  be  identical. 

It  is  to  be  noted  that  these  oscillations  were  not 
anticipated  when  the  computer  model  was  assembled  and  that, 
as  the  system  damping  was  decreased,  the  oscillations 
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increased  in  magnitude  to  the  point  where  the  system  became 
unstable o  In  an  attempt  to  determine  their  source,  the 
voltage  regulator  and  controller  models  were  isolated  from 
trie  power  system  without  noticeable  effecto  Therefore , 
these  oscillations  which  have  a  natural  frequency  of 
approximately  2  c0p0s.,  must  originate  in  the  governing 
loop  (Inner  Loop)  and/or  the  synchronizing  torque  loop 
(Outer  Loop) t  shown  in  Figure  III-l-lo  As  the  former  loop 
has  a  natural  frequency  of  approximately  0„1  c.p»s0,  it  is 
necessary  that  they  originate  with  the  latter. 


f  (S) 


Load  Power  as 
a  Function  of 


Figure  III-l-l 

For  the  system  frequency  oscillation  to  be  caused  by 
the  outer  loop,  it  is  necessary  that  the  power  absorbed  by 
the  load  vary  as  a  function  of  the  torque  angle  (  S  )  as 
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indicated  by  the  dashed  lines  in  Figure  III-l-l*  This  will 
now  be  shown  to  be  the  case  by  assuming  constant  internal 
generator  voltage  (refer  to  Figure  III-1-2.) 


Figure  III-1-2 
Letting  X  =  Xd  +  X^,  +  XL 

E1  =  E2 

The  current  absorbed  by  the  load  is  given  by  equation  III-l-l, 
IL  =  Et/(2R  +  jX)  +  Et  Z"6/(2R  +  jX)  =  2EjCOS  (  S/2)/(2R  +jX) 

(III-l-l) 

The  power  absorbed  by  the  load  is  given  by  equation  III-1-2 
and  is  shown  in  Figure  III-1-3, 


P  =  I  2R 
L  L 


4  E12R 


cos2  S/2  =  2E^2  (1+cosS) 


(4R2  +  X2)/R  (III“X-2) 


Figure  III-1-3 
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From  equation  III-1-2  and  Figure  III-1-3,  it  is  seen 
that  an  oscillation  of  the  angle  *0  will  produce  an  oscillation 
in  the  power  absorbed  by  the  load.  The  capability  of  this 
phenomena  to  produce  oscillations  in  the  system  generators 
will  now  be  investigated.  To  accomplish  this  the  power 
system  shown  in  Figure  III-1-2  will  be  converted  to  that 
shown  in  Figure  III-1-4  and  then  the  equations  relating 
the  generator  angles,  the  synchronizing  torque  angle  and  the 
turbine  torque,  will  be  developed.  Following  this  it  will 
b$  shown  that  the  system  is  stable  if  the  governor  responses, 
inertia  constants  and  damping  constants  are  equal.  Further, 
if  the  inertia  constants  are  not  equal  and  if  the  system 
damping  lies  within  prescribed  limits,  it  will  be  shown  that 
sustained  oscillations  will  result. 


Figure  III-1-4 

As  the  generator  voltages  are  assumed  constant,  the 
power  supplied  to  the  loads  Z^  and  Z is  constant,  and 
therefore  may  be  ignored  in  analyzing  the  cause  of  the 
oscillations.  Thus,  assuming  the  system  is  in  equilibrium, 
the  equations  relating  the  generator  angles  to  the  turbine 
and  load  torques  are  as  follows.  The  equations  are  expressed 
in  per  unit  quantities. 
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2h!  h 

2H2  02 


+  Pdl  01  =  aPTl 


+  Pd2  02  -  ^PT2 


Generator  Load  Torque  (III-1-3) 


Generator  Load  Torque  (III-1-4) 


The  variation  in  the  generator  load  torque,  as  a  function 
of  the  torque  angle,  will  now  be  obtained  (refer  to  Figure 
III-1-5) e  To  obtain  the  formulae,  the  Superposition  Thereom 
will  be  employed  to  obtain  the  currents  1^  and  I2  created 
by  E ^  and  E2  acting  alone „ 


180°- (9+  ) 


Figure  III-1-5 
Load  Torque  of  Generator  #1 
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cos  0  -  cos  (9  +  S  ) 


(III-1-5) 


E 
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cos  9  (1  -  cos  S  )  +  sin  9  sinSj 
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Load  Torque  of  Generator  #2 

T,  2 

cos  0  -  cos  (0  -  S )] 
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12 
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12 


cos  0  (1  -  cosS)  +  sin  9  sin  (- S)J 


(III-1-6 ) 
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Thus  equations  III-1-7  and  III-1-8  yield  the  relation¬ 
ship  of  the  generator  angles,  the  torque  angles  and  the 
turbine  torque . 


2H1^1  +  P 0^  +  sin©  sin  5 


3 


12 


=  AP^  -  E^2  cos©  (1  -  cos  & ) 


8 


12 


^2^2  +  Pd2^2  +  El^  sinQ  sin(-  8) 


8 


12 


=  £»PT2  “  Ei^  cos©  (1  -  cos  8 ) 


8 


12 


(III-1-7) 


(III-1-8) 


It  will  now  be  shown  that  the  system  is  stable  if 


Hj  —  H2  —  H,  P  —  Pd2  —  pd '  and  Pip-j^  —  PT2  ET  f  (w)  , 

by  redefining  0^  =  0 Q  +  S/2  and  02  =  0Q  "  S>/2. 


As  PT1 

and 

PT2  are  functions 

of 

w,  the 

relationships 

are  written 

as 

f ollows . 

APti  =  f1 

(  s. 

0O)  =  -F1  (  6  )  - 

G1 

(0O) 

(III-1-9 ) 

ApT2  =  f2 

(  s. 

0O)  =  -F2  (  &  )  - 

g2 

<0o> 

(III-1-10) 

Now,  as  generator  #2  sees  the  same  angle  0O  but  the  negative 
of  the  angle  &  to  that  of  generator  #1,  the  relationships 
of  to  F2  and  Gj  to  G2  are 

F1  (  8  )  =  -f2  (  i  )  (III-l-ll) 

G1  =  G2  (^o)  (III-1-12) 


Substituting  the  above  relationships  into  equations 
III-1-7  and  III-1-8,  and  after  further  manipulation  equations 
III-1-13  and  III-1-14  are  obtained. 
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cos©  (1  -  CCS  8 )  ( III-1-13) 


+  2F( 6  )  =0  ( II I- 1-14) 

73 12 

From  equations  III-1-13  and  III-1-14,  it  may  be  realized 

that  the  torque  angle  ( S  )  is  independent  of  the  datum  angle 

( 0q)  and  therefore  the  redefinition  of  0-^  and  02  to  0^  = 

S/2  +  0O  and  02  =  6/2  +  0O,  is  valid. 

From  equation  III-1-14  it  is  seen  that  any  oscillation 

in  the  torque  angle  will  damp  out  even  with  P^  =  0,  as  the 

governoring  loop  provides  damping  through  the  term  2F(S  ). 

Thus  the  forcing  function  (2E1/Z12  cos©  cos  6  )  of  equation 

I1I“1~14  becomes  a  constant  and  therefore  any  oscillation 

in  the  system  datum  angle  (0  )  will  be  damped  out  through 

o 

the  governor  response  G(0O) . 

Therefore,  with  equal  governor  responses,  inertia  and 
damping  constants,  the  synchronizing  torque  and  frequency 
oscillations,  shown  in  the  computer  results,  cannot  be 
justified. 

To  determine  the  effect  on  the  system  when  H-^  /  H2 
in  the  absence  of  system  damping,  equation  III-1-8  is 
subtracted  from  equation  III-1-7  to  obtain  equation  III-1-15. 
2H1^1  “  2H2^2  =  ~2E12  [sin©  sinS]  +  APT1  -  Apt2 


2H/)  +  2P</o  +  2G(0o)  =  2  E x2 


HS  +  +  2  E-j^  sin©  sinS 


Letting  8=0-^- 
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2H15  +  2E12 


sine  sin 6  +  APT2  -  &PT1  =  2(H2  -  H1)02 


(III-1-15) 


Assuming  a  slight  decrease  in  lodd,  equations  III-1-7 
and  III-1-8  show  that  both  generators  will  accelerate.  With 
H2  greater  than  or  visa  versa,  it  is  necessary  for  the 
generators  to  separate,  to  satisfy  equation  III-1-15.  The 
right  hand  side  of  this  equation  cannot  be  zero  when  the 
units  are  accelerating.  As  the  generators  separate,  the 
angle  6  increases  and  consequently  the  load  decreases. 

Since  the  governoring  loop  response  is  slow  with  respect 
to  that  of  the  synchronizing  torque  loop,  the  load  decreases 
more  rapidly  than  the  governor  can  cut  back  the  input  power 
from  the  turbine.  As  a  result  of  this  cumulative  action, 
the  system  becomes  unstable. 

To  make  the  system  stable,  the  energy  dissipated  by 
system  damping  must  be  equal  to  or  greater  than  the  energy 
supplied  to  the  system  by  a  decreasing  load.  The  energy 
dissipated  by  damping  and  that  supplied  by  a  decreasing 
load  are  given  by  equations  III-1-16  and  III-1-17  respectively. 
The  energy  dissipated  by  damping 

~Pd  d0x  + 

dt 


/* 


J 


Pd  ^2 
“dt 


d0. 


Letting  0^  =  02  +  & 


.  fd02 

2 

<302 

d  6 

[d  8^ 

2 

r>  1 

"L— j 

+  z. 

dt 

dt“ 

+ 

Ldt  J 

(III-1-16) 
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The  energy  supplied  by  a  decreasing  load 


E 


2 

1 


S 


12 


cos0  ( 1  -  cos  S )  dt 


(III-1-17) 


Now  as  oscillations  in  6  and  0^  are  found  to  be 
approximately  sinusoidal  in  the  computer  results,  the 
energy  supplied  to  the  system  minus  the  energy  dissipated  is 
given  by  equation  III-l-18o  The  angles  8  and  0^  are  assumed 
to  be  equal  to  A  sin(wat)  and  B  sin(wat)  respectively „  The 
equation  is  expressed  in  per  unit  quantities  with  the  magni¬ 
tudes  of  the  angles  expressed  in  radians. 
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cos©  £  1  -  cos  (A  sinwat)  J  -  waP( 
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[2B2  +  2AB  +  A2] 
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cos  w0t 

d 

_ 

dt 


(III-1-18) 


For  small  angle  variations  it  may  be  seen  from 
equation  III-1-15  that  02  ^  C 8  as  sin  (  8  )  ^  b0 
Thus  letting  0^  =  C  8  t  equation  III-1-18  is  rewritten 
below. 
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E^2  cos0  ^  1  -  cos  (A  sinwat)^J 


z  ~ 

waPd 

Wo)2 
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[2  c 


2  +  C  +  1  1  A2  cos2  w  t 
— 1  a 


dt 


From  the  coefficient  of  the  time  intergrand,  it  may  be 
realized  that,  with  (2C2  +  C  +  l)wa2  greater  than  E^2  cos©, 
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the  system  will  be  stable  and  oscillations  will  not  result. 

If  the  former  is  not  the  case  however,  it  is  not  only  un¬ 
necessary  but  also  unlikely  that  the  system  will  be  unstable. 
This  conjecture  is  supported  by  the  following  reasoning. 

The  system  damping  will  attempt  to  hold  the  system  together 
and  in  its  efforts  to  accomplish  this,  it  allows  the  governoring 
loop  more  time  to  decrease  the  turbine  power  input  to  offset 
the  decreased  load.  Further,  as  the  energy  made  available 
by  the  decreasing  load  accelerates  both  units,  provided 
or  H2  is  not  equal  to  infinity,  it  is  unnecessary  for  the  system 
damping  to  dissipate  all  the  supplied  energy  to  retain 
stability.  Further,  with  increasing  S  ,  the  energy  dissipated 
by  damping  increases  more  rapidly  than  the  energy  supplied  to 
the  system  by  decreasing  Dead.  This  latter  phenomena  is 
believed  to  produce  the  sustained  oscillations. 

Based  on  the  former  equations,  the  average  energy  supplied 
to  and  dissipated  by  the  decreased  load  and  increased  damping 
respectively  are  shown  in  Figures  III-1-6  and  III-1-7.  As 
the  average  frequency  is  constant,  the  average  energy 
supplied  to  and  absorbed  by  the  generator  inertia  is  zero 
and  therefore  this  energy  is  not  considered. 

The  figures  below  do  not  incorporate  the  stabilizing 
action  of  the  governoring  loop.  However,  it  is  felt  that  if 
the  sustained  oscillations  will  exist  in  the  absence  of  this 
action,  they  can  be  made  to  exist  with  it,  by  decreasing  the 
system  damping  (P^) *  This  reasoning  is  corroborated  by 
equation  III-1-19 ,  which  is  obtained  from  the  coefficient  of 
time  integrand  of  equation  III-1-18. 
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Figure  III-1-6  Figure  III-1-7 

Figures  III-1-6  and  III-1-7  show  that,  with 
E i2  cos©  (2  fQ) 2  <  pd  <  Ei2  cos© ( 2  fQ)2 

*12  ^Amax=lo57)  ^(2C^C+1)  Sjj  w_^  (2C^+C+1)  (III-1-19) 

the  energy  supplied  to  the  system  is  greater  than  that  dis¬ 
sipated  for  small  variations  in  the  angle  S  and  visa  versa  for 
large  variations 0  This  phenomena  leads  to  sustained  oscilla¬ 
tions,  since  when  the  oscillations  are  small,  they  will  increase 
to  a  given  magnitude  and  when  they  are  large  they  will  decrease 
to  the  same  given  magnitude „ 

It  is  to  be  noted  that  the  inertia  constants  were  intended 
to  be  equal  in  the  computer  simulation.  However,  as  the 
inertia  constants  directly  affect  the  gain  of  the  synchron¬ 
izing  torque  loop,  shown  in  Figure  III— 1—1,  any  error  made 
in  the  potentiometer  settings,  affecting  the  gain  of  this 
loop,  would  produce  the  affect  of  unequal  generator  inertia 
constant. 

In  view  of  the  above  and  on  the  basis  of  the  former  analysis, 
it  is  the  contention  of  the  author  that  these  oscillations 
have  a  justified  existence,  and  were  caused  by  unequal 
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inertia  constants,. 

Section  III-1-2  The  Time  Error  Resulting  From  a  Load  Change 
The  time  error  resulting  from  a  load  change  depends  on 
the  type  of  load  variation,,  the  magnitude  of  the  controllers' 
signals  and  the  sampling  duration „  As  the  controllers  operate 
on  a  static  rather  than  a  proportional  basis,  it  may  be 
realized  that  the  more  rapid  the  load  change  the  larger  will 
be  the  ensuing  time  error  for  a  given  magnitude  of  the  con- 
trollers’  parameters c 

Assuming  the  power  system  comprises  a  single  generator 
and  load  the  optimum  controller  parameters,  to  minimize  the 
time  error,  will  now  be  determined  for  a  step  full  load, 
a  step  quarter  load  and  two  ramp  load  changes.  The  first 
ramp  load  change  will  be  equal  to  the  maximum  response  of 
the  controller  (i,e.  the  controller  is  just  able  to  hold  the 
time  error  to  a  minimum  value).  The  second  load  change  will 
be  slower  than  the  maximum  response  of  the  unit. 

In  resume,  as  the  system  load  and  frequency  are  related 
through  the  permanent  speed  droop  of  the  governor,  a  load 
change  will  produce  a  proportional  steady-state  frequency 
error*  The  controller  integrates  this  frequency  error 
and,  acting  to  correct  the  resulting  time  error,  applies 
two  signals  to  the  governor.  The  first  is  the  effective 
speed  droop  signal  which  is  permanently  applied  to  correct 
the  frequency  error.  The  second  is  the  dashpot  frequency 
signal  which  is  applied  temporarily  to  correct  the  existing 
time  error. 

In  addition  to  these  signals,  the  load  sharing  feature 
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introduces  a  time  bias  signal  in  the  controller,  which  allows 
the  generators  to  lose  time  when  they  are  picking  up  load 
and  visa  versa  when  they  are  dropping  load.  To  the  individual 
generator,  each  of  which  has  its  own  controller,  this  signal 
acts  very  similar  to  the  dashpot  frequency  signalo  In  its 
absence,  it  would  be  necessary  to  increase  the  dashpot 
frequency  signal  to  prevent  the  units  from  losing  this  timeQ 
Thus,  to  simplify  the  analysis,  the  time  droop  and  dashpot 
frequency  signals  are  defined  as  a  damping  signal  (v) .  This 
damping  signal  is  assumed  to  act  in  the  same  manner  as  the 
dashpot  frequency  signal.  This  approximation  will  be  recon¬ 
sidered  when  the  optimum  controller  parameters  are  determined 
for  a  10  second  sampling  period  in  Section  III-1-3. 

To  facilitate  the  forementioned  investigation  let 
K  The  frequency  error  equal  to  the  permanent  speed  droop 

setting 

k  The  effective  speed  droop  signal 

v  the  damping  signal 

m  The  sampling  period 

tc  The  minimum  time  required  for  the  controller  to  raise  the 

speed  droop  from  its  full— load  to  its  no-load  position 
Step  Full  Load  Change 

It  will  now  be  shown  that  for  critically  damped  time 
and  frequency  errors,  the  dampinq  signal  (v)  must  be  related 
to  the  effective  speed  droop  signal  by  v  =  k  tc/2  m.  In 
making  this  analysis  the  response  time  of  the  governor,  the 
turbine  and  the  generator  are  assumed  to  be  negligible  in 
comparison  to  that  of  the  controller. 
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The  equations  necessary  to  determine  the  controller 

parameters  for  a  critically  damped  system,  are  developed 

as  follows  and  are  portrayed  in  Figure  III-l-80 

Approximating  the  frequency  error  (  A  f )  as  follows 

Af ,  for  0  51  t  <  m  =  K 

Af,  for  t>m  =  K-kt 

m 

the  time  error  is  given  by  equation  III-1-20 


m 

t 

t 

/ 

AT 

=  1/60 

r 

K  + 

r 

(K  -  k  t)dt  -  r 

v  dt 

o 

i 

m 

Jm 

- 

AT 

=  1/60  f Kt 

kt2  + 

km 

-  vt  +  vm 

(Ill- 

2m 

2 

To  have 

a 

critically 

damped  system  it  is 

necessary 

have  zero  time  error  when  the  frequency  error  is  nullified. 

The  minimum  time  required  for  the  controller  to  correct  the 

frequency  error  (K)  is  t  =  tQ  +  iru  Thus  the  required  damping 

signal  (v)  to  yield  a  critically  damped  system,  with  a 

given  effective  speed  droop  signal  (k) ,  may  be  found  by 

substituting  AT  =  0  when  t  =  tQ  +  m  in  equation  III-1-20. 

0  =  k  (t  +  m)  +  k  -  k  (t  +  m)  2  -  v  (t  ) 

°  2m  2tr  °  ° 

/ 

As  K  -  k  tQ ,  the  required  damping  signal  is  v  =  k  tQ 

m  _ 

2m 

The  time  error  resulting  from  this  load  change  as  well 
as;  the  effective  speed  droop  signal,  the  damping  signal  and 
the  frequency  signals  applied  to  the  governor,  are  shown  in 
Figure  III-l-80 

The  requirement  for  a  large  dashpot  frequency  signal 


. 


* ;  T 


- 


. 


■ 


Xi&npjLa  xoasupo'i'l  ioq/iasb  p'isi  &  *ro3  jn  vim  sax  ups*!  srlT 


115 


Freq  aency 


K 


m 


k 


T 


Summation  of  the 
Effective  Speed 
Droop  Signals 


time 


Frequen< 


Frequency 


K 


Frequency  Signals 
Appl Led  to  the 
Gove  trnor  by  the 
Controller 


Figure  III-1-8 


. 


116 


to  yield  a  critically  damped  system,  following  a  step 
load  change,  is  supported  by  the  computer  results.  By 
comparing  Figures  III-2-6  and  III-2-7,  it  may  be  realized 
that  the  oscillations  in  the  frequency  and  time  errors 
are  more  heavily  damped  as  the  dashpot  frequency  signal 
is  increased. 

Step  Quarter  Load  Change 

The  time  and  frequency  errors  resulting  from  a  step 
quarter  load  change  are  shown  in  Figure  III-1-9,  assuming 
that  the  controller  parameters  are  set  for  a  step  full  load 
change , 

It  is  to  be  noted  that,  in  the  interest  of  simplifying 
the  analysis  of  the  step  full  load  change,  the  controller 
was  assumed  to  raise  the  permanent  speed  droop  curve  from 
its  full-load  position  to  its  no-load  position  in  five  steps. 

In  this  analysis  it  is  assumed  that  twenty  steps  are  required 
to  exemplify  the  oscillation  in  the  time  and  frequency  errors. 
This  assumption  is  justified,  as  in  practice  approximately 
fifty  steps  are  utilized  to  provide  closer  frequency  tolerances 
on  the  system. 

It  may  be  realized  from  Figure  III-1-9  that  it  is  not 
possible  to  optimize  the  controller  parameter  settings  for  a 
step  full  and  step  quarter  load  change.  The  oscillations 
in  the  time  and  frequency  errors,  following  the  step  quarter 
load  change  are  intolerable.  The  optimum  damping  signal 
for  a  step  quarter  load  change  is  obtained  as  follows. 
Substitute  AT  =  0  when  t  =  tQ/4  +  m  in  equation  III-1-20. 
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0  =  K/4  (t  /4  +  m)  -  k/m  (to/4  +  m)  2  +  km  -  v  (t  /4) 

°  "2“  2“  ° 

/ 

v  =  kto 
8m 

Thus  from  the  calculations  made  on  the  step  load  changes, 
it  may  be  realized  that  the  controller  parameters  cannot  generally 
be  optimized  for  this  type  of  load  variation. 

Ramp  Load  Change  Equal  to  the  Maximum  Controller  Response 

A  ramp  load  change  which  is  equal  to  the  maximum  response 
of  the  controller  is  defined  as  one  that  creates  a  steady 
state  frequency  error  as  rapidly  as  the  controller  can  correct 
it.  The  ensuing  frequency  and  time  errors  for  this  type  of 
load  change  are  determined  as  follows  and  are  portrayed 
in  Figure  III-1-10. 

The  frequency  error  (  A  f )  is  approximated  by 

Af  -  N  t  for  o  <  t  <  m 

£f  «  (N  -  k/m)  t  +  N  m  -  4  -  k/2,  for  m  <  t  <  t 

The  term  k/2  is  incorporated  in  the  latter  approximation  to 
account  for  the  fact  that  the  summation  of  the  effective 
speed  droop  signals  changes  in  steps  rather  than  as  a  ramp 
function. 

The  equation  for  the  time  error  is 


AT  =  1/60 

-  m 

r  Ntdt  + 

t 

r* 

|~(N-k/m)  t  +Nm  -  k/2-vJ 

L  Jo  v- 

m 

_ 

AT  =  1/60  [  Nm2/2  +  (N-k/m)  t2/2  -  (N-k/m)  m2/2 

+  (Nm-v-k/2)t  -  (Nm-v-k/2)  nT|  (III-1-21) 

As  the  frequency  error  produced  by  the  load  change  is 
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equal  to  the  maximum  response  of  the  controller  (N  =  k/m) , 
the  time  error  may  be  held  to  a  minimum  by  letting  v  =  k/2. 
Assuming  v  =  k/2;  the  effective  speed  droop  signal,  the 
frequency  signals  applied  to  the  governor  and  the  time  error 
are  shown  in  Figure  III-1-10, 

From  the  figure,  it  is  seen  that  during  the  load  change, 
the  time  error  is  held  to  Nm2/2.  However,  after  the  load 
change  has  terminated  a  sustained  oscillation  of  magnitude 
Nitr1 ,  in  the  time  error,  results.  The  oscillation  may  be 
removed  by  making  the  controller  insensitive  to  time  errors 
less  than  or  equal  to  km/2,  which  is  equal  to  Nm2/2. 

If  a  dead  zone  equal  to  km/2  is  incorporated  in  the 
controller,  however,  the  controller  will  not  respond  at  the 
first  sampling  instant,  as  the  time  error  is  just  equal  to 
km/2,  Thus  with  v  =  k/2,  the  controller  is  unable  to 
totally  correct  the  frequency  error  of  2k,  at  the  second  sampling 
instant.  As  a  result  of  this,  a  time  error  is  accumulated 
during  the  load  change  which  will  produce  a  damped  oscillation 
following  the  load  change. 

The  maximum  rate  of  load  variation,  to  which  the  con¬ 
troller  can  respond,  is  indicated  by  equation  III-1-22. 

It  is  assumed,  of  course,  that  the  controller  is  insensitive 
to  time  errors  less  than  or  equal  to  km/2. 

Approximating  the  frequency  error  ( A  f )  by 
Af  =  Nt,  for  o  <  t  <  2m 

Af  =  2Nm  +  jjN  -  k/mj  t  -  v  -  k/2,  for  2m  <  t  <  t^ 
where 

t^  is  defined  as  the  time  when  AT  <  km/2  or  t^  2  tQ,  which- 
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ever  occurs  first.  It  is  necessary  to  limit  the  frequency 
approximation  in  this  manner  as,  if  AT  S.  km/2  at  a  sampling 
instant,  the  summation  of  the  effective  speed  droop  signal  does 
not  increase  and  the  equation  becomes  invalid. 

The  time  error  is  given  by  equation  III-1-22  . 


AT  =  1/60 
AT  =  1/60 


2m 

Ntdt 


+  P  £(N-k/m)t  +  2  Nm-v-k/2  J  dt 
2m 


2(N-k/m)nP  +  (v  +  k/2)m  -  2NnP 
+  (N-k/m) t±2  /2  +  (2Nm-v-k/2) t1 


( I II-1-22 ) 


Assuming  a  load  change  N  =  k/m,  with  v  =  k/2 ,  it  may 
be  realized  that  the  time  error  will  accumulate  until  the 
load  change  terminates.  The  controller  will  then  produce 
a  damped  oscillation  in  the  frequency  and  time  error,  as 
it  must  over  regulate  to  correct  the  residual  time  error. 
Thus  it  is  necessary  to  either  suffer  these  oscillations, 
for  load  changes  more  rapid  than  N  =  k/2m,  or  increase  the 
dead  zone  to  km  and  the  damping  signal  to  ^/  the  effective 
speed  droop  signal.  If  the  latter .method  of  operation  were 
to  be  implemented,  it  would  be  necessary  to  accept  a  system 
frequency  error  equal  to  the  average  value  of  the  dashpot 
frequency  signal.  Otherwise  a  time  error  would  not  exist 
to  activitate  the  controller. 

The  operation  of  the  controller  for  a  ramp  load  change 
(N  =  k/2m,  with  v  =  k/2)  is  shown  in  Figure  III-l-ll.  From 
the  dashed  lines  through  the  controller's  response,  it  may 
be  realized  that  the  controller  is  capable  of  correcting 
the  frequency  error  created  by  the  load  change.  Further, 
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Figure  III-l-ll 
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as  the  controller  is  able  to  adjust  its  rate  of  response  for 
slower  load  variations,  it  is  able  to  minimize  the  time  and 
frequency  errors  for  load  changes  less  rapid  than  N  =  k/2m. 

Thus  from  the  former  analysis  the  following  conclusions 
are  drawn s 

1)  that  the  controller  cannot  be  optimized  if  the  steady 
state  frequency  variation  is  more  rapid  than  the  maximum 
response  of  the  controller  (N  =  k/m) „ 

2)  that,  to  eliminate  oscillations  in  the  time  and  frequency 
errors,  it  is  necessary  to  employ  a  time  dead  zone 
equal  to  km/2 0 

3)  that,  to  minimize  the  system  frequency  errors;  the  rate 
of  load  change  should  not  exceed  one-half  the  maximum 
controller  response,  and  the  dashpot  frequency  signal 
should  be  equal  to  one-half  the  effective  speed  droop 
signal . 

The  response  rate  of  the  controller  may  be  increased  by 
increasing;  the  magnitude  of  the  effective  speed  droop  and 
damping  signals,  the  sampling  rate  or  by  employing  a  combin¬ 
ation  of  the  former.  The  optimum  magnitude  of  the  signals 
and  the  sampling  rate  will  be  discussed  in  Section  III-1-3 
and  III-1-4  respectively. 

Section  III-1-3  The  Optimum  Controller  Parameters  for  a 

10  Second  Sampling  Duration 
It  is  the  purpose  of  this  section  to  determine  the 
optimum  value  of  the  controllers'  parameters  for  a  10 
second  sampling  duration.  These  parameters  will  first  be 
determined  theoretically  and  then  will  be  verified  by  the 
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computer  results 0 

The  optimum  relative  magnitude  of  the  signals,  for  con¬ 
trollers  on  different  generators,  will  be  investigated  first 
and  then  the  optimum  magnitude  of  the  signals,  for  the 
controller  on  the  individual  generator,  will  be  determined. 
Following  this,  the  theoretical  results  will  be  compared 
to  the  computer  results » 

For  optimum  system  operation,  the  load  should  be 
correctly  divided  among  the  generating  units  with  minimum 
operation  of  the  controllers 0  As  each  generator  has  a 
separate  controller,  the  optimum  magnitude  of  the  signals 
will  vary  from  unit  to  unit,  depending  on  the  characteristics 
of  the  respective  generating  units „  It  is  expected  that  the 
larger  the  permanent  speed  droop  setting,  the  larger  should 
be  the  controller's  signals  to  correct  the  system  frequency 
error. 

It  is  assumed  that  the  governoring  loop  response  is 
rapid  with  respect  to  that  of  the  controller  and  therefore 
may  be  neglectedo  This  subject  is  investigated  in  Section 
III-l-4o  To  determine  the  optimum  relative  controller 
signals,  a  system  of  n  generators  is  assumed  and  the 
equations  are  developed  as  follows: 


Define 

Af 

the 

system  frequency  change 

=  Nt 

the 

system  load  change 

^,000  , 

the 

load  change  carried  by  each  unit 

S^.  , .  „ .  ,S  the  rate  of  frequency  change  as  related 
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to  the  load  change 

IS,  ,00.0  ,"6~n  the  percentage  load  setting  of  each  generator 
S  -p  1 ,  * . o ,S jnt  the  response  of  each  controller,  comprising 

the  summation  of  the  effective  speed  droop 
signals  and  the  damping  signal 
The  system  frequency  change  will  be  the  same  for  all 
generators,  therefore 

A  f  —  S  i  A  M  ^  S 1 1 1  ,ooooo*©ooo,  Sfn  a  Mj-^  *"*  ^  1  n  ^ 

(III-1-23) 


As  the  system  load  change  must  be  carried  by  the  generators 
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( III-1-24 ) 

It  is  desired  that  each  generator  will  absorb  a  given 
percentage  of  the  total  load  change =  Thus,  to  determine  the 
sufficient  conditions  for  this  to  be  accomplished,  assume 
that  each  unit  absorbs  a  percentage  of  some  load  AM. 
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Substituting  the  desired  load  apportionments  into  equation 
III-1-23,  the  controller  responses,  sufficient  to  achieve 

f 

these  load  divisions,  are  determined  by  equation  III-1-27. 
This  equation  is  developed  through  equation  III-1-26,, 
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Utilizing  generator  #1  as  the  reference 


N 


IS 

1  n 


Sfl  ^1  ’  Sfn  ~*n 


S11  “  Sln 


(III-1-26) 


As  the  controller  is  able  to  vary  its  response  for 
load  changes  slower  than  a  given  maximum,  it  is  only  necessary 
to  consider  the  ratio  of  the  maximum  controller  responses,  to 
the  permanent  speed  droop  curves  and  percentage  load  settings, 
to  determine  the  relative  optimum  controller  parameters. 
Equation  III-1-27,  shows  that  the  signals9  magnitude  should 
vary  directly  with  the  permanent  speed  droop  curve  and 
percentage  load  settings.  Further,  these  settings  are 
independent  of  the  size  of  the  system. 


S 


(III-1-27) 


This  hypothesis  will  now  be  corrobrated  by  the  computer 
results.  With  the  generators  on  100  per  cent  load  control, 
the  magnitude  of  the  effective  speed  droop  signals,  as 
related  to  the  permanent  speed  droop  settings  of  the  governor, 
will  now  be  considered.  When  the  effective  speed  droop 
signals  are  proportional  to  the  permanent  speed  droop 
settings  of  the  governor,  it  may  be  seen,  as  shown  in 
Figure  III-2-1,  that  following  a  step  load  change  the  load 
is  correctly  divided  between  the  generators.  However,  when 
the  former  is  not  the  case,  the  load  is  not  correctly  divided 
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following  a  step  load  change.  Thus,  as  may  be  seen  from 
Figures  III-2-2  and  III-2-3,  it  is  necessary  for  the  load 
sharing  feature  to  redisturb  the  load.  As  a  result,  the 
system  is  subjected  to  unnecessary  power  disturbances. 

The  effect  of  the  load  sharing  feature  on  the  optimum 
settings  of  the  effective  speed  droop  signals  may  be  seen 
from  Figures  III-2-4  and  III-2-5.  Figure  II-2-4  indicates 
the  power  disturbances  incurred  on  the  system,  following 
a  step  load  change,  when  the  units  are  on  100  and  60  per  cent 
load  control,  and  the  effective  speed  droop  signals  are  set 
for  100  per  cent  load  control.  It  may  be  seen  that,  following 
the  load  disturbance,  it  is  again  necessary  for  the  load 
sharing  feature  to  apportion  the  load  between  the  generators. 
Figure  I II- 2-5  shows  that  if  the  magnitude  of  the  effective 
speed  droop  signals  are  reduced  as  the  percentage  load 
setting  is  reduced,  the  units  will  pick  up  their  correct 
share  of  the  load. 

It  is  concluded  therefore,  on  the  basis  of  the  former 
analysis  and  the  computer  results,  that  equation  III-1-27 
provides  a  criterion  to  determine  optimum  relative  magnitude 
of  the  effective  speed  droop  signals.  If  the  controller's 
parameters  are  set  in  accordance  with  this  criterion,  the 
controllers  will  respond  as  one  unit.  Thus  when  the  optimum 
magnitude  of  the  signals  for  the  base  or  first  unit  is 
determined,  the  optimum  magnitude  of  the  signals  for  the 
other  units  is  known. 

The  effect  of  other  than  optimum  relative  controller 
parameters  will  now  be  considered.  It  will  be  first  assumed 
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that  all  generators  are  on  100  per  cent  load  control,  and  that 
the  effective  speed  droop  signals,  expressed  as  a  percentage 

of  the  permanent  speed  droop  settings,  are  not  constanto  As 
the  controller’s  are  capable  of  adjusting  their  response 

rates,  for  load  changes  less  rapid  than  their  maximum  response, 
the  rate  of  their  responses  will  be  equal  for  load  variations 
less  rapid  than  one-half  the  maximum  response  of  the  slowest 
unito  Thus  for  the  above  case,  the  relative  magnitude  of 
the  controller  signals  is  not  significants  Therefore,  only- 
load  variations  in  excess  of  the  f orementioned  value  will 
be  considered,,  Assuming  that  the  controller  response  of 
one  unit  is  one-half  that  of  the  other  (or  others)  and  that 
a  load  change  occurs  which  is  equal  to  the  maximum  response 
of  the  slower  unit,  this  unit  will  pick  up  its  share  of  the 
load  but  will  not  do  its  share  in  correcting  the  frequency 
error  accumulated  before  the  time  error  exceeds  the  dead 
zone.  Thus  the  maximum  rate  of  load  change  is  limited, 
by  the  controller’s  ability  to  minimize  the  time  and  frequency 
errors,  to  one-half  the  average  value  of  the  controller 
responses o  If  the  response  rate  of  any  unit  is  less  than 
one-half  that  of  the  other  unit  (or  units) ,  the  maximum 
rate  of  load  change  is  limited  by  the  controller’s  ability 
to  have  the  respective  generators  carry  their  correct  share 
of  the  load,  or  by  the  former  criteria,  whichever  is  the 
lesser.  Assuming  the  former  limitation  dictates  and  that 
the  load  change  exceeds  the  maximum  response  of  the  slower 
unit,  this  unit  will  not  absorb  its  correct  portion  of  the 
load  during  the  load  change »  Consequently  a  disturbance  is 
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applied  to  the  system  at  the  termination  of  the  load  change, 
when  the  load  sharing  feature  re-distributes  the  load» 

It  will  now  be  assumed  that  all  of  the  generators  are 

not  on  100  per  cent  load  control.  It  is  to  be  noted  that 
the  load  sharing  feature  is  able  to  apply  the  full  time 
droop  bias  signal  (5/8  seconds) ,  as  a  ramp  function,  in 
165  seconds.  Also  note  that  the  lower  the  percentage  load 
setting  on  a  particular  unit,  the  lower  is  the  units  required 
response  since  it  assumes  a  smaller  portion  of  the  load. 

Thus,  if  the  per  unit  load  variation,  which  any  unit  is 
required  to  absorb,  does  not  occur  in  less  than  165  seconds, 
little  advantage  is  to  be  gained  by  reducing  the  controller 
signals  in  accordance  with  the  percentage  load  settings  0 
Therefore,  although  equation  III-1-27  provides  the 
criteria  to  determine  the  optimum  relative  magnitude  of  the 
controller  signals,  it  is  the  opinion  of  the  author  that 
it  is  unnecessary  to  vary  the  magnitude  of  the  controller 
signals  in  accordance  with  the  percentage  load  settings. 

It  will  be  shown  in  Section  III-1-4  that  the  optimum  maximum 
controller  response  is  500  seconds.  Therefore,  on  account 
of  the  ensuing  frequency  oscillations,  it  is  highly  undesirable 
to  have  a  per  unit  load  change  on  any  unit,  occur  in  165 
seconds.  The  rate  of  per  unit  load  change  on  the  generators 
may  be  decreased  by  placing  more  units  on  the  system,  under 
the  regulation  of  the  controllers.  In  conclusion,  maximum 
performance  may  be  obtained  if  the  relative  controller  signals, 
expressed  as  a  percentage  of  the  permanent  speed  droop  settings. 


are  constant. 
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The  optimum  magnitude  of  the  effective  speed  droop 
signals  are  determined  oy  the  characteristics  of  the  power 
system  rather  than  by  that  of  the  controller.  A  limit  is 
placed  on  the  magnitude  of  the  effective  speed  droop  signal, 
as  the  larger  its  magnitude,  the  larger  the  power  disturbance 
created  when  the  controller  adjusts  the  permanent  speed  droop 
curve.  If  only  one  station  adjusts  its  speed  droop  curve 
at  a  given  instant,  it  will  assume  a  higher  percentage  of 
the  system  load*  Consequently  a  power  disturbance  is  fed  to 
the  other  stations.  The  magnitude  of  this  disturbance  is 
equal  to  the  sum  of  the  effective  speed  droop  and  dash pet 
frequency  signals,  expressed  as  a  percentage  of  the  permanent 
speed  droop  setting. 

The  larger  the  magnitude  of  the  effective  speed  droop 
signal,  however,  the  more  rapid  is  the  response  of  the  con¬ 
troller,  as  may  be  seen  by  comparing  Figures  III-2-7  and 
111-2-8.  Consequently,  to  handle  rapid  system  load  variations 
large  effective  speed  droop  signals  are  desired. 

As  the  power  oscillations  of  the  generator  stations 
vary  between  1  and  5  per  cent  of  the  station  capacity,  it 
is  felt  that  the  disturbances  imposed  by  the  controller 
should  be  limited  to  2  or  3  per  cent  of  the  station  capacity. 
To  accomplish  this  the  magnitude  of  the  effective  speed 
droop  signal  must  be  limited  to  2  per  cent  of  the  permanent 
speed  droop  setting,  assuming  a  dashpot  frequency  signal 
equal  to  one-half  the  effective  speed  droop  signal.  It 
will  be  shown  subsequently  that  this  is  the  optimum  magnitude 
of  the  latter  signal. 
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Thus,  in  the  interest  of  minimizing  the  power  disturbances 
«nd  in  maximizing  the  controller  response,  an  effective  speed 
droop  signal  equal  to  2  per  cent  of  the  permanent  speed  droop 
setting  is  considered  to  be  optimum. 

Accepting  the  above  as  the  optimum  effective  speed 
droop  signal,  the  optimum  dashpot  frequency  signal  will 
now  be  determined.  It  was  found  in  Section  III-1-2  that 
the  magnitude  of  the  damping  signal,  which  comprises  the 
dashpot  and  time  droop  basis  signals,  should  be  equal  to 
one-half  the  magnitude  of  the  effective  speed  signals,  if 
the  system  frequency  and  time  errors  were  to  be  minimized. 

Thus  to  obtain  the  optimum  magnitude  of  the  dashpot  frequency 
signal,  it  is  necessary  to  determine  the  contribution  of 
the  time  droop  bias  signal  to  the  damping  signal. 

In  resume,  when  a  load  variation  occurs,  the  Kelvin 
Balance  activates  a  reversing  motor  to  vary  the  quadrant's 
position  according  to  the  load  absorbed  by  the  respective 
generator.  If  the  load  variation  is  very  rapid  the  mechanism 
is  unable  to  follow,  as  it  takes  165  seconds  for  the  quadrant 
to  travel  from  its  no— load  to  its  full— load  position.  However, 
unit  ramp  load,  or  slower  load  changes,  occuring  on 
the  individual  generator  in  165  seconds  or  less,  the  load 
sharing  feature  is  able  to  respond  immediately.  Thus  the 
damping  provided  by  the  time  droop  bias  signal,  in  c.p.s., 
is  given  by  equation  III-1-28, 

U  =  Y  (5/8)  (60) 


(III-1-28) 
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where 

U  damping  signal  provided  by  the  time  droop  basis  signal 
in  cop0s. 

y  per  unit  load  change  of  the  individual  generator 
capacity  (  &  Mn/Mnr) 

percentage  load  setting  of  the  generator 
tx  time  required  for  the  load  to  change 

It  may  be  realized  from  equation  III-1-28,  that  the 
damping,  provided  by  the  time  droop  signal  decreases  with 
slower  and  smaller  load  changes 3  It  would  appear  that  the 
damping  would  increase  as  the  percentage  load  setting 
increases.  This  is  not  the  case,  however,  as  the  load  change 
absorbed  by  this  unit  decreases. 

The  average  damping  provided  by  the  dashpot  frequency 
signal  also  decreases  with  slower  and  small  load  changes 0 
For  a  ramp  load  change,  equal  to  one-half  the  maximum 
response  of  the  controller,  this  signal  is  applied  only 
during  one-half  the  time  interval  of  the  load  change  as  may 
be  seen  from  Figure  III-l-ll,  Thus  the  average  damping 
provided  by  the  dashpot  frequency  signal  is  given  by 
equation  III-1-29,  in  c0p0s0 
Average  Dashpot  Freqc  Sige  =  y  t 

V  (II I- 1-29) 


where 

t  the  minimum  time  required  for  the  controller  to  raise 
the  speed  droop  curve  from  its  no-load  to  its  full-load 
position  (  C  m  =  500  seconds) 
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The  damping  signal,  determined  in  Section  III-1-2,  was 
assumed  to  have  the  same  shape  as  the  dashpot  frequency  signal. 
Therefore  the  average  damping  signal,  required  to  minimize 
the  system  frequency  error,  is  equal  to  y(tQ/t1)  (k/2) . 

The  required  dashpot  frequency  signal,  to  yield  this  value 
of  damping  during  the  load  change,  is  obtained  by  equating 
the  former  values  and  is  given  by  equation  Ill-l-30o 

V  =  Sp  “  (5/8)  (60) 

V  =  Sp  -  0.075  (III-1-30) 

ToxT 

To  minimize  the  system  frequency  error  it  would  appear 
from  equation  III-1-30,  that  with  a  2  or  4  per  cent  permanent 
speed  droop  setting,  the  dashpot  frequency  signal  should  be 
made  negative.  However,  when  one  reconsiders  the  controllers' 
operation,  it  may  be  seen  that  this  is  not  the  case. 

The  signal  actuating  the  controller  is  the  actual 
time  error  minus  the  time  droop  signal.  Therefore,  for  a 
ramp  load  change  equal  to  one-half  the  maximum  response  of 
the  controller,  the  frequency  error  must  equal  and  remain  at 
0.0375  c.p.s.  Otherwise  a  time  error  will  not  exist  to 
actuate  the  controller.  If  the  dashpot  frequency  signal 
is  made  negative,  it  will  act  to  increase  the  frequency 
error  and  will  thereby  provide  an  activitating  time  error. 

The  controller  will  respond  to  reduce  the  frequency  error 
to  0.0375  c.p.s.  and,  as  the  dashpot  frequency  signal  is 
removed  when  the  load  change  terminates ,  a  frequency  error 
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will  not  remain  on  the  system,  It  is  assumed,  of  course, 
that  the  time  error  is  equal  to  zero  when  the  load  change 
terminates.  On  the  other  hand,  if  the  dashpot  frequency  s 
signal  is  made  positive,  a  frequency  error  of  0,0375  c,p,s, 
will  remain,  at  the  termination  of  the  load  variation. 

It  is  to  be  noted  that  in  either  of  the  above  cases, 
a  step  frequency  error  of  0,0375  c,p,s,  is  applied  to  the 
controller.  Therefore,  no  advantage  is  to  be  gained  through 
a  negative  dashpot  frequency  signal.  Further,  if  this 
signal  is  made  negative  and  a  time  error  should  be  incurred 
when  the  load  is  constant,  the  system  will  become  unstable 
as  the  controller  is  negatively  damped.  Therefore  the 
dashpot  frequency  signal  should  be  positive. 

The  criteria  for  establishing  the  magnitude  of  the 
dashpot  frequency  signal  are;  a  large  signal  is  desired  to 
maximize  the  damping,  a  small  signal  is  desired  to  minimize 
the  system  frequency  error  and  power  disturbances.  The 
latter  criteria  governs  as  a  large  damping  signal  is  not 
required,  provided  the  load  change  does  not  exceed  one-half 
the  maximum  controller  response.  Further,  the  magnitude 
of  the  system  frequency  error,  created  by  this  signal  is 
not  large.  In  the  interest  of  providing  some  damping  and 
in  minimizing  the  power  disturbances  to  3  per  cent  of  the 
generator  capacity,  it  is  concluded  that  the  dashpot  frequency 
signal  should  equal  one-half  the  effective  speed  droop  signal. 
It  is  possible  to  increase  the  dashpot  frequency  signal  and 
decrease  the  effective  speed  droop  signal  without  adversely 
affecting  the  power  disturbances,  however,  this  would  degrade 
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the  maximum  controller  response 0 

It  would  appear  possible  to  incorporate  another  signal 
to  counteract  the  frequency  error  incurred  through  the  time 
droop  bias  signalG  However,  if  this  signal  were  to  be  imple¬ 
mented,  it  would  cause  the  generators  to  drop  load  when  the 
load  is  increasing  and  visa  versa.  As  this  effect  is  un¬ 
desirable,  it  must  be  weighted  against  the  effect  of  a 
step  frequency  disturbance  in  the  controller. 

The  step  frequency  error  created  by  the  time  droop  bias 
signal  does  not  significantly  degrade  the  operation  of  the 
controller.  For  a  load  change  equal  to  one-half  the  maximum 
response  of  the  controller,  this  frequency  error  is  equal 
to  3  and  1  1/2  per  cent  of  a  2  and  4  per  cent  permanent  speed 
droop  setting  respectively.  As  this  disturbance  does  not 
exceed  that  imposed  by  the  controller  signals  it  is  con¬ 
sidered  acceptable. 

In  summary ,  it  is  concluded  that  the  effective  speed 
droop  and  dashpot  frequency  signals  should  be  equal  to  2 
and  1  per  cent  of  the  permanent  speed  droop  setting,  when 
the  units  are  on  100  per  cent  load  control.  It  was  found 
unnecessary,  though  possibly  preferable,  to  reduce  the 
magnitude  of  the  controller  signals  as  the  per  cent  load 
setting  is  reduced. 

Section  III-1-4  The  Optimum  Sampling  Period 

The  response  rate  of  the  controllers,  for  a  given 
magnitude  of  the  effective  speed  droop  signals,  may  be 
increased  by  decreasing  the  sampling  period.  It  will 
now  be  shown  theoretically,  however,  that  the  maximum 
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sampling  rate  (minimum  sampling  period  is  limited  by  the 
transient  response  of  the  governoring  loop. 

The  controller  response  comprises  the  load  sharing 
feature  as  well  as  the  dashpot  frequency  and  effective 
speed  droop  signals.  The  load  sharing  feature  will  be 
considered  first,  although  it  is  not  altered  by  the  sampling 
periodo  It  will  then  be  determined  if  this  feature  has 
a  detrimental  affect  on  the  governoring  loop  response , 

Following  this,  the  effect  of  the  dashpot  frequency  and 
effective  speed  droop  signals  will  be  considered.  The 
rate  of  summation  of  the  effective  speed  droop  signals 
is  directly  related  to  the  sampling  rate. 

To  determine  the  effect  of  the  controller,  on  the  tran¬ 
sient  response  of  the  governoring  loop,  an  impulse  load 
change  is  applied.  The  governoring  loop  and  controller 
responses  are  then  taken  individually.  By  comparing  these 
responses,  the  effect  of  the  controller  on  the  governoring 
loop  is  determined.  The  results  of  this  analysis  will  then 
be  compared  to  the  computer  results. 

The  load  sharing  feature  is  effected  through  the 

* 

quadrant  movement,  by  the  signal  £>TD,  as  shown  in  Figure  III-1-12 


Figure  III-1-12 


* 
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The  relationship  of  the  quadrant  movement  to  the 
generator  output  power,  for  rapidly  varying  loads ,  is 
approximately  given  by  equation  III-1-31. 

AT  =  AM  /M 

nx  nr 

+  Tds  ( III-1-31) 

As  the  time  droop  signal  subtracts  from  the  actual  time 

srror ,  the  load  sharing  feature  reduces  the  time  error  which 
activates  the  controller*  The  relationship  of  &Mn/Mpr, 

AT  and  a  Td  are  shown  in  Figure  111-1-13. 

From  observation  of  Figure  III- 1-13,  it  is  concluded 
that  the  load  sharing  feature  does  not  degrade  the  transient 
response  of  the  governoring  loop*  This  may  be  realized 
intuitively  as  the  load  sharing  feature  reduces  the  commanding 
signal,  through  which  the  controller  has  the  potential  of 
degrading  the  system  performance. 

The  effective  speed  droop  and  dashpot  frequency  signals 
are  shown  in  Figure  111-1-14,  for  the  same  impulse  load  change. 

A  "rapid  sampling  rate  or  a  short  sampling  period  is  assumed. 

From  Figure  III-l— 14,  it  rnay  be  realized  that  the  summa¬ 
tion  of  the  effective  speed  droop  signals  is  in  phase  with 
the  frequency  error  existing  in  the  governoring  loop.  As  the 
magnitudes  of  the  former  signals  are  fixed  by  the  response 
rate  of  the  controller,  this  method  of  operation  would 
result  in  sustained  frequency  oscillations.  It  is  possible, 
with  low  system  damping  and  a  sufficiently  rapid  sampling 
rate,  to  create  system  instability. 
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Time  Error 


Figure  III-1-14 
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To  avoid  the  above  mentioned  phenomena,  it  is  necessary 
to  make  the  controllers  insensitive  to  the  natural  frequency 
oscillations,  This  may  be  accomplished  in  either  of  two 
ways,  by  limiting  the  sampling  rate  or  by  increasing  the 
controller  dead  zone  such  that  it  is  insensitive  to  the 
time  errors o  The  former  method  will  receive  first 
consideration^. 

If  the  sampling  period  is  made  equal  to  or  greater 
than  the  natural  period  of  oscillation,  the  controller  is 
insensitive  to  the  natural  frequency  oscillations  of  the 
governoring  loop  as  may  be  seen  by  Figure  III-1-15,  From 
Figure  III-1-15,  it  would  appear  that  a  sampling  duration 
greater  than  the  natural  period  of  oscillation  would  allow 
the  summation  of  the  effective  speed  droop  signals  to  in¬ 
crease  unchecked 0  This  is  not  the  case,  however,  as  in 
practice  the  frequency  error  comprises  both  the  natural 
frequency  oscillation  and  the  controller  signals 0  Thus 
if  the  summation  of  the  effective  speed  droop  signals  tends 
to  increase,  it  produces  a  frequency  error  which  creates 
a  time  error »  The  time  error  in  turn,  corrects  the  summa¬ 
tion  of  the  effective  speed  droop  signals. 

By  increasing  the  dead  zone,  the  controller  can  also 
be  made  insensitive  to  these  oscillations.  The  amount  by 
which  it  would  be  necessary  to  increase  the  dead  zone 
(seconds) ,  assuming  a  sinsoidal  oscillations  of  period 
10  seconds,  is  given  by  equation  111-1-32, 

(Peak  magnitude  of  the  frequency 

oscillation  in  c,p0s0)  (0,637)  (5/60)  (III-1-32) 


- 


' 
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As  may  be  calculated  from  equation  III-1-32,  it  would 
be  necessary  to  increase  the  dead  zone  by  0,055  and  0.0275 
seconds  if  a  1  and  a  1/2  c.p.s. ,  oscillation,  respectively, 
occured  on  the  system.  Although  the  magnitude  of  the  system 
frequency  oscillations  are  dependent  on  the  characteristics 
of  the  system,  it  is  realizable  that  an  oscillation  of  the 
latter  magnitude  would  occur.  As  it  is  undesirable  to  in- 

f  i 

crease  the  dead  zone  by  these  magnitudes,  in  the  interest  of 
accurate  load  apportionment,  as  will  be  shown  in  Section 
III-1-5,  it  is  concluded  that  the  minimum  sampling  duration 
should  be  limited  to  10  seconds. 

The  computer  results  do  not  exhibit  oscillations  of 
this  nature  as  it  was  necessary  to  damp  the  govemoring 
loop  response,  above  the  values  normally  found  in  practice, 
to  limit  the  oscillations  in  the  synchronizing  torque  angle. 

As  may  be  seen  from  any  of  the  results,  the  frequency  oscil¬ 
lations  following  a  load  change  are  over-damped  and  therefore, 
the  controller  is  unable  to  create  this  type  of  oscillation. 

In  summary,  it  is  concluded  that  the  sampling  duration 
should  be  equal  to  or  greater  than  the  natural  frequency 
of  the  governoring  loop.  As  this  period  of  oscillation, 
for  a  waterwheel  generator,  is  known  to  be  approximately  10 
seconds,  a  sampling  duration  of  10  seconds  is  optimum.  (20) 
Section  III-1-5  The  Effect  of  Dead  Zone  and  Backlash 

It  was  shown  in  Section  III-1-3  that  it  is  necessary 
to  incorporate  a  dead  zone  equal  to  or  greater  than  km/2, 
to  eliminate  oscillations  in  the  time  and  frequency  errors. 
Utilizing  the  computer  results,  it  will  now  be  shown  that 
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backlash,  which  was  assumed  inherent  in  the  controller,  has 
the  potential  of  creating  sustained  time  and  frequency 
error  oscillations.  These  oscillations  may  be  removed  by 
increasing  the  dead  zone  such  that  it  is  equal  to  or 
greater  than  km/2  plus  the  backlash. 

After  it  is  established  that  oscillations  are  created 
by  backlash,  a  detailed  investigation  of  the  controller 
will  be  made  to  justify  the  physical  existence  of  backlash 
and  to  indicate  where  dead  zone  may  be  employed  to  eliminate 
its  affect.  Following  this  the  detrimental  effects  of  the 
required  dead  zone,  on  the  time  keeping  capability  and  load 
sharing  feature  of  the  controller,  will  be  determined. 

When  the  controller  model  was  simulated,  a  backlash 
signal  of  0.02  seconds  was  assumed,  which  affected  .the 
operation  of  the  controller,  in  the  absence  of  dead  zone, 
as  follows.  With  a  negative  time  error  the  backlash  in¬ 
creases  the  signal  activating  the  controller  by  0.02 
seconds.  When  the  controller  corrects  the  latter  signal, 
the  sign  of  the  former  is  reversed  making  the  activating 
signal  0.02  seconds  positive. 

As  backlash  introduces  a  disturbance,  it  may  be  con¬ 
sidered  as  an  energy  source  in  the  controller.  In  the 
absence  of  dead  zone  this  phenomena  produces  sustained 
oscillations  in  the  time  and  frequency  errors,  as  shown 
in  Figures  III-2-11  to  III-2-13  of  the  computer  results. 

The  effective  speed  droop  signal  is  unable  to  correct 
a  time  error  either  existing  or  introduced  in  the  system. 
Rather,  its  purpose-  is  to  correct  a  system  frequency  error. 
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As  the  backlash  is  a  time  error  rather  than  a  frequency 
error,  it  must  be  counteracted  by  the  dashpot  frequency 
signalo  Further,  the  magnitude  of  the  oscillations  in 
the  time  and  frequency  errors  will  decrease  from  a  large 
value,  or  increase  from  a  small  value,  until  the  time  error 
dissipated  by  the  dashpot  frequency  signal  is  equal  to 
that  introduced  by  back  lash » 

The  dashpot  frequency  signal,  existing  in  the  controller, 
was  0o042  and  0.082  c0p0s0  for  the  lead  changes  shown  in 
Figure  111-2-11  and  Figures  III-2-12  and  III-2-13  respec¬ 
tive  ly o  As  this  signal  is  applied  continuously  during 
any  half-cycle  period  of  the  oscillation,  the  time  in 
seconds  required  to  dissipate  the  backlash  is  given  by 
equation  III-l-33o 

t  =  (Bl)  (60)  (III-1-33) 

V 

Thus  for  the  backlash  simulated  in  the  computer  model 
(0»02  seconds) ,  approximately  30  and  15  seconds  are  required 
to  dissipate  the  energy  introduced  by  the  respective  dashpot 
frequency  signals 0  The  former  period  of  oscillation  (60 
seconds) ,  which  is  incurred  with  a  10  second  sampling  duration, 
agrees  with  the  average  period  of  oscillation  found  in  the 
computer  results 0  However,  the  computer  results  show  that 
the  latter  period  of  oscillation  is  40  rather  than  30  seconds. 
This  difference  is  believed  to  be  caused  by  two  factors. 

The  dashpot  frequency  signal  is  applied  through  the  governor 
dashpot,  therefore,  on  account  of  the  dashpot  relaxation 
time  constant  (Tr) ,  its  affect  is  neither  applied  to 
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or  removed  from  the  system  when  the  controller  applies 
the  dashpot  frequency  signal  to,  or  removes  it  from,  the 
governor o  Second,  as  the  governoring  loop  response  is  over- 
damped,  the  system  frequency  does  not  follow  (instantaneously) 
the  summation  of  the  effective  speed  droop  signals 0  Thus  a 
longer  period  of  time  is  taken  to  counteract  the  backlash 
signalo  Both  of  these  factors  have  a  larger  affect  with 
the  more  rapid  sampling  rates;  in  the  former  case  the  time 
taken  to  apply  or  remove  the  dashpot  frequency  signal  is 
greater  with  respect  to  the  sampling  duration,  in  the  latter 
case  the  response  rate  of  the  controller  increases  with  the 
more  rapid  sampling  rateQ 

If  the  dead  zone  is  made  equal  to  or  greater  than  km/2 
plus  the  backlash,  the  time  and  frequency  error  oscillations 
are  eliminatedo  This  may  be  seen  by  observation  of  Figures 
III-2-1  to  III-2-10  of  the  computer  results,  where  the  dead 
zone  is  equal  to  0.05  seconds „ 

To  determine  where  backlash  may  originate  in  the  con¬ 
troller  and  how  its  affect  may  be  eliminated  by  dead  zone, 
the  controller  is  sectionalized  as  shown  in  Figure  III-1-16. 
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Section  A-A  incorporates  the  integrator  and  the  load 
sharing  feature  of  the  controller  The  integrator  is  a  6  . 
RoPoMo  telechron  synchronous  motor  while  the  load  sharing 
feature  is  composed  of?  the  Kelvin  Balance ,  a  reversing 
motor ,  and  the  gear  train  which  drives  the  quadranto 

A  telechron  motor  runs  at  3,600  R«P0MOJ  however,  its 
output  speed  is  reduced  through  a  gear  trainc  Any  back¬ 
lash  between  the  system  frequency  and  the  rotating  arm  of 
the  controller  must  exist  in  the  gear  train  of  the  motor„ 

A  synchronous  motor  may  not  be  separated  from  the  system 
frequency  by  more  than  90  electrical  degrees,  which  is 
equivalent  to  a  time  error  of  approximately  4  x  10~  ~  seconds 
at  the  rotating  arm  of  the  controllers  Assuming  a  tolerance 
of  1/1000  in  the  gear  train  of  the  motor,  a  backlash  of  0„01 
seconds  could  be  incorporated  in  the  controller  at  this 
point,  as  it  takes  10  seconds  for  the  arm  to  make  one 
revolution  0 

The  block  schematic  of  the  load  sharing  feature  is 
shown  in  Figure  III-1-170 
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Any  play  between  the  gear  drive  of  the  load  bias  reversing 
motor  and  quadrant,  on  which  the  raise  and  lower  relay  con¬ 
tacts  are  mounted,  would  allow  backlash  in  this  section 
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of  the  apparatus 0  Assuming  2/1000  tolerance  in  the  gear 
train,  the  amount  of  backlash  incorporated  would  be  2/1000 
of  the  maximum  time  droop  bias  signal  (5/8  seconds) .  This 
is  equivalent  to  a  time  backlash  of  0.00125  seconds „ 

Section  B-B  of  the  controller  apparatus,  shown  in 
Figure  III-1-18,  will  now  be  consideredo  This  incorporates 
the  controller  action  from  the  load  bias  and  electrical 
time  signals,  to  the  command  signals  given  to  Section  C-C 
of  the  controller. 

This  section  of  the  apparatus  will  not  contain  backlash, 
as  it  is  composed  of  relays  which  are  solidly  attached  to 
the  quadrant  which  have  little  resilience  or  inertia.  It 
is  possible,  however,  to  incorporate  the  desired  insensitivity 
or  dead  zone  in  this  section  of  the  apparatus. 


Figure  III-1-18 

Assuming  that  the  electrical  time  and  load  bias  signals 
are  correct,  the  master  clock  signal  should  initiate  and 
terminate  immediately  after,  and  immediately  before,  the 
controller  opens  and  closes  the  raise  and  lower  relay 
contacts.  If  the  contacts  are  under  set,  false  time  errors 
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will  be  generated  and  if  they  are  over  set,  dead  zone  will 
be  incorporated  in  the  controller. 

Section  C-C,  which  incorporates  the  controller  operation 
from  the  comparator  to  the  point  at  which  the  signals 
are  applied  to  the  governor,  will  now  be  considered0 

The  method  by  which  the  effective  speed  droop  signal 
is  applied  to  the  governor  is  shown  in  Figure  III-1-19. 
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Figure  III-1-19 

The  backlash  in  the  above  apparatus  is  insignificant 
as  the  selsyn  motor  responds  to  the  raise  and  lower  signals 
(MCR  and  MCL) ,  applied  through  the  reversing  motor,  and 
thereby  removes  the  force  on  the  gear  train 0  The  play  in 
the  gear  train  will  decrease  the  magnitude  of  the  first 
raise  or  lower  signal  applied  to  the  governor.  This  affect 
is  not  significant,  however,  as  long  as  the  gear  tolerances 
are  small. 

The  method  by  which  the  dashpot  frequency  signal  is 
applied  to  the  governor  is  shown  in  Figure  III-l-20o 
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The  backlash  in  this  apparatus  is  also  thought  to  be 

insignificant  due  to  the  damping  in  the  dashpot  of  the 

governor.  When  a  raise  or  lower  signal  (GR  or  GL)  is 

applied  to  or  removed  from  the  governor,  the  inherent 

damping  of  the  dashpot  prohibits  any  rapid  variation  in 

the  magnitude  of  the  signal.  Thus  backlash  may  not  exist. 

From  the  former  analysis ,  it  is  concluded  that  only 

the  backlash  in  the  integrator  and  load  sharing  feature 

can  affect  the  operation  of  the  controller.  Further,  based 

on  the  assumptions  of  a  time  error  backlash  equal  to  0o01 

and  0,00125  seconds  in  the  integrator  and  load  sharing 

feature  of  the  apparatus,  the  required  dead  zone  for  optimum 

controller  operation  is  given  by  equation  III-1-34, 

Required  Dead  Zone  =  0,01  +  0,00125  +  km 

(2)  (60) 

=  0,01  +  0,00125  +  S  p 

650  (III-1-34) 

Assuming  a  2  and  4  per  cent  permanent  speed  droop  setting, 
the  required  dead  zone  is  calculated  to  be  0,01325  and  0,01525 
seconds  respectively.  For  ease  of  calculations,  the  required 
dead  zone  is  approximated  as  i0,015  seconds.  This  required 
dead  zone  (insensitivity  of  the  controller)  incorporates 
an  error  in  the  load  division  between  generating  units 
and,  of  course,  allows  a  time  error  of  *0,015  seconds  on 
the  system. 

As  the  apportionment  of  load  is  based  on  the  time 
droop  bias  signal  (5/8  seconds  for  full-load) ,  it  is  expected 
that  the  error  in  load  division  will  increase  as  the  individual 
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generator  load  decreases e  The  maximum  possible  value  of 
this  error  is  given  by  equation  III-1--35  and  is  shown  in 
Figure  III-1-210 

Possible  Error  in  Machine  Load  = 

M  M 

nr  0  o  0 30  =  4  o  8%  (III-1-35) 
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Figure  III-1-21 


As  it  was  found  necessary  to  increase  the  absolute  value 
of  the  dead  zone  to  0o030  seconds  f  to  make  the  controller 
insensitive  to  backlash,  the  controller  will  not.  respond 
to  a  frequency  error  created  by  a  load  change  until  the 
activating  time  error  exceeds  the  f ore-mentioned  value „ 

At  this  time  a  frequency  error  will  exist  on  the  system 
which  will  act  to  increase  the  time  error0  Considering 
a  ramp  load  change  within  the  maximum  response  of  the 
controller,  the  largest  frequency  error  will  be  accumulated 
for  the  most  rapid  load  change c  This  is  the  case  as  the  time 
error  is  the  integral  of  the  frequency  error. 

The  frequency  error  which  the  controller  must  correct, 
existing  when  the  activating  time  error  is  equal  to  the 
dead  zone,  will  not  significantly  degrade  the  operation  of 
the  controller  if  it  is  not  larger  than  3  or  4  per  cent  of 
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the  permanent  speed  droop  setting*  For  a  ramp  load  change 
equal  to  one-half  the  maximum  response  of  the  controller, 
assuming  that  a  positive  time  error  of  0o015  seconds  existed 


when  the  load  change  began,  the  frequency  error  on  the  system 
when  the  time  error  is  equal  to  the  dead  zone  is  given  by 


equation  III-l-360 
Absolute  Value  of  the 
Dead  Zone 


0*030 


1/2  Sp 


1 


t 


2 


Frequency  Error  (°/0)  =  t  \l 36/ SD  (III-1-36) 

10  p 

Thus  for  a  2  or  4  per  cent  permanent  speed  droop  setting, 
a  step  frequency  error  of  505  and  3*9  per  cent  respectively 
is  applied  to  the  controller  when  the  time  error  is  equal 
to  the  dead  zone  (-0*015  seconds) 0  It  is  to  be  realized 
that  these  values  represent  the  worst  case  as,  if  a  zero 
initial  time  error  were  assumed,  they  would  reduce  to  3*9 
and  207  per  cent  respectively 0 

With  a  2  per  cent  permanent  speed  droop  setting,  it  may 
be  realized  that  it  is  desirable  to  decrease  the  backlash 
below  _0«015  seconds 0  However,  as  may  be  seen  from 

III  — 1—21,  it  is  desirable  to  minimize  the  dead  zone 
regardless  of  the  permanent  speed  droop  setting  in  the 
interest  of  accurate  load  division  between  the  generating 
units o  This  latter  factor  is  thought  to  be  the  dictating 
requirement,  as  otherwise  the  controller9 s  20  per  cent 
load  setting  would  have  little  meaning,, 


. 
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In  summary f  it  is  desirable  to  minimize  the  dead  zone 
and  accordingly  the  backlash  to  provide  accurate  load  divi¬ 
sion  among  the  generating  units  and  yet,  prevent  sustained 
oscillations  on  the  systemu 

It  has  now  been  determined  that  the  controller  per¬ 
formance  may  be  optimized  for  load  changes  equal  to  or 
less  rapid  than  one-half  the  maximum  response  of  the  controller. 

Section  III-1-6  Conclusion 

The  sustained  oscillations  in  the  frequency  and  syn¬ 
chronizing  torque  angle  were  found  to  have  a  justified  existence 
if  the  generator  inertia  constants  are  unequal.  Although  this 
was  not  the  intention  in  the  computer  simulation,  it  is  ex¬ 
pected  that  slight  inaccuracies  in  the  computer  potentiometer 
settings  resulted  in  unequal  inertia  constants  which  produced 
the  oscillations o 

The  controllers  could  not  be  optimized  for  load  varia¬ 
tions  more  rapid  than  one-half  their  maximum  response.  As 
it  is  necessary  to  incorporate  a  time  dead  zone,  to  eliminate 
sustained  time  and  frequency  oscillations,  a  frequency  error 
is  accumulated  before  the  controllers  are  activated  by  a 
time  error.  To  correct  this  frequency  error  and  thereby 
prevent  an  increasing  time  error,  the  load  change  must  not 
exceed  the  former  limit.  The  maximum  response  rate  of  the 
controllers,  which  is  their  maximum  ability  to  correct  the 
frequency  error  created  by  the  load  change,  was  found  to  be 
dependent  on  the  magnitude  of  the  effective  speed  droop 
signals  and  the  sampling  rate. 
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With  all  generators  on  100  per  cent  load  control,  it 
was  found  that  the  magnitude  of  the  controller's  signals, 
as  related  to  the  respective  speed  droop  settings  of  the 
individual  generators,  should  be  constant 0  However,  it  was 
found  unnecessary,  although  possibly  preferable,  to  reduce 
the  magnitude  of  the  controller  signals,  on  a  given  unit, 
as  the  percentage  load  setting  is  reduced  on  that  unit. 

This  was  found  to  be  the  case  as  a  lower  percentage  load 
setting  reduces  the  required  controller  response. 

The  optimum  magnitude  of  the  controller's  signals  were 
found  to  be  dependent  on  the  characteristics  of  the  power 
system  since  the  controllers  apply  a  disturbance  to  the  system. 
The  power  output  of  the  generators  varies  in  accordance  with 
the  frequency  error  seen  by  the  governors.  In  turn,  the 
frequency  error  seen  by  the  governor  varies  in  accordance 
with  the  position  of  its  respective  permanent  speed  droop  curve. 
As  the  effective  speed  droop  signal  varies  the  position  of 
the  permanent  speed  droop  curve  and  the  dashpot  frequency 
signal  applies  a  frequency  error  to  the  governor,  the  con¬ 
troller  applies  a  power  disturbance  to  the  system,  via  the 
governor,  when  it  responds  to  a  time  error.  To  limit  the 
magnitude  of  this  disturbance,  it  is  the  opinion  of  the  author 
that  the  magnitude  of  the  effective  speed  droop  signal  should 
be  limited  to  2  per  cent  of  the  permanent  speed  droop  setting. 

As  it  was  found  that  the  dashpot  frequency  signal  should  be 
equal  to  one-half  the  effective  speed  droop  signal,  the  power 
disturbances  will  not  exceed  3  per  cent  of  the  station  or 


. 
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generator  capacity . 

It  was  found  that  the  sampling  period  should  be  equal 
to  or  greater  than  the  natural  frequency  of  oscillation  of 
the  governoring  loop,  otherwise  the  controllers  are  able  to 
create  sustained  frequency  oscillations  on  the  system.  This 
period  of  oscillation  is  known  to  be  approximately  10  seconds 
from  the  literature  published  (20) „  Therefore,  in  the 
interest  of  maximizing  the  controller  response,  a  10  second 
sampling  duration  is  optimum. 

As  the  controller  is  made  up  of  mechanical  elements, 
it  was  assumed  that  backlash  was  inherent  in  its  operation. 

It  was  found  that  this  backlash  has  the  potential  of  creating 
sustained  frequency  and  time  error  oscillations.  To  eliminate 
its  affect,  it  is  necessary  to  increase  the  dead  zone  to 
encompass  the  backlash.  The  dead  zone,  however,  was  found 
to  have  a  detrimental  affect  on  the  load  sharing  feature, 
as  the  latter  apportions  the  load  on  a  time  droop  basis. 

The  possible  error  in  the  load  carried  by  the  individual 
unit,  assuming  a  dead  zone  of  i0„015  seconds,  was  found  to 
be  4.8  per  cent  at  full-load  and  to  vary  inversely  with 


load. 
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Section  III-2  The  Computer  Results 

In  presenting  the  computer  results,  the  following  symbols 
are  defined: 


Af 

AT 

AT' 


T1 


T2 


the  system  frequency  expressed  in  c.p.s. 
the  time  error  on  the  system  expressed  in  seconds 
the  activating  time  error  of  unit  #1  expressed 
in  seconds 

the  summation  of  the  effective  speed  droop  signals 
of  unit  #1,  expressed  as  percentage  of  its  respective 
permanent  speed  droop  curve 

the  summation  of  the  effective  speed  droop  signals 
of  unit  #2,  expressed  as  a  percentage  of  its 
respective  permanent  speed  droop  curve 
the  turbine  output  power  of  unit  #1  expressed  in 
per  unit 

the  turbine  output  power  of  unit  #2  expressed  in 
per  unit 

the  synchronizing  torque  angle  between  machines 
expressed  in  degrees  (assumed  positive  when  unit 
#1  is  leading  unit  #2) 
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Figure  III-Z.-3  ((test 
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Figure  III  ~Z  “4-  (TEST  13) 
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Section  III-1-3  General  Conclusion 

The  generator ,  governor,  hydro-system,  load  and  controller 
models  were  simulated  on  the  computer  as  developed  in  Sections 
1  to  5  of  Chapter  II .  The  voltage  regulator  model,  developed 
in  Reference  4,  was  accepted  and  incorporated  with  the  gen¬ 
erator  model.  The  R.M.S.  magnitude  of  system  voltage,  the 
variation  from  system  frequency,  as  well  as  the  ether  quantities 
were  represented  by  proportional  D.C.  quantities.  The  system 
simulated  consisted  of  two  generating  stations,  approximately 
150  miles  apart,  with  a  load  concentrated  half-way  between 
the  stations. 

The  voltage  variations  in  the  generator  and  voltage 
regulator  models  were  not  found  to  be  adversely  affected  by 
the  operation  of  the  controller,  for  the  sampling  periods 
tested  (2  to  10  seconds).  This  result  was  expected  as  the 
natural  frequency  of  this  loop  is  approximately  0.5  to  1.0 
c.p.s.  Further  for  the  same  reason,  the  controllers  did 
not  affect  the  power  oscillations  between  the  generating  units 
(synchronizing  torque  loop) .  The  natural  period  of  this 
oscillation  was  found  to  be  from  1  to  3  seconds. 

However,  it  was  found  that  the  sampling  periods,  of 
the  individual  generators,  should  be  equal  to  or  greater  than 
the  natural  period  of  oscillation  of  the  governoring  loop 
(10  seconds)  (20).  The  governoring  loop  consists  of  the 
generator,  the  turbine  and  the  governor.  If  the  sampling 
rate  is  increased  above  this  value,  it  is  possible  for  the 
controller  to  produce  sustained  oscillations  in  the  system 
frequency.  This  may  be  realized  as  an  oscillation  in  the 
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system  frequency,  caused  by  a  power  disturbance,  creates  a  time 
error.  The  controllers  respond  to  this  time  error  and  attempt 
to  adjust  the  speed  droop  curves  to  correct  the  frequency 
error.  Thus  the  frequency  error,  seen  by  the  governor,  increases 
and,  as  a  result,  the  oscillation  in  the  turbine  output  power 
also  increases.  This  in  turn  increases  the  frequency  error 
on  the  system.  Therefore,  to  maximize  the  controller's  re¬ 
sponse  but  prevent  sustained  frequency  oscillations  on  the 
system,  the  sampling  period  should  be  equal  to  10  seconds. 

With  a  sampling  period  equal  to  or  greater  than  the 
natural  period  of  the  governoring  loop,  the  controller  does 
not  affect  the  transient  response  of  the  power  system.  How¬ 
ever,  when  the  controllers  adjust  the  permanent  speed  droop 
curves  of  a  power  station,  a  disturbance  is  applied  to  the 
system.  The  magnitude  of  this  disturbance  is  directly  pro¬ 
portional  to  the  magnitude  of  the  effective  speed  droop  and 
dashpot  frequency  signals  as  they  determine,  via  the  governor, 
the  change  in  the  turbine  output  power.  Therefore,  to  minimize 
the  system  power  disturbances,  it  is  desirable  to  limit  the 
magnitude  of  the  controller  signals.  However,  there  are  also 
requirements  for  large  effective  speed  droop  and  dashpot 
frequency  signals.  The  former  determines  the  maximum  rate 
of  the  controller  response,  for  a  10  second  sampling  duration, 
and  the  latter  determines  the  controller  damping.  In  the 
interest  of  the  foregoing  conflicting  criteria,  it  is  the 
opinion  of  the  author  that  the  effective  speed  droop  and  dash- 
pot  frequency  signals  should  be  equal  to  2  and  1  per  cent  of 
the  permanent  speed  droop  setting  respectively. 
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With  all  generators  on  100  per  cent  load  control  it  was 
found  that  the  relative  magnitude  of  the  controller  signals, 
expressed  as  a  percentage  of  respective  speed  droop  settings 
of  the  units  on  which  the  controllers  operate,  should  be 
equal.  Otherwise,  following  a  system  load  change,  the  con¬ 
trollers  with  the  smaller  signals  must  respond  more  rapidly. 

Thus  the  maximum  response  of  the  controllers  is  unnecessarily 
limited.  Further,  if  a  per  unit  ramp  load  change  does  not 
occur  in  less  than  165  seconds  on  any  generator,  it  was 
found  unnecessary  to  reduce  the  magnitude  of  the  controller 
signals  as  the  percentage  load  setting  is  reduced.  This  was 
found  to  be  the  case  as  a  lower  percentage  load  setting  re¬ 
duced,  rather  than  increased,  the  required  controller  response. 

former  restrictions  were  imposed  on  the  controller 
by  the  system  characteristics.  As  the  controller  is  also  a 
mechanical  device,  it  has  inherent  characteristics,  which 
could  be  detrimental  to  the  power  system  performance.  The  con¬ 
troller  signals  are  of  constant  magnitude  and  therefore  the 
system  time  error  can  only  be  corrected  within  the  tolerances 
of  the  signals,  thus  it  is  necessary  to  make  the  controller 
insensitive  to  time  errors  below  this  magnitude  or  accept  small 
but  sustained  time  and  frequency  oscillations.  Further, 
assuming  a  per  unit  tolerance  of  1/1,000  and  2/1,000  in  the 
gear  drives  of  the  comparator  and  load  sharing  feature 
respectively,  a  backlash  of  approximately  0.01  seconds  results 
from  the  inertia  of  the  components.  This  backlash  has  the 
potential  of  increasing  the  magnitude  of  the  former  oscillations. 
To  prevent  these  sustained  oscillations  from  occuring,  it 
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is  necessary  to  make  the  controller  insensitive  to  time  errors 
less  than  or  equal  to  approximately  +0.015  seconds.  As  these 
oscillations  are  undesirable,  it  is  assumed  that  this  dead 
zone  is  implemented. 

The  dead  zone  has  two  main  detrimental  affects.  Since 
the  load  is  divided  among  the  units  on  a  time  droop  basis, 
the  error  in  load  division  varies  inversely  proportional 
to  load,  and  is  4.8  per  cent  at  full  load.  Second,  a  frequency 
error  must  exist  on  the  system,  which  is  proportional  to  the 
dead  zone  and  rate  of  load  change,  before  a  time  error  is 
created  to  which  the  controller  can  respond.  To  correct 
this  frequency  error  in  minimum  time,  the  rate  of  load  change 
must  not  exceed  one-half  the  maximum  response  of  the  controller. 
The  maximum  response  is  defined  as  the  maximum  rate  at  which 
a  steady  state  frequency  error  can  be  corrected.  With  an 
effective  speed  droop  signal  equal  to  2  per  cent  of  the  per¬ 
manent  speed  droop  setting  and  a  sampling  duration  of  10  seconds, 
it  takes  the  controller  a  minimum  of  500  seconds  to  raise  the 
speed  droop  curve  from  its  full-load  to  its  no-load  position. 

The  speed  droop  curve  will  be  raised  as  a  ramp  function.  Thus 
the  controller  may  be  optimized  (or  can  minimize  the  time  and 
frequency  errors)  for  per  unit  ramp  load  changes,  occuring 
on  the  individual  generator,  as  rapidly  as  1,000  seconds. 

If  the  rate  of  load  change  exceeds  this  value,  in  a  manner 
which  creates  a  time  error  at  the  termination  of  the  load 
change,  damped  frequency  and  time  error  oscillations  will 
result. 
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Therefore,  to  optimize  the  controllers  performance  on 
the  system,  it  is  necessary  to  have  sufficient  generation 
under  the  regulation  of  the  controller  such  that  the  load 
variations  do  not  exceed  the  former  limit. 
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APPENDIX  A 


DEFINITION  OF  SYMBOLS 

Section  1-2 

The  symbols  utilized  in  the  method  of  system 
simulation 

Id'  Xq 

the  direct  and  quadrature  axes  currents 

Vd'  vq 

the  direct  and  quadrature  axes  voltages 

w 

the  electrical  frequency 

© 

the  generator's  internal  torque  angle 

91'  92 

the  internal  torque  angles  of  generators 
#1  and  #2 

^  12 

the  synchronizing  torque  angle  between 
generators 

Section  II-l 

The  symbols  utilized  in  the  generator  and 
voltage  regulator  models 

a 

the  turns  ratio 

Ef 

the  internal  generator  voltage 

Eq 

the  voltage  behind  the  quadrature  reactance 

Fa 

the  magnetomotive  force  created  by  the 
current  in  phase  A 

Fda 

the  component  of  Fa  along  the  direct  axis 

F 

qa 

the  component  of  Fa  along  the  quadrature  axis 

Fds 

the  stator  magnetomotive  force  resulting 
from  the  three  phase  currents  along  the 
direct  axis 

Fqs 

the  stator  magnetomotive  force  resulting 
from  the  three  phase  currents  along  the 
quadrature  axis 

^a'  '  -^c 

the  instanteous  phase  currents 

Ia'  Xbf  Ic 

the  maximum  phase  currents 
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IBr 

the  rotor  base  current 

X  a  X  -  I 
a9  q 9  xo 

the  direct,  quadrature  and  zero  sequence 
currents 

Ikd'  Ikq 

the  direct  and  quadrature  axes  currents  in 
the  damper  windings 

*f 

the  field  current 

*f  o 

the  field  current  required  to  generate  rated 
voltage  on  the  air  gap  line 

Laao 

the  total  constant  term  of  the  self  of 
phase  A 

Laa'  Lbb '  Lcc 

the  self  inductances  of  the  phases 

^ab '  ^bc'  -^ac 

the  mutual  inductances  between  phases 

Laf f 'Lbf f f 

the  mutual  inductances  between  the  rotor  and 
stator  phases 

Laf  =  Lgd 

the  direct  axis  mutual  inductance  between  the 
stator  phases  and  rotor  as  seen  by  the  latter 

Lfl 

the  field  leakage  inductance 

Lgaa 

the  self  inductance  of  phase  A  due  to  the 
air  gap  flux 

Lg° 

the  constant  term  of  L 

gaa 

Lg2 

the  amplitude  of  the  second  harmonic  of 
the  variable  term  in  L 

9" 

Lgq 

the  quadrature  axis  mutual  inductance  between 
the  stator  phases  and  the  rotor  as  seen  by 
the  latter 

Lkd'  Lkq 

the  direct  and  quadrature  axes  damper  winding 
leakage  reactances 

N 

a 

the  effective  turns  ratio 

Pgd'  Pgq 

the  permeance  coefficients  along  the  direct 
and  quadrature  axes 

PL 

the  generator  output  power 

R 

a 

the  armature  resistance 

^d'  Rkq 

the  direct  and  quadrature  axes  damper  winding 
resistances 

' 
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R. 


'u 


V 


act 
VBr 

Vd'  va-  V, 


the  field  resistance 
the  electrical  air  gap  torque 
the  generator  terminal  voltage 
the  rotor  base  voltage 

the  direct,  quadrature  and  zero  sequence 
voltages 


vf 

the 

Vref 

the 

vt 

the 

w 

O 

the 

W  =  w 

the 

&  W 

the 

Xad '  Xaq 

the 

xadr 

the 

Xads 

the 

xal 

the 

the 

xBr 

the 

xd>  xq 

the 

Xkd '  Xkq 

the 

the  mutual  reactance  referred  to  the  stator 


0 


gba 


^gda'  ^gqa 


0gd<  % 


^a,  Ab, 


Ttj  9s. 

d  t  q  t  o 


the  mutual  flux  in  phase  B,  resulting  from 
the  flux  in  phase  A 

the  air  gap  fluxes,  along  the  direct  and 
quadrature  axes,  resulting  from  the  flux 
in  phase  A 

the  direct  and  quadrature  fluxes,  resulting 
from  the  currents  in  the  three  phases 

the  per  phase  flux  linkages 

the  direct,  quadrature  and  zero  sequence 
flux  linkages 


■ 
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dr. 

^  qr 

kd. 

Tv, 

kq 

s 


the  direct  and  quadrature  axis  flux  linkages 
referred  to  the  rotor 

the  direct  and  quadrature  axes  damper  winding 
flux  linkages 

the  synchronizing  torque  angle  between  machines 
the  generator’s  internal  torque  angle 
the  rotor  position  expressed  in  radians 


Section  II  2  The  symbols  utilized  in  the  Speed  Governorinq 

Model 

B  the  effective  gain  of  the  flyball  assembly 

as  viewed  in  terms  of  the  electrical  frequency 

fe  the  electrical  frequency 


A  f. 


eo 


G 

H 

K 

K' 


the  deviation  from  the  electrical  frequency 
the  rated  electrical  frequency 
the  turbine  wicket  gate  position 
the  generator  and  turbine  inertia  constant 
the  gain  of  the  governor 

the  governor  gain  between  the  pilot  valve 
and  the  servo  motor  valve 

the  sum  of  the  self  regulation  of  the  turbine 
and  load 

the  load  power  at  rated  frequency 

the  load  power  at  rated  frequency  plus  the 
change  in  power  resulting  from  the  voltage 
regulator  frequency  relationship 

the  self  regulation  of  the  load 


the  load  power 

the  turbine  power  at  rated  frequency 


the  self  regulation  of  the  turbine 
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PT 

the  turbine  power 

R.PoMo 

the  revolutions  per  minute 

tl 

the  load  torque 

tt 

the  turbine  torque 

Tr 

the  dashpot  relaxation  time  constant 

Trr 

the  rated  machine  torque 

Wm 

m 

the  mechanical  frequency 

AW 

m 

the  deviation  from  rated  mechanical  frequency 

W 

*mo 

the  rated  mechanical  frequency 

WR2/g 

the  moment  of  inertia 

XP 

the  pilot  valve  position 

*v 

the  servo  motor  position 

8P 

the  effective  permanent  speed  droop 

S'p 

the  actual  permanent  speed  droop 

St 

the  effective  temporary  speed  droop 

S't 

the  actual  temporary  speed  droop 

Section 

II-3  The  Symbols  Utilized  in  the  Hydro  Model 

a 

the  acceleration  of  the  water  in  the  conduit 

A 

P 

the  area  of  the  conduit  or  pipe 

AAp 

the  change  in  area  of  the  conduit  or  pipe 

AS 

the  area  of  the  surge  tank 

A 

r 

the  area  of  the  reservoir 

A^jn 

the  area  of  the  turbine 

b 

• 

the  velocity  of  sound  in  water 

C 

the  equivalent  electrical  capacitance 
of  the  conduit 
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D 

P 

D 


s 


E 


1 


the  discharge  coefficient 

the  equivalent  electrical  capacitance  of 
the  surge  tank 

the  diameter  of  the  conduit  or  pipe 

the  diameter  of  the  surge  tank 

the  equivalent  voltage  to  pressure  head 
variations 

the  sending  voltage 

the  receiving  voltage 


A  f 


F 


g 

G 

AG 


H 

AH 

Ho 

H 

r 

A  H 


H 


t 

AH, 


AH 


Ln 


the  electrical  frequency  deviation 

the  coefficient  of  pipe  friction 

the  force  on  the  water  column 

the  force  of  gravity 

the  wicket  gate  position 

the  change  in  wicket  gate  position  that 
will  deliver  rated  power  to  the  turbine 

the  piezometric  head 

the  change  in  the  piezometric  head 

the  rated  head 

the  piezometric  head  at  the  reservoir 

the  change  in  the  piezometric  head  at  the 
reservoir 

the  piezometric  head  at  the  turbine 

the  change  in  the  piezometric  head  at 
the  turbine 

the  change  in  the  piezometric  head  over  a 
segiment  of  the  conduit 

the  current  which  is  equivalent  to  the  flow 
rate  in  the  conduit 


the  receiving  current 


. 
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I 


s 


K 


K 


1 


K 


2 


L 


A  L 
L 

e 


L* 


L' 


1 


L' 


2 


m 


P/15 


P 


t 

A  P 


t 


Q 


A  Q 


Qr 

A  Qr 


the  sending  current 

the  bulk  modulus  of  elasticity 

the  approximate  pipe  friction  coefficient 

the  approximate  throttling  loss  coefficient 

the  length  of  the  conduit 

the  change  in  length  of  the  conduit 

the  effective  length  of  the  conduit  as 
seen  by  the  turbine 

the  inductance  which  is  equivalent  to  the 
inertia  of  the  water  in  the  conduit 

the  length  of  the  conduit  from  the  turbine 
to  the  surge  tank 

the  inductance  equivalent  to  the  inertia  of 
the  water  from  the  turbine  to  the  surge  tank 

the  length  of  the  conduit  from  the  surge 
tank  to  the  reservoir 

the  inductance  equivalent  to  the  water  inertia 
from  the  surge  tank  to  the  reservoir 

the  line  propagation  constant 

the  pressure  energy 

the  turbine  power  at  rated  frequency 

the  change  in  turbine  power  at  rated  frequency 

the  flow  rate 

the  change  in  flow  rate 

the  flow  rate  that  yields  rated  turbine  power 

the  flow  rate  in  the  conduit 

the  change  in  flow  rate  in  the  conduit 

the  flow  rate  in  the  conduit  prior  to 
a  disturbance 

the  flow  rate  at  the  reservoir 

the  change  in  flow  rate  at  the  reservoir 
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- 
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183 


O 


wt 

*Qi 
A  Q 


t3 


Q 


Ln 


R 

R. 


R. 


s 

t' 


T 

A-T 


T' 


w 


U 


U 

P 


u 


t 


V 

p 

V 

c 


w 


the  flow  rate  in  the  surge  tank 

the  flow  rate  in  the  turbine 

the  change  in  flow  rate  in  the  turbine 

the  third  harmonic  component  of  the  flow 
rate  at  the  turbine 

the  change  in  flow  rate  over  a  segment  of 
the  conduit 

the  radius  of  the  conduit 

the  resistance  representing  the  pipe  friction 
losses 

the  resistance  representing  the  throttling 
losses 

the  Laplace  operator 

the  thickness  of  the  conduit 

the  hoop  stress  on  the  conduit 

the  change  in  hoop  stress  on  the  conduit 

the  water  starting  time 

the  water  starting  time  from  the  turbine 
to  the  surge  tank 

the  water  velocity 

the  water  velocity  which  will  yield  rated 
turbine  power 

the  water  velocity  in  the  conduit 

the  water  velocity  in  the  conduit  prior 
to  the  application  of  a  disturbance 

the  water  velocity  in  the  turbine 

the  water  velocity  in  the  surge  tank 

the  volume  of  the  conduit 

the  volume  of  the  water 

the  electrical  frequency 

the  water  compressibility 
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X 


y 

% 


gTT 
8 'ft'  3 


8'^  3 


the  distance  from  the  turbine  toward  the 
reservoir 

the  shunt  admittance 

the  potential  energy  per  unit  weight 
the  characteristic  impedance  of  the  line 
the  series  impedance 

the  series  impedance  of  the  line  to  third 
harmonic  variations 

the  shunt  impedance  of  1/2  the  line 

the  shunt  impedance  of  1/2  the  line  to 
third  harmonic  variations 

the  efficiency 

the  specific  weight  of  water 

the  natural  frequency  of  the  water  hammer 
oscillations 

the  density  of  water 


Section  II- 4  The  symbols  utilized  in  the  Load  Model 


(Note : 


I 


d' 


•I 


f 


I 


kd 


Barred  symbols  represent  vector  quantities) 

the  voltage  behind  the  quadrature  reactance 
the  direct  and  quadrature  axis  currents 
the  field  current 

the  direct  axis  current  component  in  the 
damper  windings 


I 


11 


I 


12 


1 12d2 '  1 12q2 


I 


22 


the  current  flowing  in  generator  #1  when 
is  acting  alone 

the  current  flowing  in  generator  #2  when 

E  _  is  acting  alone 
ql 

the  current  generated  by  unit  #1  along  the 
direct  and  quadrature  axes  of  unit  #2,  when 
E  ^  is  acting  alone 

the  current  flowing  in  generator  #2  when 
E^  is  acting  alone 


' 


. 
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I2 1 

the  current  flowing  in  generator  #1  when 

E^2  is  acting  alone 

glext'  Z2ext 

the  impedances  between  the  generator 
terminals  and  the  load 

ell'  022'  e12 

the  angles  associated  with  the  impedances 

«•, ,  ,  Soo  ,  and  S,  0 

11  2^  12 

0±,  02 

the  machine's  datum  angle  or  angular  position 

6 

the  synchronizing  torque  angle  between  machines 

Section  II-5 

The  symbols  utilized  in  the  Controller  Model 

B 

the  effective  gain  of  the  flyball  assembly 
in  terms  of  the  electrical  frequency 

C 

the  spring  constant  of  the  flyball  assembly 

Z> 

the  summation  of  the  effective  speed  droop 
signals 

Zk- 

the  summation  of  the  effective  speed  droop 
signals  applied  to  the  flyball  assembly 

f 

the  electrical  frequency 

6  f 

the  deviation  from  the  electrical  frequency 

fo 

the  rated  electrical  frequency 

A  fi,  A  f2 

the  deviation  from  the  no-load  to  the  load 
frequency  when  the  units  are  carrying  loads 

M-^  and  M2 

F ' 

the  force  applied  to  the  dashpot  of  the 
governor  by  the  dashpot  frequency  signal 

m 

the  sampling  duration 

M 

the  system  load 

Mi 

the  load  carried  by  unit  #1 

m2 

the  load  carried  by  unit  #2 

pt 

the  turbine  power  at  rated  frequency 

P1 

the  load  power  a  rated  frequency 

' 


V 
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q 

the  movement  of  the  governor  dashpot  in 
response  to  the  equivalent  dashpot  frequency 
signal 

Sfl'  Sf  2 

the  rate  of  change  of  frequency,  as  related 
to  the  load,  on  generators  #1  and  #2 
respectively 

Stl,  St2 

the  rate  of  change  of  the  time  droop,  as 
related  to  the  load,  on  generators  #1  and 
#2  respectively 

A  T 

the  time  error 

A  T’ 

the  activating  time  error  as  seen  by  the 
controller 

V 

the  effective  dashpot  frequency  signal 

XV 

the  servo  motor  position 

«’t 

the  temporary  speed  droop 

Section  III- 

-1  The  symbols  utilized  in  the  Analysis  of 

Results 

bl 

the  magnitude  of  the  backlash  signal  in  seconds 

El'  E2 

the  voltage  behind  the  synchronous  reactance 
of  generators  #1  and  #2  respectively,  assuming 

a  round  rotor  machine 

fo 

the  rated  system  frequency 

a  f 

the  frequency  error  on  the  system 

V  f  2 

the  functional  relationship  between  the 
turbine  output  power  of  generators  #1  and 
#2  respectively,  the  synchronizing  torque 
and  frequency  datum  angles 

Fl'  F2 

the  functional  relationship  between  the 
turbine  output  power  of  generators  #1  and 
#2  respectively  and  the  synchronizing  torque 
angle 

Gl'  G2 

the  functional  relationship  of  the  turbine 
output  power  of  generators  #1  and  #2 
respectively  and  the  frequency  datum  angle 

H 

the  generator  and  turbine  inertia  constant 

' 
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Hx,  H2 

the  generator  and  turbine  inertia  constants 
of  units  #1  and  #2  respectively 

xr  x2 

the  output  currents  of  units  #1  and  #2 
respectively 

XL 

the  current  absorbed  by  the  load 

k 

the  magnitude  of  the  effective  speed  droop 
signals 

K 

the  frequency  error  equal  to  the  permanent 
speed  droop 

m 

the  sampling  duration 

&  Miji 

the  system  load  change 

AM„ 

the  load  change  absorbed  by  the  nth  generator 

M 

nr 

the  rated  capacity  of  the  nth  generator 

N 

the  rate  of  change  of  system  frequency 
caused  by  the  load  change 

pd 

the  system  damping  coefficient 

PL 

the  power  absorbed  by  the  load  simulated 
on  the  computer  model 

A  PT 

the  turbine  output  power 

A  Pm  l  - 

T 1 1 

a  p.,_  ■Hhe  turbine  outout  oov/er  on  units  #1  and 

#2  respectively 

R 

the  system  load  resistance 

Sfn 

the  rated  change  of  frequency,  as  related 
to  load,  on  the  nth  generator 

Sln 

the  controller  response  of  the  n^h  unit 

s 

the  laplace  operator 

t 

the  time  in  seconds 

the  minimum  time  required  for  the  controller 
to  raise  the  permanent  speed  droop  curve  from 
its  full-load  to  its  no-load  position 

n 

the  time  interval  from  the  start  of  the 
load  change  to  until  the  controller  corrects 
the  time  error 

. 
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A  T 
ATd 

A  T' 
U 


V 

v* 

w 

W 

a 

y 


X 


0O 

&1*  02 
© 


the  time  error 

the  time  droop  bias  created  by  the  load 
sharing  feature 

the  activating  time  error  signal  (AT  -ATd) 

the  equivalent  damping  signal  created  by 
the  time  droop  bias 

the  dashpot  frequency  signal 

the  damping  signal  (V  +  u) 

the  system  frequency 

the  natural  frequency  of  the  synchronizing 
torque  loop 

the  per  unit  load  carried  by  the  individual 
generator 


Xd  +  XT  +  XL 

the  synchronous  reactance  of  the  generators 
the  transformer  reactance 
the  line  reactance 


the  synchronizing  torque  angle 

the  permanent  speed  droop 

the  frequency  datum  angle 

the  angles  of  generators  #1  and  #2 
respectively ,  relative  to  ground 

the  angle  associated  with  the  mutual 
impedance  between  the  generators 
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APPENDIX  B 


VALUE  OF  THE  MODEL  PARAMETERS 

Section 

II-l  Generator  and  Voltage  Regulator  Parameters 

Xad 

lol  per  unit  Base  MVA  =  40.0 

Xad  saturated  0.66  per  unit  Base  KV  =  138.0 


xal 

0.15  per  unit 

xf  Id 

0.15  per  unit 

XKLD 

0.04  per  unit 

Rfd 

0.0017  per  unit 

rkd 

0.05  per  unit 

rkq 

0.05  per  unit 

T'd 

1.5  seconds 

H 

1.5  KW-seconds/KVA 

Machine 

Rating  40.0  MVA 

K 

1.00  per  unit 

K2 

5.00  per  unit 

K3 

0.2  per  unit 

T 

0.1  seconds 

Section 

II-2  Governor  Parameters 

K 

5.00  per  unit 

s  T 

0.214  per  unit 

Tr 

3.00  seconds 

Pd 

3.92  per  unit 

Variable 
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Section  II-3 

Hydro-System  Parameters 

T 

w 

Variable 

Section  II-4 

Load  Parameters 

8L 

2.0,/  36°  per  unit  for  light  load 
loO  zl  1°  per  unit  for  heavy  load 

XT 

OolO  per  unit 

XL 

104  ohm's  or  0.218  per  unit  for  150 
miles  of  line.  It  is  assumed  that 
there  are  75  miles  of  line,  operating 
@  138  KV,  between  each  generator  and 
the  loado  Further  it  is  assumed  that 
the  line  is  capacitively  compensated 

X  =  8/3  R 

R 

38.9  ohm's  or  0.084  per  unit  for  150 
miles  of  line 

Section  II-5 

Controller  Model  Parameters 

Quadrant  Time 
Constant 

165  seconds 

Backlash 

1/50  seconds 

Dead  Zone 

Variable 

Effective  Speed 
Droop  Signal 

Dashpot  Frequency 
Signal 

Percentage  Load 
Setting 


Sampling  Duration 
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pei  V 

■ 


' 


c/o 

^A) 

C/O 

10 

c 

r< 

cr 

0 

c 

<J 

7> 

Ul 

o 

V0 

<r 

LU 

cr 

LU 

cr 

iLl 

5- 

LU 

UQ 

\A 

O 

o 


</0 

£ 

V) 

N 

<T 

LU 


lu 


C/O 

</) 

o 

u 


C/O 

c 

v> 

c< 

o- 


?M<X>9  W<X> 


Signal 


